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The gravitational couplings of matter are studied in the presence of Lorentz and CPT violation. 
At leading order in the coefficients for Lorentz violation, the relativistic quantum hamiltonian is 
derived from the gravitationally coupled minimal Standard-Model Extension. For spin-independent 
effects, the nonrelativistic quantum hamiltonian and the classical dynamics for test and source bodies 
are obtained. A systematic perturbative method is developed to treat small metric and coefficient 
fluctuations about a Lorentz- violating and Minkowski background. The post-newtonian metric and 
the trajectory of a test body freely falling under gravity in the presence of Lorentz violation are 
established. An illustrative example is presented for a bumblebee model. The general methodology 
is used to identify observable signals of Lorentz and CPT violation in a variety of gravitational 
experiments and observations, including gravimeter measurements, laboratory and satellite tests of 
the weak equivalence principle, antimatter studies, solar-system observations, and investigations of 
the gravitational properties of light. Numerous sensitivities to coefficients for Lorentz violation can 
be achieved in existing or near-future experiments at the level of parts in 10'^ down to parts in 
10^^. Certain coefficients are uniquely detectable in gravitational searches and remain unmeasured 
to date. 
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I. INTRODUCTION 

General Relativity (GR) is known to provide an accu- 
rate description of classical gravitational phenomena over 
a wide range of distance scales. A foundational property 
of the gravitational couplings of matter in GR is local 
Lorentz invariance in freely falling frames. The realiza- 
tion that a consistent theory of quantum gravity at the 
Planck scale mp ~ 10^^ GeV could induce tiny man- 
ifestations of Lorentz violation at observable scales [1] 
has revived interest in studies of Lorentz symmetry, with 
numerous sensitive searches for Lorentz violation being 
undertaken in recent years [2] . 

Gravitational signals of Lorentz violation are more 
challenging to study than ones in Minkowski spacetime 
for several reasons, including the comparative weakness 
of gravity at the microscopic level and the impossibility 
of screening gravitational effects on macroscopic scales. 
Both for purely gravitational interactions and for matter- 
gravity couplings, Lorentz violations can be classified and 
enumerated in effective field theory [3]. Several searches 
for purely gravitational Lorentz violations in this con- 
text have recently been performed [4-6] using a treat- 
ment in the post-newtonian regime [7]. These results 
enlarge and complement the impressive breadth of tests 
of GR performed in the context of the parametrized post- 
newtonian (PPN) formalism [8]. 

Although the coupling between matter and gravity has 
historically been a primary source of insights into the 
properties of the gravitational field, a general study of 
matter-gravity couplings allowing for Lorentz violation 
in the context of effective field theory has been lacking 
to date. In this work, we address this gap in the literature 
and investigate the prospects for searches for Lorentz vi- 
olation involving matter-gravity couplings. Our goal is 
to elucidate both theoretical and experimental aspects 
of the subject. We seek a post-newtonian expansion for 



the equation for the trajectory of a test body moving un- 
der gravity in the presence of Lorentz violation, allowing 
also for Lorentz-violating effects from the composition of 
the test and source bodies and for effects from possible 
additional long-range modes associated with Lorentz vi- 
olation. We also seek to explore the implications of our 
analysis in a wide variety of experimental and observa- 
tional scenarios, identifying prospective measurable sig- 
nals and thereby enabling more complete searches using 
matter-gravity couplings. 

Despite the current lack of a satisfactory quantum 
theory of gravity, established gravitational and particle 
phenomena at accessible energy scales can successfully 
be analyzed using the field-theoretic combination of GR 
and the Standard Model (SM). This combination there- 
fore serves as a suitable starting point for a comprehen- 
sive effective field theory describing observable signals 
of Planck-scale Lorentz and CPT violation in gravity 
and particle physics [9]. The present paper adopts this 
general framework, known as the gravitational Standard- 
Model Extension (SME) [3], to analyze Lorentz violation 
in matter-gravity couplings. Each term violating Lorentz 
symmetry in the SME Lagrange density is a scalar den- 
sity under observer general coordinate transformations 
and consists of a Lorentz-violating operator multiplied 
by a controlling coefficient for Lorentz violation. Under 
mild assumptions, GPT violation in effective field theory 
comes with Lorentz violation [10], so the SME also de- 
scribes general breaking of CPT symmetry. This feature 
plays a crucial role for certain signals of Lorentz violation 
in what follows. 

In this work, our focus is on gravitational Lorentz vi- 
olation in matter-gravity couplings, both with and with- 
out CPT violation. These couplings introduce operator 
structures offering sensitivity to coefficients for Lorentz 
violation that are intrinsically unobservable in Minkowski 
spacetime. In fact, comparatively large gravitational 
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Lorentz violation in nature could have remained unde- 
tected in searches to date because gravity can provide 
a countershading effect [11], so this line of investigation 
has a definite discovery potential. Several searches for 
gravitational Lorentz violation have led to constraints on 
SME coefficients for Lorentz violation with sensitivities 
down to parts in 10^ [4-6, 12], and additional constraints 
can be inferred by reanalysis of data from equivalence- 
principle tests [11]. 

The nature of the Lorentz violation plays a crucial role 
in determining the physics of matter-gravity couplings. 
In Riemann geometry, externally prescribing the coef- 
ficients for Lorentz violation as fixed background con- 
figurations is generically incompatible with the Bianchi 
identities and hence problematic [3]. This issue can be 
avoided via spontaneous Lorentz breaking [1] , in which a 
potential term drives the dynamical development of one 
or more nonzero vacuum values for a tensor field. This 
mechanism implies the underlying Lagrange density is 
Lorentz invariant, so the coefiicients for Lorentz violation 
are expressed in terms of vacuum values and can there- 
fore serve as dynamically consistent backgrounds satis- 
fying the Bianchi identities. The presence of a potential 
driving spontaneous Lorentz violation implies the emer- 
gence of massless Nambu-Goldstone (NG) modes [13] as- 
sociated with field fluctuations along the broken Lorentz 
generators [14] . If the potential is smooth, massive modes 
can also appear [15]. Some features of the NG and mas- 
sive modes are generic, while others are specific to details 
of the model being considered. In any case, the nature of 
these modes plays a key role in determining the physical 
implications of spontaneous Lorentz violation. 

For the purposes of the present work, the presence of 
NG modes is particularly crucial because they can cou- 
ple to matter and can transmit a long-range force. A 
careful treatment of these modes is therefore a prereq- 
uisite for studies of Lorentz violation in matter-gravity 
couplings. In what follows, we develop a methodology to 
extract the dominant Lorentz- violating effects in matter- 
gravity couplings irrespective of the details of the un- 
derlying model for spontaneous Lorentz violation. In ef- 
fect, the NG modes are treated via a perturbation scheme 
that takes advantage of symmetry properties of the un- 
derlying Lagrange density to eliminate them in favor of 
gravitational fluctuations and background coefficients for 
Lorentz violation. This treatment allows leading Lorentz- 
violating effects from a large class of plausible models to 
be handled in a single analysis. 

The portion of this paper developing theoretical issues 
spans Sees. II-V. It begins in Sec. II with a review of the 
SME framework. We present the field-theoretic action, 
describe its linearization, and discuss observability issues 
for the coefficients for Lorentz violation. We also describe 
the two notions of perturbative order used in the subse- 
quent analysis, one involving Lorentz and gravitational 
fluctuations and the other based on a post-newtonian ex- 
pansion. Section III concerns the relativistic and nonrel- 
ativistic quantum mechanics arising from the quantum 



field theory. One technical issue is extracting a mean- 
ingful quantum theory in the presence of gravitational 
fiuctuations. We resolve this issue via a judicious field 
redefinition, which yields a hamiltonian that is hermi- 
tian with respect to the usual scalar product for wave 
functions and that reduces correctly to known limiting 
cases. We construct the relativistic quantum hamilto- 
nian at leading order in Lorentz violation and gravity 
fluctuations. For the spin-independent terms, we per- 
form a Foldy-Wouthuysen transformation to obtain the 
nonrelativistic hamiltonian. Section IV treats the classi- 
cal dynamics corresponding to the quantum theory. The 
point particle action is presented and the structure of test 
and source bodies is discussed. The equations of motion 
for a test particle and the modified Einstein equations are 
derived. We describe the methodology for handling coef- 
ficient and metric fluctuations. Combining the results de- 
termines the Lorentz- violating trajectory of a test body. 
The results are illustrated in Sec. V in the context of a 
special class of models of spontaneous Lorentz violation 
known as bumblebee models. 

The remaining research sections of the paper. Sees. 
VI-XI, concern implications of our theoretical analysis 
for experiments and observations. Section VI contains 
basic facts concerning frame choices and outlines the 
canonical Sun-centered frame used for reporting mea- 
surements. We also consider the sensitivities to coef- 
ficients for Lorentz violation that can be attained in 
practical situations. Section VII treats laboratory tests 
near the surface of the Earth using neutral bulk mat- 
ter, neutral atoms, or neutrons. The theoretical descrip- 
tion of these tests is presented to third post-newtonian 
order, and some generic features of the test-body mo- 
tion are discussed. A wide variety of gravimeter and 
equivalence-principle tests is analyzed for sensitivities to 
coefficients for Lorentz violation that are presently un- 
constrained. Satellite-based searches for Lorentz viola- 
tion using equivalence- principle experiments are studied 
in Sec. VIII. A generic situation is analyzed, and the 
results are applied to major proposed satellite tests. Sec- 
tion IX treats gravitational searches using charged par- 
ticles, antihydrogcn, and particles from the second and 
third generations of the SM. Estimated sensitivities in 
future tests are provided, and illustrative toy models for 
antihydrogen studies are discussed. Searches for Lorentz 
violation using solar-system observations are described 
in Sec. X. We consider measurements of coefficients for 
Lorentz violation accessible via lunar and satellite rang- 
ing and via studies of perihelion precession. Section XI 
addresses various tests involving the effects of gravita- 
tional Lorentz violation on the properties of light. We 
analyze the Shapiro time delay, the gravitational Doppler 
shift, and the gravitational redshift, and we consider 
the implications for a variety of existing and proposed 
searches of these types. Finally, in Sec. XII we summa- 
rize the paper and tabulate the various estimated actual 
and attainable sensitivities to coefficients for Lorentz and 
CPT violation obtained in the body of this work. 
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Throughout the paper, we fohow the conventions of 
Rcfs. [3] and [7]. In particular, the Minkowski metric is 
diagonal with entries (— 1, + !,+!, +1). Greek indices are 
used for spacetime coordinates, while Latin indices are 
used for local Lorentz coordinates. Appendix A of Ref. 
[3] provides a summary of most other conventions. Note 
that parentheses surrounding index pairs in the present 
work denote symmetrization with a factor of one-half. 



II. FRAMEWORK 

The focus of this work is the study of relativity viola- 
tions in realistic matter-gravity interactions. The basic 
field theory of relevance concerns a single fermion field ip 
coupled to dynamical gravity and incorporating Lorentz 
and CPT violation. In this section, we summarize the 
action for the model, describe the linearization proce- 
dure, discuss conditions for the observability of effects, 
and present the perturbation scheme developed for the 
analysis to follow. 



A. Action 

The theory of interest is a special case of the gravi- 
tationally coupled SME [3]. The action can be written 

as 

S = Sg + S^ + S'. (1) 

The first term in this expression is the action Sq contain- 
ing the dynamics of the gravitational field, including any 
coefficients for Lorentz violation in that sector. The geo- 
metric framework is a Ricmann-Cartan spacetime, which 
allows both the Riemann curvature tensor R'^Xfiiy '^''^'^ 
torsion tensor T^^^. To incorporate fermion-gravity in- 
teractions, the vierbcin formalism [16] is adopted, with 
the vierbein e'^^ and the spin connection w^"*' taken as 
the fundamental gravitational objects. In the limit of 
zero torsion and Lorentz invariance, Sq reduces to the 
Einstein-Hilbert action of General Relativity, 

Sg ^ f d'^x (eR-2eA), (2) 
2k J 

where k = SttGn, e is the vierbein determinant, and A 
is the cosmological constant. 

The second term in Eq. (1) is the action for the 
fermion sector of the SME. In this work, we limit at- 
tention to terms in this sector with no more than one 
derivative, which is the gravitationally coupled analogue 
of the minimal SME in Minkowski spacetime. In this 
limit, the action 5^ for a single Dirac fermion ip of mass 
m can be written as 

S^= [ d'^xi^iee^JfT" D^i)- ei^Mip). (3) 



In the present context, the action of the covariant deriva- 
tive Dfj, on ip is 

D^V = + izc^/VabV- (4) 
It is convenient to introduce the symbol 

(?:D^) = a^v^-i*^/V^„6 (5) 

for the action of the covariant derivative on the Dirac- 
conjugate field ip. The action (3) contains the covariant 
derivative in a combination defined by 

xr« D^i, = xr"£>^v - (xS^)r"V. (6) 

The symbols F" and M appearing in the action (3) are 

defined by 

-e^.e"" - - ^5A;,.e'^«e\e'^,a''= (7) 

and 

M = m + a^e\^'^ + 6^e'',757'^ + \ii ^.e\e\o''\ (8) 

The first term of Eq. (7) leads to the usual Lorentz- 
invariant kinetic term for the Dirac field, while the first 
term of Eq. (8) leads to a Lorentz-invariant mass. A term 
of the form im^j^ could also appear in M, but here we 
suppose it is absorbed into m via a chiral field redefini- 
tion. The coefficient fields for Lorentz violation a^, 6^, 
c^iy, d^^, e^, ffj., gx^iv, H^^ typically vary with space- 
time position. The coefficient field is antisymmet- 
ric, while ^A/ni/ is antisymmetric in A/x. Note the use of 
an uppercase letter for Hfj,,^, which avoids confusion with 
the metric fluctuation /i^,^. The CPT-odd operators for 
Lorentz violation are associated with the coefficient fields 
a^L, ''m' and gxf,^. 

The form of the action (3) implies the torsion T^^^ 
enters the fermion action only via minimal coupling. This 
coupling has the same form as that of the coefficient field 
b/j,, so the effects of minimal torsion can be incorporated 
into a matter-sector analysis by replacing with the 
effective coefficient field 

(6eff)^ = b^ + iT«^Te„^^^. (9) 

Note that nonminimal torsion couplings can be incorpo- 
rated into the more general coefficient fields appearing 
in the full SME. Nonminimal torsion couplings and their 
experimental constraints are discussed in Ref. [17]. 

The final term in Eq. (1) is the action S' containing the 
dynamics associated with the coefficient fields for Lorentz 
violation. Addressing possible contributions from this 
sector is the subject of Sec. IV C. 

B. Linearization 

For the purposes of this work, it suffices to consider 
weak gravitational fields in a Minkowski-spacetime back- 
ground. Under these circumstances, the Latin local in- 
dices can be replaced with Greek spacetime indices, so 
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the weak-field forms of the metric, vierbein, and spin 
connection can be written as 

+5 {Tx^'^ + T^Xi^ - Tux^d ■ (10) 

The quantities Xtiv contain the six Lorentz degrees of 
freedom in the vierbein. 

The coefficient fields for Lorentz violation are expected 
to acquire vacuum values through spontaneous Lorentz 
breaking. An arbitrary coefficient field tx^^... can there- 
fore be expanded about its vacuum value txixv..., 

_ ^ 

tXuv... = tXijiv...+tx^v..., (11) 

where the fluctuation tx^j^u... includes massless NG modes 
and massive modes [14, 15]. The vacuum value txuv... = 
{txtiu...} is called the coefficient for Lorentz violation. One 
can instead choose to expand the contravariant coefficient 
field i^"''-, 

where f'^'''' " is related to txtiv... by raising with rj'^" . The 
reader is cautioned that the relation between tx^iv... and 
the index-lowered version txnu... of t^i^"--- involves terms 
containing contractions of ixfiv... with h^j^v This paper 
uses the expansion (11) and gives expressions in terms of 

^Xixv...- 

To provide a smooth match between our analysis and 
previous work on the matter sector of the SME in 
Minkowski spacetime and on the gravitational sector in 
asymptotically Minkowski spacetime, we make two as- 
sumptions about the coefficients for Lorentz violation. 
First, we assume they are constant in asymptotically in- 
ertial cartesian coordinates, 

^Jx^.u... = 0. (13) 

This preserves translation invariance and hence 
energy-momentum conservation in the asymptoti- 
cally Minkowski regime. It also ensures that our barred 
coefficients correspond to the usual coefficients for 
Lorentz and CPT violation investigated in the minimal 
SME in Minkowski spacetime [18]. Second, we assume 
that the vacuum values txfiv... are sufficiently small to be 
treated perturbatively. This is standard and plausible, 
since any Lorentz violation in nature is expected to 
be small. These two assumptions suffice for most 
of the analysis that follows. To obtain the leading 
Lorentz-violating corrections to hf^i, without specifying 
a dynamical model for the coefficient fields for Lorentz 
violation, one further assumption is required, which is 
presented in Sec. IV C. 



C. Observability 

A given coefficient for Lorentz violation can lead to 
observable effects only if it cannot be eliminated from 
the Lagrange density via field redefinitions or coordinate 
choices [3, 18-25]. In this subsection, we outline some 
implications of this fact relevant to the present work. 

1. Field redefinitions 

One result of key interest here is that matter-gravity 
couplings can obstruct the removal of some coefficients 
that are unphysical in the Minkowski-spacetime limit. 
For example, in the; single-fermion theory in Minkowski 
spacetime, the coefficient = for Lorentz and CPT 
violation in Eq. (8) is unobservable because it can be 
eliminated by the spinor redefinition 

V'(a;) exp[i/(a;)]V'(a;) (14) 

with f{x) = a^x^. However, in Riemann or Riemann- 
Cartan spacetimes we have = af^ + a^, so this redefi- 
nition typically leaves the four components of the fluctu- 
ation Sjii in the theory. Instead, the redefinition (14) with 
an appropriate f{x) can be used to move one component 
of the coefficient field into the other three, unless is 
constant or the total derivative of a scalar [3] . Note that 
in the presence of gravity this freedom may be insufficient 
to eliminate any components of the coefficient because 
the components of the fluctuation can depend on all 
four components of through the equations of motion. 
This line of reasoning shows that gravitational couplings 
provide a unique sensitivity to the coefficient a^, which 
can be exploited in various experiments [11]. 

Another type of field redefinition can be written in the 
generic form [3] 

i){x)^[l + v{x)-r]i,{x), (15) 

where F represents one of 7°, 757", cr"'' and v{x) is a 
complex function carrying the appropriate local Lorentz 
indices. This can be viewed as a position-dependent com- 
ponent mixing in spinor space. Field redefinitions of this 
type can be used to demonstrate the leading-order equiv- 
alence of observable physical effects due to certain coeffi- 
cients for Lorentz violation. An example relevant in the 
present context is a redefinition involving a real vector 
Va{x). Together with assumption (13), this redefinition 
can be used to show that at leading order in Lorentz vi- 
olation the coefficients and always appear in the 
combination 

{acs)ti = a^- me^, (16) 

up to derivatives of fluctuations. Combining this result 
with the above discussion of the redefinition (14) shows 
that observables involving gravitational couplings offer 
the prospect of measuring (oeff)/*. 
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Related ideas can be used to simplify the weak-field 
limit of the theory (3). In particular, the antisymmet- 
ric part x^ll' of the vierbein can be removed everywhere 
except for possible contributions to fluctuations of the 
coefficient fields, by applying the field redefinition 

« (1 - itx^^.a'^'' - hx^uXafi'^^^'j'^^) ^{x). 

(17) 

Note that this redefinition takes the form of a Lorentz 
transformation on tf) but that the other fields in the La- 
grange density remain unaffected. Note also that in the 
absence of Lorentz violation Xtiv can be removed entirely, 
a fact compatible with the interpretation in Ref. [14] of 
the role of Xfiv in Lorentz- violating theories. In the re- 
mainder of this work, we assume the redefinition (17) has 
been performed on the Lagrange density, so that quan- 
tities such as r", M, and are understood to acquire 
no contributions from Xiiv except possibly through the 
fluctuations of the coefiicient flelds for Lorentz violation. 



2. Coordinate choices 

The observability of certain combinations of coeffi- 
cients for Lorentz violation is also affected by the freedom 
to make coordinate choices. Intuitively, the key point is 
that any one sector of the SME can be used to define the 
scales of the; four coordinates, to establish the meaning 
of isotropy, and to set the synchronization scheme. The 
freedom therefore exists to choose the sector in which the 
effective background spacetime metric takes the form of 
the Minkowski metric r]^,^. This implies that in any ex- 
perimental configuration there are always 10 unobserv- 
able combinations of coefficients for Lorentz violation. 

As an illustration, consider the SME restricted to the 
single- fermion and photon sectors [18, 22]. In the fcrmion 
sector, the 10 relevant coefficient components are the vac- 
uum values c^^ of the coefficient fields c^^ in Eq. (7) be- 
cause these coefficients enter S-^, in the same way as the 
metric. In the photon sector, the SME Lagrange density 
contains a term 

£photonD-ie(fci.),A;..F"^-F''", (18) 

and the 10 relevant coefficient components can be shown 
to be the trace (fci?)"^„^. At leading order, the coordi- 
nate transformation 

^ x''' =x^- l^T^^^x" (19) 

redefines the background metric to take the form r]^v in 
the photon sector. The effective metric in the fermion 
sector is then also changed, with the observable coeffi- 
cient combination becoming c^j, + (fcF)"^ai//2. The or- 
thogonal combination 2c^y — {kF)"fj_au is thus unobserv- 
able in any experiments involving only these two sectors. 



Alternatively, one could perform the coordinate transfor- 
mation 

_^ 2,m' =x^' + c^'^x", (20) 

which instead redefines the fermion-scctor background 
metric to be rjni, and at leading order produces the ef- 
fective photon-sector coefficient 2c^y -|- {kF)°'^av This 
coordinate choice is equally valid for analysis, and as be- 
fore the orthogonal combination 2c^^ — (fcF)"^a,y is un- 
observable. 

Similar results apply for experimental searches for 
Lorentz violation involving comparisons of different 
fcrmion species. Labeling the species by w, each has 
a coefficient c™^. Then, for example, the effective metric 
for any one species X can be reduced to 77^^ by a co- 
ordinate transformation with of the form (20). The 
resulting effective coefficients for the remaining species 
involve the differences c^^, — c^i,, which in this coordinate 
scheme become the relevant observable combinations of 
coefficients. 

In the gravity sector, the situation is more involved 

because a geometrically consistent treatment of Lorentz 
violation generically requires the incorporation of effects 
from the NG modes to ensure the Bianchi identities are 
satisfied [3]. It turns out that the 10 relevant coefficient 
components in the gravity sector are the vacuum values 
s^^ of the coefficient fields s^i, in the gravity-sector SME 
Lagrange density 

.^gravity ^ "77. 7-* 65^^ ^uiy (^l) 

IottGat 

We find that at leading order the transformation (19) 
generates an accompanying shift s^v s^v — (fcF)"^^,^, 
while the transformation (20) produces the shift s^v 

One way to obtain these results is to consider the 
leading-order effect of a metric shift on the equations of 
motion, and then to match to the known results [7] for 
observable effects in a post-newtonian expansion. Con- 
sider, for example, the restriction of the SME to the 
Einstein-Hilbert action and the single-fermion action 
with nonzero c^^ only. At leading order, the coordinate 
transformation (20) removes from the fermion action 
at the cost of introducing a metric shift gi^^ — )• Qfj^v — 20^^ 
in the Einstein-Hilbert term. The resulting equations 
of motion involve the Einstein tensor Gfj,^{g — 2c) with 
shifted argument, which can be written in terms of the 
Einstein tensor G^vig) for the original metric and an ef- 
fective energy-momentum tensor ^^j, . We find 

Gij,u{g~2c) = Gf,^{g) - (f>%, 

+r^Ro.t.uf3). (22) 

The trace-reversed form of matches the post- 
newtonian term <l?^j^ arising from Eq. (21) and given 
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explicitly as Eq. (24) of Ref. [7] with the combination 
Sfiv replaced by 2c^y. The transformation (20) there- 
fore produces the shift s^i/ — >■ Sy^jy + 2c^^, as claimed. 
A similar line of reasoning verifies the claimed shift 
Spi/ — > s^i/ — {kF)°'u,av for the coordinate transformation 
(19). 

To keep expressions compact throughout this work, we 
choose to work with coordinates satisfying 

(M%a. = 0. (23) 

To obtain results valid for arbitrary coordinate choices, 
the following substitutions can be applied throughout: 

r'^ r"- + KTT^^" 

^ s^,.-(MV-^- (24) 

We emphasize that all coordinate choices are equivalent 
for theoretical work or for data analysis, with the choice 
(23) adopted here being purely one of convenience. 

D. Perturbation scheme 

In this work, we are interested in experimental searches 
for Lorentz violation involving gravitational effects on 
matter. Many of these searches involve test particles 
moving in background solutions to the equations of mo- 
tion for gravity and for the coefficient fields. Since the 
gravitational fields involved are weak and since no com- 
pelling evidence for Lorentz violation exists to date, any 
effects are expected to be small. We can therefore focus 
attention on perturbative modifications to the behavior 
of test particles. This subsection describes the scheme we 
use to track perturbative orders in the construction of the 
relativistic quantum hamiltonian and in the subsequent 
developments. 

Perturbative effects on physical observables can arise 
through modifications to the background coefficient fields 
t\fiv... for Lorentz violation and to the background gravi- 
tational field 3^1^, or directly through modifications to the 
equation of motion for the test particle. The analysis of 
these effects is simplified by introducing an appropriate 
notion of perturbative order. Several ordering schemes 
arc possible. In this work, we adopt a scheme that tracks 
the orders in the coefficients for Lorentz violation txfj,v... 
and in the metric fluctuation /i^^. The overall pertur- 
bative order of a given term in an equation is denoted 
as 0(m,n), where m represents the order in t\nv... and 
n the order in h^^^. Within this scheme, the fluctua- 
tions t\ij,v... of the coefiicient fields for Lorentz violation 
are viewed as secondary quantities that are determined 
via their equations of motion in terms of the coefficients 
for Lorentz violation and the gravitational field. There is 
also a subsidiary notion of order associated with the usual 
post-newtonian expansion of /i^j, itself. We denote a pth- 
order term in this latter expansion as PNO(p). However, 
performing an explicit post-newtonian expansion at the 



initial stage would complicate the ensuing analysis, so in 
what follows we often write results in terms of h^i, while 
commenting as needed on the post-newtonian counting. 

To preserve a reasonable scope, this work focuses on 
dominant perturbative effects involving both Lorentz vi- 
olation and gravity. We next discuss the relevant pertur- 
bative orders required to achieve this goal. 

Consider first contributions from the fluctuation ^Api/... 
of the coefficient fields. The detailed structure of tx^iv... 
depends on the nature of the action S" in Eq. (1). In the 
scheme adopted here, tx^i,,,, can be viewed as a series in 
txixv... and h^v of the form 

t - ,7(0,1) ^1,0) , , /2^^ 

For spontaneous Lorentz breaking, txixv... includes mas- 
sive modes and massless NG modes [15, 26]. In this work, 
we suppose the massive modes either arc frozen or have 
negligible degree of excitation. Incorporating their possi- 
ble effects into the analysis of matter-gravity phenomena 
would be of potential interest but lies beyond our present 
scope. In contrast, the massless NG modes play a key 
role in what follows. Their fate can include identification 
with the photon in Einstein-Maxwell theory [14], with the 
graviton in GR [27, 28], or with a new force [11, 29, 30], 
or they can be absorbed in the torsion via the Lorentz- 
Higgs effect [14]. In some models the NG modes can be 
interpreted as composite photons [31] or gravitons [32]. 
We consider below the perturbative orders required for 
the various possibilities. 

Suppose first the NG modes correspond to photons, or 
more generally to a known force field other than grav- 
ity. The term ^^^J then contains conventional Lorentz- 
invariant terms describing this field in Minkowski space- 
time, while 'V'y^l^J contains conventional leading-order 
gravitational interactions with the field. Effects from 
both these terms are therefore part of the conventional 
description of the force. The term 'i^'^^^J describes possi- 
ble Lorentz violations in Minkowski spacetime involving 
the known field, many of which are tightly constrained by 
experiments [2] . For the purposes of this work, which fo- 
cuses on Lorentz violation involving gravity, we can take 
this term as experimentally negligible. The dominant 
term of interest is therefore V-^^^J , which lies at 0(1,1). 

If the NG modes correspond to gravitons as, for ex- 
ample, in the cardinal model [27], then the expansion 
(25) contains no terms at O(m,0). The term 'i^^^J cor- 
responds to the gravitational fiuctuations h^v, and its 
effects are part of the conventional description of gravity. 
The dominant term of interest is therefore again 't^^^^J ■ 

If instead the NG modes correspond to a presently un- 
observed force field, then and ^xlja! describe un- 
observed Lorentz-invariant effects in Minkowski space- 
time and in leading-order gravitational couplings. These 
modes must therefore be eliminated from the analy- 
sis prior to interpretation of observations, via solving 
the equations of motion or otherwise. In what follows. 
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we suppose this elimination has been performed where 
needed. The term describes possible Lorentz vi- 

olation in Minkowski-spacetime involving the unknown 
field. For present purposes, wc take this term to be exper- 
imentally negligible, although in principle it might offer a 
novel way to access certain types of presently unobserved 
interactions at exceptional sensitivities. The remaining 

term t\^J displayed in Eq. (25) describes the dominant 
Lorentz-violating gravitational effects involving the un- 
known interaction. As discussed in Sec. II CI, certain 
Lorentz-violating effects are observable only in the pres- 
ence of gravity, and so under suitable circumstances ob- 
servable experimental signals from t^^^^^ could arise [11] 
despite the tight existing experimental constraints [33] 
on the direct observation of additional interactions due 
to the lower-order terms t^^^f] and Vf^^^ . In scenarios 
with an unobserved force, we must therefore also allow 
for 0(1,1) effects involving i^^^^^ • 

The remaining possibility is that NG modes are ab- 
sorbed into the torsion. To handle this case, note that 
the matter-sector role of minimal torsion can be treated 
in parallel with the coefficient field &^ for Lorentz vio- 
lation according to Eq. (9). Existing experimental con- 
straints on minimal torsion components are tight, lying 
below 10-2^-10-31 GeV [17, 34, 35], so eflfects involving 
minimal torsion can be treated in the same way as those 
corresponding to a presently unobserved force field. We 
can therefore limit attention to 0(1,1) effects as before. 
In principle, any nonminimal torsion components can be 
treated in a similar fashion because they play a role anal- 
ogous to coefficient fields in the nonminimal SME, but 
effects of this type lie beyond our present scope. 

Consider next the metric fluctuation h^^. For appli- 
cations to gravitational tests with matter, can be 
treated as a background field obtained by solving the 
appropriate equation of motion, which is the modified 
Einstein equation in the presence of Lorentz violation. It 
can therefore be viewed as the sum of a Lorentz-invariant 
piece /i^y^'' with a series of corrections of increasing per- 
turbative order in t\^i,,,,, 

V = + + h^^' + ■■■■ (26) 

When we specify the perturbative order of an expression 
containing /i^^ in what follows, it is understood that the 
correct terms from the above series are included. 

For a given expression, establishing the relevant per- 
turbative order for our analysis typically involves a com- 
bination of experimental restrictions and theoretical con- 
siderations. As an illustration, we outline here the rea- 
soning establishing the appropriate perturbative orders 
in the construction of the relativistic quantum hamilto- 
nian. First, note that terms quadratic in /i^j^ involve 
PN0(4) and higher. Since the sensitivity of current lab- 
oratory and solar-system tests lies at the PN0(4) level, 
we must keep these terms but can discard terms cu- 
bic in h^j^u and ones involving the product of coefficients 
for Lorentz violation with terms quadratic in h^v The 



Lorentz-invariant part of the hamiltonian can therefore 
be truncated at 0(0,2), while the Lorentz-violating part 
can be truncated at 0(1,1). To maintain consistent post- 
newtonian counting, the 0(0,2) terms must be limited to 
PN0(5), while the 0(1,1) terms are limited to PN0(3). 
Next, note that for laboratory and solar-system tests 
the variations in over the experimental scale L are 
small compared to h^i,, \dahfii,\ <C For exam- 

ple, the typical value of the gravitational acceleration on 
the surface of the Earth is 5 ~ IQ-'^^ GeV, which is 
tiny compared to the ratio of the gravitational poten- 
tial — IQ-^ and the size of a typical laboratory 
experiment L ~ IQ^^ GcV^i. Terms in the relativis- 
tic hamiltonian proportional to derivatives of /i^j/ can 
therefore be limited to 0(0,1). Finally, note that prod- 
ucts of Lorentz-violating terms lead to higher-order ef- 
fects with operator structures matching ones already ap- 
pearing in the fermion sector of the Minkowski-spacetime 
SME. These are already accessible in nongravitational ex- 
periments. It therefore suffices for our purpose to restrict 
attention to terms at leading order in Lorentz violation. 
To summarize, the construction of the perturbative rela- 
tivistic quantum-mechanical hamiltonian can be limited 
to terms at perturbative orders 0(0,1), 0(1,0), 0(1,1), 
and 0(0,2), except for terms involving derivatives of the 
gravitational fields, which can be limited to 0(0,1). 



III. QUANTUM THEORY 

This section studies the quantum mechanics associated 
with the fermion action in Eq. (3). We begin in Sec. 
Ill A by addressing the issue of the unconventional time 
dependence arising from the Dirac equation derived from 
S^. The relativistic quantum- mechanical hamiltonian H 
is then obtained in Sec. IIIB, and the relevant parts of 
the nonrelativistic limit JInr are extracted in Sec. IIIC. 



A. Time dependence 

In the weak-field limit, the Lagrange density for 
the action (3) takes the schematic form 

£ - i2[^(7^ + C^)a^V-(5^V^)(7^ + C^^)^]-?^V-', (27) 

where and D represent spacetime-dependent opera- 
tors without derivatives acting on ij). These operators 
satisfy the conditions 

(70(7'')t = (70£))t = (28) 

and C is perturbative. 

The Eulcr-Lagrange equations obtained from Eq. (27) 
yield a Dirac equation with unconventional time depen- 
dence, 

i{f + C°)do^p = [-i(7^ + C')dj - \id^C^ + D],p. (29) 
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This equation differs from the standard Dirac form by 
the presence of C'^ , which impedes the interpretation of 
the operator acting on (/; on the right-hand side as the 
hamiltonian. In this subsection, two approaches address- 
ing this issue at first order in C are discussed. 

1. Field-redefinition method 

One method for constructing the hamiltonian has been 
developed in the context of the SME in Minkowski space- 
time [36- 38]. It uses an appropriate field redefinition 
at the level of the action to ensure the Dirac equation 
emerges with conventional time dependence. In typical 
applications, the field redefinition is defined perturba- 
tively at the desired order in C^. 

For present purposes, it suffices to work at first order 
in C^. The appropriate field redefinition is 

^ = Ax, A=l-i70c°. (30) 

The resulting hamiltonian can be written as 

= H^°'> + Hll\ (31) 

where H^^^ is the hamiltonian in the absence of C^, given 

by 

= -i7V5j+7°£'- (32) 

The correction H^^ is first order in C and takes the 
form 

41) = _i/(c^- - 1CO7V + 57Vc°)a,- 

-'^j°{C°j°D + D^°C'). (33) 

The subscript x serves as a reminder that the operator 
acts on the spinor x- Note that the hamiltonian is 
hermitian with respect to the usual scalar product in flat 

space, 

(Xi,X2)f = J d^xx{x2- (34) 

This implies, for example, that energies can be calculated 
in the usual way. 

Some physical insight into the field-redefinition method 
is obtained by noting that the combination 7'' -|- C takes 
the generic form 

^^^ + C^ = E^^-y", (35) 

where E'^^ can be interpreted as an effective inverse vier- 
bein. It reduces to the conventional inverse vierbein e''„ 
in the purely gravitational case but includes coefficients 
for Lorentz violation when Lorentz symmetry is broken. 
This suggests that the field-redefinition method can be 
interpreted as transforming the problematic situation of 
a fermion on an effective manifold with vierbein com- 
ponents into the physically equivalent but tractable 
theory of a different fermion field on a manifold with vier- 
bein components E'^^ = S'^^, in which the hamiltonian is 
hermitian with respect to a conventional scalar product. 



2. Parker method 

Another method has been presented by Parker [39] in 
the context of field theory in curved spacetime. It in- 
volves multiplying the Dirac equation by a suitable fac- 
tor that removes the unconventional time dependence to 
the desired order in C^. The resulting hamiltonian is 
hermitian with respect to a modified scalar product. 

Applying this method at first order in requires left- 
muhiplying both sides of Eq. (29) with 7°(1-C°7°). The 
ensuing hamiltonian can be written as 

i?^=i?(o)+i?W^ (36) 

where H'-°^ is given in Eq. (32) and the first-order cor- 
rection in C is 

(37) 

The subscript tp indicates an operator acting on the orig- 
inal spinor ip. 

In this method, the modified Dirac equation implies a 

modified continuity equation and hence requires a mod- 
ified scalar product. At first order in C, the continuity 
equation is 

5o[V'^(l + 7°C°)V'] + 5,[V'^°(7'' + C^'>] = 0, (38) 

and the probability density can be identified as the com- 
bination ^1'(1-|-7°C°)V'. The corresponding scalar prod- 
uct is 

(V'i,V'2)p = J d^x Vl(l + 7°C°)V'2. (39) 

Provided H^^ is time independent, the hamiltonian H-ip 
is hermitian with respect to this modified scalar prod- 
uct and so quantum-mechanical calculations can proceed. 
When C*^ is time dependent, hermiticity with respect to 
the product (39) can be restored by adding an extra term 
[40]. We thereby obtain the hermitian hamiltonian 

J?v = 5n%C'0 + ^^ (40) 

at first order in C^. 



3. Comparison 

The hamiltonians in Eq. (31) and H,^ in Eq. (40) 
typically have different forms. For our present purposes, 
however, they are physically equivalent because they give 
rise to the same eigenenergies at first order. 

To demonstrate this, first note that the difference Ail 
between the two hamiltonians can be written as 

Ail = ilW _ijW ^^(o)^_^^(o)^ (41^ 

where A is given by Eq. (30) . This can be shown by ma- 
nipulation of the field redefinition and the Dirac equation 
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or verified by direct calculation. The physical quantities 
of interest are the eigenenergies. These are obtained as 
perturbations -B^^l to the unperturbed values E^'^^ cal- 



culable at first orijfer as 



■0,X 



(42) 



where tl)^°^ are solutions to the unperturbed Schrodinger 
equation H^^^ip^^^ = E^^^'ip^^\ However, the expectation 
value of AH is zero in this scalar product, so the first- 
order perturbations to the energies are identical for both 
hamiltonians. 

The above first-order result suffices for the present 
work, although we anticipate equivalence also holds at 
higher orders in a complete analysis. For our purposes, 
the field-redefinition method proves technically and con- 
ceptually easier because the hamiltonian is hermitian 
with respect to the usual scalar product. We therefore 
adopt the field-redefinition method in the remainder of 
this work, and all references to H implicitly refer to i/^^. 



B. Relativistic hamiltonian 

The relativistic hamiltonian for the action S.^ in Eq. 

(3) can be obtained via the field-redefinition method. At 
the appropriate perturbative order, we find the operator 
A of Eq. (30) takes the form 



A = 1-1/ 



■7 



_ 3g(0,l)g(0,l)^0 



|e(o.i)rO-|v (r°7V + 7Vr° 



(43) 



Implementing the field redefinition and varying the trans- 
formed Lagrange density results in a Dirac equation from 
which the hamiltonian H can be identified. 

The hamiltonian H can be split into pieces according 
to perturbative order, 

H = if (O'O) + if(O'i) + if (I'O) + /f(i'i) + . (44) 

The component H^'^'*^^ is the conventional hamiltonian 
for the Lorentz-invariant Minkowski-spacetime limit of 
the theory. The first-order Lorentz-invariant piece is 

if(°'i) = ^iihjk + hooVjkh°j''d^ + ihjod' - 

+ \i{dihoo + d%k)l''l'- (45) 



It represents the first-order correction to the conventional 
hamiltonian H^^'^^ arising from gravitational and inertial 
effects. 



The first-order correction to the conventional hamilto- 
nian arising from Lorcntz violation can be written 

ij(i'O) = ao — meo + 2ic(jo)9^ — (mcoa — iejd^)j^ 

+i-bo - 2zrf(,o)a^)75 + {iHoj + 2gj^k0)d''h' 
-\bj+i{djkd'' + doodj) - im5™e,fei]757V 

-ie,7m(5'°°r?'" + |5'"')5fe757^'- (46) 

This result matches the one previously obtained for the 
SME in Minkowski spacetime [37] when the change in 
metric signature is incorporated. When minimal torsion 
is included in the analysis, its background value enters 
Eq. (46) through the replacement — > (6off),i specified 
by Eq. (9). It can be constrained through a reinterpreta- 
tion of experiments searching for nonzero [17, 34, 35]. 

An interesting issue is the extent to which the gravita- 
tional and inertial effects in Eq. (45) mimic the Lorentz- 
violating effects in Eq. (46). For example, in a rotat- 
ing frame of reference the term h^'^ ejki'y5j^j'' / i in Eq. 
(45) contains a coupling of the rotation to the spin of 
the particle with the same operator structure as the 
term —bjjsj^'y^ in Eq. (46). At this order, a frame 
rotation can therefore mimic potential signals arising 
from a nonzero coefficient bj for Lorentz and CPT vi- 
olation. This effect has been observed in tests with a 
spin-polarized torsion pendulum [41]. The same term in 
Eq. (45) also contains gravitomagnetic effects that are in 
princnplc! observable in tests searching for 6^ if siifficient 
sensitivity is reached. Certain Lorentz- violating effects 
can be separated from gravitational and inertial effects 
because the former generate time-varying signals due to 
the motion of the Earth and can have flavor dependence, 
but a complete separation may be problematic. 

The 0(1,1) contribution to H can be separated as 



^(1,1) = 



(1.1) 



(1,1) 



b ^ ^^c ^ ^^d 
(1>1) _L fj-(l,l) _L E7(1>1) 



(47) 



Here, the term H^'^^ arises from Lorentz-violating cor- 
rections to the metric fluctuation h^^. It is given by 

+¥ih^'^+h^o'^nok)l°l^d=- (48) 

The other terms in iJ*^^'^) are labeled according to the 
type of coefficient for Lorentz violation involved. The 
contributions involving the four-component coefficients 



and b^ of mass dimension one are 



(49) 



and 



^b'"'^ = (-^o + &^M75 

+ (6, - p^oo - ^^'/i,fe)7^57^ (50) 



10 



respectively. In the latter equation, effects from minimal 
torsion are included via the replacement 6^ — > (bcs)^ 
given in Eq. (9). 

The contributions involving the dimensionless coeffi- 
cients Cf_ii, and dij,i, are 

iJ^i'i) = i[coo??jfc +Ckj + 5C00 {hmVik - hjk) 
+ 2c(;o)/i'°'7ifc + jCko^jo — jCjohko 
- \cij {hoov'k + h\) - CMh\]-i^-i^cP 

-to(coo + ^Coo/loo + 2c(jo)/i-'°)7° (51) 

and 

+1000 + + 2d^ko)h^°)l'^-i5-i^ dj 

+i\dkj + jdkohjo — jdjohko — ^doohjk 

+m\>fj^k - -4j/i^°)757' 
-i2mrf-'"/i'=0e,fezV, (52) 

respectively. The dimensionless four-component coeffi- 
cients and generate the expressions 

-meo + me^hjo + jmeo/ijo7°7'' (53) 

and 

(54) 

while the dimensionless coefficient g\^v leads to 

-im(r'=" +5'""/^^'™ +r""/imo)eifea°757' 
+ im%fco/i'°7°7^' + >r'°/i'°e,fc;75. (55) 

Finally, the antisymmetric coefficient i/^^ of mass dimen- 
sion one contributes 

-i{HjO + Whjk +Hjkh''°)-f'. (56) 

In the above expressions, the fluctuations of the various 
coefficient fields appear in H only at perturbative order 



0(1,1). For most purposes, it is necessary to find expres- 
sions for these fluctuations prior to using the hamiltonian 
H in a given analysis. This issue is addressed further in 
Sec. IV C, where the spin- independent coefficient fluctu- 
ations (acs)ti and c^^ are considered in more detail. 

The remaining piece of the hamiltonian H lies at per- 
turbative order 0(0,2) and represents the second-order 
Lorentz-invariant contribution from gravitational and in- 
ertial effects. It takes the form 

= i{\hQohoorijk + \hiQh!-°r}jk-\hoohjk 
-lh,ih\)f^^d^-ihkoy''dj 
-m{\hoohoo + i/i,o/i^°)7°- (57) 

C. Nonrelativistic hamiltonian 

Most experimental tests of interest in this work 
are nonrelativistic. In this section, we use a Foldy- 
Wouthuysen transformation [42] to extract from the rel- 
ativistic hamiltonian H the parts of the nonrelativistic 
hamiltonian i?NR relevant for the subsequent analyses. 

The Foldy-Wouthuysen transformation is a system- 
atic procedure for determining the nonrelativistic content 
of certain relativistic quantum-mechanical hamiltonians. 
For a massive four-component Dirac fermion, the trans- 
formation generates a series expansion in powers of the 
fermion momentum. In the present case, the transfor- 
mation can be implemented as usual, but care must be 
taken to keep track of both the order in momentum and 
the perturbative order 0{m,n) in coefficients for Lorentz 
violation and in gravitational fluctuations. 

Performing the Foldy-Wouthuysen transformation for 
the complete hamiltonian H of Eq. (44) is cumbersome 
and also unnecessary for our scope because most attained 
sensitivities to spin couplings are unlikely to be improved 
by studying the suppressed effects from gravitational cou- 
plings. However, only limited sensitivity currently exists 
to spin-independent effects controlled by the coefficients 
ttfi and Bfj, because these are unobservable for baryons and 
charged leptons in Minkowski spacetime. In the remain- 
der of this work we focus on general spin-independent ef- 
fects, which are associated with the coefficients a^, Cp,^, 
and e^. Since the minimal torsion coupling also involves 
spin, this focus implies also disregarding nonrelativistic 
effects due to torsion, effectively restricting attention to 
the limiting Riemann geometry. Although beyond our 
current scope, a Foldy-Wouthuysen analysis incorporat- 
ing spin-dependent effects could lead to additional tor- 
sion sensitivities beyond those obtained via searches for 
[3, 17, 34, 43]. 

In the relativistic quantum theory, the upper two com- 
ponents of the four-component wave function describe 
the particle while the lower two describe the antiparti- 
cle. The hamiltonian H can be separated into an odd 
part O containing terms that mix the upper and lower 
components and an even part £ that involves no mix- 
ing. The idea of the Foldy-Wouthuysen method is to 
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find a momentum-dependent unitary transformation S 
in the Hilbert space such that the 4x4 hamiltonian 
H = e^'^ife"'^ js 2x2 block diagonal. The leading 
2x2 block of H then represents the desired nonrela- 
tivistic hamiltonian -ffNR- The full transformation .S* is 
obtained at the desired level of accuracy via an iterated 
series of incremental transformations reducing the off- 
diagonal content to the appropriate order. 

We proceed by separating the hamiltonian H into an 
odd part 0^™'"^ and an even part ^q™'"-* at each per- 
turbative order 0(m,n). A subscript is used to specify 
the iteration number of the transformation, with corre- 
sponding to the zeroth iteration. The relativistic hamil- 
tonian H can therefore be written as 



H. = H 



^^U ^ 0(0,0) ^ 0,.,., ^ 0,.,>v ^ 0^ 



(0,1) 



i(i,o) 



1(1,1) 



(0,0) 



(0,1) 



(1,0) 



(1,1) 



(58) 



The Foldy-Wouthuysen sequence is then defined itera- 
tively as 

= IZfci [^Sn^:_[iS^ Ho] •••]], (59) 



fe=0 



k commutations with iSn 



where 



Sn 



2m 



(60) 



At each stage, the sum on k is truncated once the ap- 
propriate order in momentum and small quantities is 
reached. The iteration continues until the hamiltonian 
is even at the desired order. Here, we proceed to 0(1,1) 
in the small quantities and to second order in the momen- 
tum, which requires three iterations and yields a hamil- 
tonian H^. 

The desired spin-independent contributions to the non- 
relativistic hamiltonian i^NR can be separated according 
to perturbative order and origin as 



_ o-(0,0) 
= 111 



NR 



Hi 



(0,1) 
NR 



(1.1) 

NR.aett 



, H-(l.O) 

r(i.i) 



^ -"NR,c ^ -"1 



(1,1) 
NR,h 



^ -"NR • 



(61) 



Here, H^^^ is the conventional Minkowski-spacetime 

hamiltonian. The conventional Lorentz-invariant contri- 



butions H, 
ten as 



(0,1) 
NR 



due to the metric fluctuation can be writ- 



er, 



(0,1) 
NR 



'Ofc 



1 
4m 



2m 



(62) 



The leading-order perturbation H^-^^ due to Lorentz 

violation and independent of /i^^ is identical to the 
Minkowski-spacetime result given as Eq. (4) of Ref. [37]. 



The corrections of primary interest for our purposes lie 
at perturbative order 0(1,1). The contribution from U/j, 
and Cfj, can be written in terms of the effective coefficient 
(cteff)(u introduced in Eq. (16), and it takes the form 

t(1,1) 
'NR.aett 



Ht^l = (Seff)o + (aeff)fe/l°'= - -(aeff)^/l,fe/ 



1 



+—{(acs)]-^{acs)jhoQ)p^. (63) 
The 0(1,1) contribution from c^i, can be written 



(1,1) 

NR,c 



Ok\ 



-TO (coo + ^Coo/loo + 2c(fe0)/l 

-2 (c(jo) +c^jk)h°'' - C(ofe)/ij*') f' 

~~il^O'njk +Cjk + jCoohooVjk + C((o)/l°'??jfe 

- ^Cjkhoo - ^coohjk - 2c(ji)h\)p' p'^ , 

(64) 

while the 0(1,1) contribution from Lorentz- violating ef- 
fects on the metric fluctuation is 



H-(l.l) - _lrr,/,(l'l) J-/,(1>1)„2 

2to, ■> 



(65) 



The remaining contribution to f/^NR is the 0(0,2) con- 
tribution involving quadratic products of /i^u^. This can 
be written as 

^nr' = -hm{hojh°^ + lhoohoo)+hojh^''pk 
- — {jehoohoorjjk + ^hoih°^rjjk 

- Ihoohjk - \hjih\)pipK (66) 



IV. CLASSICAL THEORY 

For many analyses of Lorentz violation in matter- 
gravity couplings, a classical description suffices. This 
section considers the classical limit of the quantum the- 
ory discussed above, focusing on the limit involving the 
coefficient flelds (aeff)ju and c^j/. A suitable classical rela- 
tivistic action for a point particle is presented, and its ap- 
plication to modeling test and source bodies is described. 
The modifled Einstein equation and the equation for the 
trajectory of a test body are obtained. We also discuss 
the treatment of the coefficient fluctuations (oeft)^, c^i/ 
and the procedure for determining the background grav- 
itational fleld in the presence of Lorentz violation. 



A. Particle action 

The classical action 5c corresponding to the action S 
of Eq. (1) can be written as 

S^ = Sg + Su + S'. (67) 
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As before, Sq describes the gravitational dynamics, while 
S' contains the dynamics associated with the coefficient 
fields for Lorentz violation. The partial action Su is the 
classical relativistic point-particle limit of the action 
for the fermion sector. In this subsection, we discuss Su 
and extend it to describe test and source bodies. 

1. Point particle 



sion. This follows from the adoption of {acs)fi with lower 
index as the coefficient field, which implies that the con- 
traction {aca)^u'^ in Eq. (68) contains no metric. Work- 
ing with (oeff)'' instead is possible but less convenient. 
It would produce a contribution to T^"^ along with cor- 
responding changes in the contributions to the energy- 
momentum tensor T''^'^ associated with S', leading to 
the same physical results. 



At leading order in Lorentz violation, we find 

Su = j dX (^-m^J-{g^,„ + 2c^j.)u''u'' - (aeff),,u^ 

(68) 

In this expression, the particle path = x^{X) is 
parametrized by A, and = dx^^ / d\ is the four- velocity 
of the particle. As usual, a gauge choice for A is required 
to fix the path-reparametrization invariance and to de- 
fine the proper time of the particle on shell. We adopt 
here the conventional proper-time interval 



dr = —g^^ydxi^dx^ 



(69) 



The leading-order form (68) of the classical action can 
be deduced in several ways. At the intuitive level, the 
term involving {acs)ii has the same structure as the usual 
coupling of a classical relativistic particle to an electro- 
magnetic 4-potential, and this is consistent with the cou- 
pling of afjt in the field-theory action (3). Similarly, the 
coefficient c^i, enters Eq. (68) as a shift in the metric, 
which is compatible with the way it appears in the field 
theory (3) . In a different vein, the contributions from c^v 
to the relativistic particle action have previously been 
discussed in the context of the photon sector [22], where 
the appearance of c^^ as a metric shift is related to the 
coordinate choices discussed in Sec. II C 2. The valid- 
ity of the action (68) can also be verified by extracting 
the leading-order terms from the all-orders expression ob- 
tained by construction of the exact relativistic dispersion 
relation [44]. In the present context, we can demonstrate 
explicitly that the action (68) reproduces the correspond- 
ing terms in the nonrelativistic hamiltonian iJNR gener- 
ated from the Foldy-Wouthuysen transformation as Eq. 
(61). This involves expanding the action (68) to the 
appropriate orders in velocities and Lorentz violation, 
extracting the conjugate 3-momentum, constructing the 
corresponding hamiltonian, and matching it to I^nr in 
Eq. (61). These methods all confirm that Eq. (68) is 
the correct leading-order form of the relativistic classical 
action. 

The energy-momentum tensor T^" for the point parti- 
cle can be derived from the action (68) by variation with 
respect to the metric, as usual. We obtain 



J 



^_mu^u^6\x-x'ir)) 

where the proper-time interval is given by Eq. (69). Note 
that no contributions from (aeff)^ appear in this expres- 



2. Test and source bodies 

The experiments and observations considered in this 
work involve bodies B acting as test bodies T or as 
sources S. Many of these bodies consist of atoms or 
macroscopic matter rather than individual particles. It is 
therefore useful to extend the point-particle action (68) 
to an action for a body B. This requires consideration 
of several issues. 

One issue arises because the interactions involved in 
binding electrons, protons, and neutrons into atoms 
and macroscopic matter contribute additional Lorentz- 
violating effects. This issue appears also in the study 
of fermion-sector SME coefficients in the Minkowski- 
spacetime limit [45]. However, for the gravitational tests 
of interest here, it is reasonable to assume that these in- 
teraction effects arc small compared to the propagation 
effects. 

Another issue arises from the spacetime dependence 
of the coefficient fields (ftes)^ and c^^, which implies 
Lorentz-violating effects may vary over the region filled 
by the body. Most of the test bodies we consider are 
small, so it is reasonable to approximate the coefficient 
fields as constant across the extent of the body. This 
corresponds to the usual approximation of constant met- 
ric fluctuation h^i, across a test body. However, some of 
the source bodies we consider are comparatively large, so 
some variation of the coefficient fields over the source is 
plausible. This could produce Lorentz-violating effects of 
various types, including possible dependence on the mass 
moments of the body. In what follows, we suppose that 
the variation of the coefficient fields is sufficiently mild 
and smooth that these effects can be neglected for the 
bodies we consider. A more comprehensive treatment of 
this issue would be of potential interest but lies outside 
our present scope. 

With the above assumptions and for most purposes in 
this work, a given body B can be modeled as a composite 
particle with constituents located at a single spacetime 
point and having the same 4-velocity, held together by 
binding energy. The body B can then be assigned an 
effective mass mP , expressed in terms of its constituent 
particles as 



m 



/B 



(71) 



Here, w ranges over the particle species forming the body 
B. For example, w can be taken to include the electron e. 
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the proton p, and the neutron n whenever B is an atom or 
made of ordinary macroscopic matter. The symbol 
denotes the number of particles of type w in the body, 
and m'^ represents the contribution to the mass from the 
binding energy. In practice, the exact values of are 
readily obtained for test bodies on the atomic or molec- 
ular scale, while estimating for macroscopic test and 
source bodies in the laboratory is straightforward. When 
considering the Earth as the source body, wc adopt the 
estimates = A^^ « = 1.8 x 10^^ based on recent 
studies of the bulk Earth composition [46]. The differ- 
ence A^^ — ~ 10**^ is primarily due to the iron core. 
The radial variation in neutron content is neglected in 
what follows, although it might be of interest in more 
detailed studies. 

Similarly, the Lorentz-violating properties of B can 
be represented via effective coefficient fields (a^g)^ and 
for the body. These can be viewed as the sum of 
vacuum values and coefficient fluctuations, 

(aeff)/x = (ae'ff)M + i^ce)^, (c^)m- = + (f'')^., 

(72) 

in paralkil with the point-particle case. The form of the 
action (68) implies that the coefficient field (a^^)^ takes 
the form 

w 

where {a'^)i^ is a possible coeflicient field associated with 
the binding energy that contributes to {afg)^. Also, ex- 
panding the action 5„ for small Lorentz violation shows 
that at leading order the body coefficient field (c^)/ni/ can 
be taken as 




where (c'^)pt/ is associated with the binding energy. 

The two contributions {a'^)^ and {c'^)ij,i, describe 
Lorentz violation arising from the particles associated 
with the forces binding together the body B. These par- 
ticles are primarily gravitons, gluons, or photons and are 
associated with boson fields, for which the CPT-violating 
terms axe expected to be small or zero. In the mini- 
mal SME, no such terms exist for gravitons, while for 
photons and gluons they can reasonably be assumed to 
vanish [3]. Also, the relevant photon coefficient {kApY 
is constrained well below levels relevant for this work [2] . 
Possible CPT-violating contributions from other sea par- 
ticles largely cancel due to particle-antiparticle pairings 
or are suppressed in loops involving weak interactions. 
We therefore approximate the contributions from {a'^)fj, 
as negligible, 

Kl)/. ^ 0, (75) 

in this work. In contrast, all the force fields have CPT- 
even terms that can be expected to contribute to (c'^)^i^. 



so the resulting size of (c'^)^^ may well be of the same 
order as (c"")^!, and cannot be neglected. 

Given the above discussion, wc conclude that the 
leading-order approximation to the classical action 
for a body B can be written in the simple form 

-{a^n)^u'^), (76) 

where m^, (a^fj)^, and {(^)^v are given by Eqs. (71), 
(73), and (74), respectively. In this expression, u'^ is 
the 4- velocity of the body B, which follows a world line 
parametrized by A. This form is convenient for calcula- 
tional purposes. Note, however, that the derivation es- 
tablishes validity of this form of the action only at leading 
order in (c^);^,.. 

The model action (76) for a body B suffices for most 
situations of interest in this work. In a few cases where 
the body acts as a gravitational source S, it is also useful 
to incorporate dominant effects arising from its rotation. 
For this purpose, we treat S as rigid at leading order and 
assume that the distribution of electrons, protons, and 
neutrons is approximately uniform throughout it. For the 
bodies we consider, this assumption is good to within an 
order of magnitude. The density /; of S can bo taken as 
the mass per unit volume and approximated as uniform. 
For large source bodies such as the Earth, some results 
could in principle also depend on spherical moments of 
inertia [7], but these effects are neglected here. The an- 
gular velocity u of rotation is defined in the frame at rest 
relative to S with origin at the center of mass, which can 
be identified with the location of S. 



B. Equations of motion 

The primary experimental observables arising from the 
classical theory involve the relative motion of particles. 
To investigate the motion of a test particle in the pres- 
ence of gravitational sources, the modified Einstein equa- 
tion must be solved for the background metric and the 
equation for the particle trajectory must be found. In 
this subsection, we derive the equations of motion from 
the action (67) in terms of the metric fluctuation /i^^ 
and the coefficient fluctuations (a^fj)^, ((?)ixv The issue 
of expressing these fluctuations in terms of the vacuum 
values (aj^ff)^, (c^)^jy for a given distribution of matter 
is addressed in the following subsection. Sec. IV C. We 
conclude the present subsection with comments about 
the implications of Lorentz violation for the equivalence 
principle. 
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1. Modified Einstein equation 



written as 



Varying ttie action (67) in Riemann spacetime with re- 
spect to the metric yields the modified Einstein equation 



-r ^ v^v^^ 

^ (0,1) a/3" " 



(77) 



where G^"^ is the Einstein tensor and the terms on the 
right-hand side form the energy-momentum tensor. The 
contribution T^'^ arises from Lorcntz violation in the 
pure-gravity sector. The energy-momentum tensor T^"^ 
for the matter is given in Eq. (70) . The remaining energy- 
momentum contribution T'a"^ arises from the dynamics 
of the coefficient fields for Lorentz violation and is deter- 
mined by S'. 

Taking the covariant divergence of Eq. (77) and us- 
ing the Bianchi identities shows that the geometry re- 
quires the total energy-momentum tensor to be locally 
conserved. The theory can be consistent only if this re- 
sult is compatible with the explicit form of the energy- 
momentum tensor. This requires careful accounting of 
contributions from the massless NG modes arising from 
the spontaneous Lorentz breaking [3] . In the general case, 
these modes are contained in the fluctuations t\^iv...- 

For the pure-gravity sector, the relevant analysis is 
given in Refs. [3, 7] and can be subsumed as needed 
in the present context. For the matter sector, the NG 
modes produce no relevant contribution to the energy- 
momentum tensor Tj^^ at the post-newtonian order ap- 
propriate for the tests considered here. The; key point 
is that the leading Lorentz-violating effects of coefficient 
fluctuations tx^^y,,, arise at PN0(2) or beyond, as shown 
in Sec. 11 D. Since these fluctuations are accompanied 
by an additional factor of Gjv in the modified Einstein 
equation (77), they affect the metric only at PN0(4) or 
beyond. However, for the tests considered below it suf- 
fices to work at PN0(3) for Lorentz-violating terms, so 
the coefficient fluctuations (Sg^)^ and (J^)^i/ appearing 
in T^'' can be neglected in Eq. (77). 

In contrast, the contributions to T'^^ arising from the 
dynamics of the coefficient fluctuations t\p,v... are of po- 
tential relevance in solving the modified Einstein equa- 
tion for the metric. The specific effects associated with 
the coefficients and ^^jxv are derived in Sec. 

IV C. 



+2(cT)(„^)r(o + d^i:^\^u'-u^ 

_ J_[aA'(ST^)„ - v'^^d^{^l,)0]u-, (78) 

where each dot on represents a derivative using the 
proper-time interval (69). The superscript T denotes 
quantities associated with the test particle. The first 
term on the right-hand side is the usual geodesic contri- 
bution, where F^g i)^a0 linearized Christoffel sym- 
bol. A Christoffel symbol with subscript (1,1) also ap- 
pears in Eq. (78). It is defined as the linearized Christof- 
fel symbol with hf^, replaced by /i^V^- This introduces 
matter-sector coefficients associated with the gravita- 
tional source, along with any gravity-sector coefficients 
that may be included in the analysis. 

Once the forms of (ftcff)/^, and c^,^ arc established, 
Eq. (78) can be used to determine to 0(1,1) the motion of 
a classical test particle in a curved but asymptotically flat 
spacetime with nonzero coefficients for Lorcntz violation 
Ufi, e^, and c^i^. Obtaining expressions for /i^j^, (aeff)/i) 
and c^i/ is the subject of Sec. IV C. 

Although unnecessary for the present work, we can 
comment in passing about the effects of nongravitational 
interactions on the particle trajectory. Any such interac- 
tions can be viewed as introducing an additional contri- 
bution to the right-hand side of Eq. (78). Using the 
perturbation scheme of Sec. II D, this additional acceler- 
ation can be expanded as a sum over terms a'^^ , 
one at each perturbative order 0(m, n). Notice that, al- 
though the interaction itself is nongravitational, contri- 
butions to ^ with n ^ can be induced from grav- 
itational couplings in the interaction sector. Similarly, 
Lorentz-violating contributions to a^*^ can originate 
from coefficients for Lorentz violation in the interaction 
sector. If we also expand as 



_ 4^ -I- -V^'^ -I- 

~ (0,0) ^ '^(0,1) (1,0) ^ (1,1) 



(79) 



then we obtain the following additional terms for the par- 
ticle 4- acceleration x^: 



2. Particle trajectory 

The equation of motion for a classical test particle T is 
obtained by varying the action (67) with respect to the 

particle position 4- vector x^^ . In the absence of Lorentz 
violation, this is the geodesic equation. However, in the 
presence of Lorentz violation, the trajectories of test par- 
ticles T no longer match the geodesies of the spacetime. 
Expanding to 0(1,1), the equation of motion can be 



^(0,0) "(0,0)' 
^(0,1) °^(0,1)' 

^(1,0) ^ "(1,0) - 2< (c )(a/3)a(o,o) 
-2(cT)(„^)a^„,)u'5«^ 



^(1,1) ^ "(1,1) )(a«"(0,l) 

^2ic^)^^f,)af,^,yu'^ + 2/i^"(cT)(„^)aJ^ gj 
-2r?''«(JT)(„^)af,_„^. (80) 
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The trajectory at 0(1, n) is affected both directly by 
a^j^ and indirectly by combinations of a'^^ with the 
coefficients for Lorentz violation. The origin of the indi- 
rect terms can be traced to the additional factor of 2c^,y 
in the action (68) relative to the conventional proper-time 
interval (69). 

3. Implications for the equivalence principle 

The deviations from geodesic motion implied by Eq. 
(78) can be species dependent because the couplings to 
the coefficient fields (acn )/j and c^,/ can vary with parti- 
cle flavor. This leads to apparent violations of the weak 
equivalence principle (WEP), which stipulates that the 
motion of uncharged test particles is independent of in- 
ternal structure or composition [8]. 

One implication of this observation is that experiments 
designed to test the WEP are also sensitive to {aes)^ and 
c^i,. Since all the WEP violations implied by Eq. (78) 
arc accompanied by effects associated with the breaking 
of rotation and boost symmetries, the experimental sig- 
natures associated with {aes)^ and c^^ typically differ 
from those in other scenarios for violations of the WEP. 
The latter portion of this work discusses in some detail 
the role that experiments testing the WEP can play in 
searches for Lorentz violation. 

The flavor dependence of the coefficient couplings leads 
to the philosophical question of whether spontaneous 
Lorentz violation in the matter sector violates the WEP 
or merely mimics violations of the WEP. The issue hinges 
on the interpretation of the term 'uncharged test parti- 
cle.' In models with spontaneous Lorentz violation in 
the matter sector, the NG modes couple to test parti- 
cles and so mediate an interaction. This interaction can 
be identified with Einstein-Maxwell electrodynamics [14], 
GR gravity [27], an effect on torsion [14], or a new force 
[11, 29, 30]. If the term 'uncharged' is taken in the re- 
strictive sense to mean that the test particle is unaffected 
by standard forces such as electrodynamics, then the tra- 
jectory deviations of 'uncharged' test particles caused by 
nonzero (flcfr)^ and c^^ coefficients represent violations 
of the WEP. If instead the term 'uncharged' indicates 
the test particle cannot have nongravitational couplings 
of any kind, then no violations of the WEP occur. How- 
ever, in this latter case 'uncharged' test particles may 
be nonexistent in the matter sector of the SME, where 
generically all particles experience nonzero (Seff)^ and 
c^j, coeflacients. We emphasize that the above discussion 
is a matter of philosophical classification only, without 
impact on the practical issue of using tests of the WEP 
to search for Lorentz violation via the deviations from 
geodesic motion described by Eq. (78). 

The WEP is subsumed in certain other equivalence 
principles, such as the Einstein equivalence principle or 
the strong equivalence principle. These incorporate also 
aspects of local Lorentz invariancc and local position 
invariance. Since nonzero coefficients (aeff)/*) c^i/ cor- 



respond directly to local Lorentz violation, and since 
Lorentz violation can be position dependent, the devia- 
tions from geodesic motion described by Eq. (78) can rep- 
resent violations of these broader equivalence principles 
arising in more than one way. Several related philosoph- 
ical issues remain open, including classifying violations 
of various equivalence principles according to properties 
of the coefficients for Lorentz violation and identifying 
implications for relations such as the Schiff conjecture 
[8, 47]. We note also in passing that comments analo- 
gous to those above bear on the philosophical issue of 
whether theories with matter-sector couplings to spon- 
taneous Lorentz violation constitute metric theories of 
gravity. 

C. Coefficient and metric fluctuations 

To solve the equation of motion (78) for the trajectory 
of a test particle, explicit expressions for the metric ffuc- 
tuation /i^^ and the coefficient ffuctuations (Sefj)^, '^^v 
are required. Within a specific model with known ac- 
tion S' for the coefficient fields, these expressions can be 
obtained by direct calculation. An illustration of this is 
provided in Sec. V. However, in the interest of general- 
ity, it is useful to establish results valid for a large class 
of models. In this subsection, we outline a procedure to 
obtain expressions for h^v and for the generic coefficient 
fluctuations txfn^... when S' is largely unknown, and we 
obtain explicit results for /i^j^, (ueff) fi, and c^j^ applicable 
to the equation of motion (78). These results are used 
in later sections of this work in establishing experimental 
signatures for Lorentz violation. 

1. Methodology 

Consider first the metric fiuctuation h^i,. In the per- 
turbation scheme of Sec. II D, the expansion of hf^v takes 
the form (26). To determine the tcst-particle trajectory 
at order 0(1,1) via Eq. (78), it is necessary to obtain 

exphcit expressions for /i^V^ and /i^V'- 

The Lorentz-invariant contribution h^/^^^"^ can be ob- 
tained in the usual way as the leading-order solution 
of the Einstein equation, taking the Lorentz-invariant 
part of the energy-momentum tensor as the source. To 
PN0(3), the standard solution can be written in har- 
monic coordinates as 

= 2c/, 4■■'^ = -4y^ i^t'^^^up^m 

where U and are the usual post-newtonian potentials 
defined as 
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In these expressions, the density p{x\t) and the 3- 
velocity (x', t) are properties of the source in the chosen 
asymptotic inertial frame. As described in Sec. IV A 2, it 
suffices in this work to use the approximation (76) for the 
source-body action , so the energy- momentum tensor 
takes the generic form (70). 

The Lorentz-violating component can conve- 

niently be viewed as a sum over individual contributions 
arising from each coefficient field for Lorentz violation, 

h^^u'^ = {h^:''\. + {h'^''\. + ... + {h^^''\.. (83) 

This means each coefficient can be treated in turn. How- 
ever, the procedure for determining a particular contri- 
bution can take different paths depending on the type of 
gravitational coupling of the coefficient field. 

A simple case arises for any coefficient field tx^i,,,, that 
is minimally coupled to gravity. Then, only the vac- 
uum values tx^u... in the expansion (11) contribute to the 
energy-momentum tensor T^j^" of the source at the rele- 
vant order. The key point is that the solution for /iJiV^ 
at PN0(3) arises from the combination kT^^'^, which it- 
self already lies at PN0(2). However, as discussed in 
Sec. II D, the coefficient fiuctuations of interest here are 
'^X/M^... PN0(2). These fluctuations there- 

fore cannot contribute to below PN0(4). As a re- 

sult, ft^y^^ can be found directly by solving the modified 
Einstein equation with attention limited to the vacuum 
values txf,,y.... 

We note in passing that this procedure is consistent 
with the no-go result for explicit Lorentz violation in 
gravity [3] even though tx^iv... can be interpreted as a 
coefficient for explicit breaking. This is because we are 
working to 0(1,1) and PN0(3), for which Datxixv... ~ 
0(/it) ^ 0(1,1). As a result, the covariant deriva- 
tive of kT^H;" is compatible with the Bianchi identities 
D^G^'^ = at this pcrturbative order. In effect, the 
comparatively low perturbative order implies that a con- 
stant vacuum value remains consistent with the geometry 
of spontaneous Lorentz breaking. It is also noteworthy 
that an independent contribution from T'^'^ may exist 
that satisfies local conservation and hence is compatible 
with the Bianchi identities to the relevant perturbative 
order. This would also represent a consistent theory, al- 
beit a different one. The two theories involving the co- 
efficient field (ccff)^ with minimal and with nonminimal 
gravitational couplings provide an illustration of this, as 
is discussed in the next subsection. 

In the simple case with minimal coupling, once 
has been found, it remains only to determine the direct 
contributions to the equation of motion (78) arising from 

the fiuctuations i^xiiv...- ^'^'^ ^^^^ purpose, we can apply 
the requirement that the system of the source S and the 
test body T conserves the total 4- momentum P^. For a 
two-body system, this implies the force law must be an- 
tisymmetric upon exchange of S and T. Otherwise, the 
forces on each body due to the other would violate New- 
ton's third law, and the system would self-accelerate. 



At PN0(2) and in the absence of Lorentz violation, 
the relevant force between S and T can be directly iden- 
tified as mx^ . At higher order and in the presence of 
Lorentz violation, it is simpler to impose conservation of 
the total 4-momentum, dP^/dt = 0. In principle, can 
be found by adding the conjugate momenta of S and T 
obtained from the two-body action. In practice, for the 
perturbative order to which we work, it suffices to obtain 
the conjugate momentum for T alone and require anti- 
symmetry of its time derivative under the exchange of S 
and T. The constraints fixing the fiuctuations ^xiiv arise 
from the /i = j components of dP^/dt — 0, all at PN0(3) 
except for terms at PN0(4) involving the velocities of 
both S and T. We remark in passing that extending this 
treatment to higher perturbative orders requires also in- 
corporating back-reaction effects on the metric, including 
gravitational radiation. 

The above procedure holds for coefficient fields that 
are minimally coupled to gravity. If nonminimal curva- 
ture couplings also occur in 5', then additional terms 
involving the coefficient fields txfj.u... can appear in the 
energy-momentum tensor and hence can affect the mod- 
ified Einstein equation. The curvature couplings inter- 
twine the kinetic contributions from h^^ and txp.i,..., so 
txfiu... can contribute to the solution for h^jju^^ at PN0(3). 
To proceed without specifying S' , we therefore need ad- 
ditional information about txuv...- 

In the present work, the necessary information can 
be extracted from the general structure of the equation 
of motion for txixv... and the symmetries of the theory. 
When linearized, this equation of motion can be writ- 
ten as the sum of a differential operator acting on txuv... 
and a source term at most linear in /i^,^. The differential 
operator can involve arbitrary powers oitx^lv... but is in- 
dependent of h^v The contributions ^xiw interest are 
at 0(1,1) and hence are linear in both tx^^v... and h^^. At 
this order, the solution for 'i^^^^J can therefore be writ- 
ten as a sum of terms, each containing up to one power 
of h^u along with some number n of powers of tx^„... in 
the numerator and n — 1 powers of txixv... in the denom- 
inator. This expansion oiT^'^^J in terms of txiiv... and 
h^v is constrained by two requirements. One arises from 
the restriction of tx^iv... to NG modes, which must main- 
tain the extremum of the action. The solution for 'i^xi^! 
must therefore obey the NG conditions at 0(1,1). The 
second is the requirement that tx^iv... must transform as 
expected under diffeomorphisms, as a consequence of the 
spontaneous nature of the symmetry breaking and the re- 
quirement of observer general coordinate invariance [48] . 
It turns out that these two restrictions suffice to express 

^Xixv... terms of tx^u... and h^v in the cases of interest 
here. 

Once the expression for t^x/l^ ^as been found, hfj,^, can 
be obtained by combining information from the modi- 
fied Einstein equation and the trajectory equation. The 
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modified Einstein equation yields directly the piece of 
h^i, arising from tx^^,,, in the energy- momentum tensor. 

Inserting this result and the expression for ^^^^J in the 
trajectory equation and imposing conservation of the to- 
tal 4-momentum of the source and test body as before 
determines the missing piece of /z^^ arising from t^^^ . 

In both minimal and nonminimal cases, the net re- 
sult of the above procedure is a form of the trajectory 
equation in which h^v and fA/xj/... can be replaced with 
specified gravitational potentials and the vacuum values 
t\ij,v...- The solution of the trajectory equation can then 
proceed. In the next two subsections, we apply these 
methods to obtain the relevant contributions to the fluc- 
tuations from the coefficient fields (aes)^ and 

2. Fluctuations and c^v 

The treatment of the coefficient field c^^ provides an 
example involving the comparatively simple case of min- 
imal coupling to gravity. Although nonminimal curva- 
ture couplings to c^j/ could be considered, these are of 
lesser interest in the context of searches for Lorentz vi- 
olation because direct signals from c^j^ already appear 
for minimal coupling. We therefore neglect nonminimal 
couplings to c^v here. In this subsection, the relevant 

contributions to the metric fluctuation {h^c'^^)iiv are ob- 
tained to third post-newtonian order, and effects from 
Cfjiv at 0(1,1) are considered. 

Following the procedure outlined in the previous sub- 
section, we begin with the modiflcd Einstein equation 
(77). The relevant energy- momentum tensor for the 
source S is given by the expression (70) expanded to lead- 
ing order in Lorentz violation and with c^v replaced by 
Solving for the metric fluctuation to PN0(3), we 
obtain in harmonic gauge 

(/i(i'i))oo = 2(cS)oot/ + 4(cS)(,.o)y^ 

(/i(i'i))o, = -4(cS)ooT/^ 

{h(y^),k = 2ic'')o,UP\ (84) 

A consistent expansion to PN0(3) requires only PNO(l) 
terms in {hi^'^^)oj and none in {h'^'^^)jk, but we dis- 
play PN0(3) terms in (/ii^'^')oj and PN0(2) terms in 
{hi^'^^)jk because they are useful in part of the analysis 
to follow. 

The next step is to examine the contributions to the 

equation of motion (78) from c^,^ at 0(1,1). The conju- 
gate momentum can be extracted from the action (68) 
with (tteff)^ set to zero. Conservation of the total 4- 
momcntum Pf^ of the system is ensured by the require- 
ment that its time derivative be antisymmetric under the 
exchange of the source S and the test body T. We find 
that is conserved to PN0(3) without contributions 
fromJ^y. This establishes the PN0(3) result 

dx^l^V^ = 0, (85) 



showing that the NG modes associated with c^^ play no 
role at this perturbative order. 

The results (84) and (85) complete the determination 
of the trajectory equation for the coefficient field c^^. 
For a given source S, the potentials U and can be 
calculated explicitly. The effects of (c^)^iy and {c^)^v 
on the trajectory of the test body T can therefore be 
investigated in various regimes of experimental interest. 
This line of reasoning is pursued beginning in Sec. VI. 

3. Fluctuations and {aeB)fi 

For the coefficient field (oes)^, the case of minimal cou- 
pling to gravity is of lesser interest. The modified Ein- 
stein equation is unaffected because (oeff )^ is absent from 
the energy-momentum tensor (70). Also, only (&eff) ij, en- 
ters the trajectory equation. Since it is indistinguishable 
from an electromagnetic field, it provides no relevant con- 
tributions at 0(1,1). We therefore expand the treatment 
to the case of nonminimal curvature couplings, for which 
(acff)^* becomes measurable [11]. In effect, the fluctua- 
tions {actf )fj. become observable by virtue of their non- 
minimal gravitational couplings. 

Following the procedure of Sec. IV CI, the first step 
is to obtain an expression for the fluctuations (aeff)Ai^'^'' 
originating from the source S using the NG condition 
and the requirement of diffeomorphism covariance. At 
0(1,1), the NG condition can be written as 

(Seff)l'''Ha?ff)'' = i(a?ff)^/i''^(a,^ff).. (86) 

We find that the contributions to {aes)'jl'^^ consistent 
with this equation and with diffeomorphism covariance 
take the form 

M^''^ = lah^Mf^r - ia(a,^ff)^/i^ + 5^* (87) 

in harmonic coordinates. Here, the constant a is calcu- 
lable but varies with the specifics of the theory, typically 
being determined in terms of the coupling constants that 
control the nonminimal couplings. The function ^ con- 
tains effects proportional to /i^^ and (a^jj)^ that are un- 
physical by virtue of the discussion in Sec. II C 1, so it is 
disregarded in what follows. 

At this stage, the result (87) can be combined with the 
modified Einstein equation and the trajectory equation 
to determine the contributions to ft-ji^'' proportional to 
(Ogfj)^. Working at PN0(3), we find these contributions 
can be written in harmonic gauge as 

(''a-%k = ^[~a{4s)oU6^''+a{al^)oW>^]. (88) 
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Paralleling the case of the coefficient field c^j^, we have 
kept here PN0(3) terms in (/ii^'^^)oj and PN0(2) terms 
in (h'"a'^^)jk as a convenience for the analysis to follow. 
In Eq. (88), U and are the post-newtonian potentials 
defined in Eq. (82). Additional potentials U^'' and 
also appear, defined by 



jjjk 



G 



N 



G 



N 



f ^3^, p{x',t){x-x'y{x-x')'' ^ 

f ^3^/ p{^,t)vk{S',t){x-xy{x-x')'' 
J \x — x'\^ 



(89) 



The results (87) and (88) fix the form of the contribu- 
tions involving the coefficient field (aeft)^ to the equation 
of motion (78). Modifications of the trajectory of a test 
body T arising from nonzero values of the coefficients 
(a^fj)p and (ajfj)^ can therefore be studied at third post- 
newtonian order. The resulting experimental signals are 
discussed starting in Sec. VI. 



V. EXAMPLE: BUMBLEBEE MODEL 

In this section, wc examine a specific model and 
demonstrate how it fits into the general theory developed 
above. This discussion is included solely for illustrative 
purposes and is inessential to the development of the pa- 
per. In particular, the analyses of experimental signals 
in subsequent sections are independent of this specific 
model, so the reader can proceed directly to Sec. VI if 
desired. 



A. Bumblebee model 

Bumblebee models are theories in which spontaneous 
Lorentz violation is induced by a potential V{B^) for a 
vector field [49] . As an illustration of the general the- 
oretical treatment of Sec. IV, we consider here a specific 
and comparatively simple bumblebee model and study its 
matter-gravity couplings. A discussion of generic mod- 
els of vacuum- valued vectors coupled to gravity including 
references to the substantial early literature can be found 
in Sec. Ill A of Ref. [15], while some more recent papers 
are listed in Ref. [50]. Discussions of various stability 
issues with these models are given in Ref. [51]. 

The action Sb for the specific bumblebee model of in- 
terest here can be written as 

Sb = So + Sbu + Sg 

= J d'^x eCc + J dr Cbu + j d'^x eC^. 

(90) 

The form of this action corresponds to that of the general 
action (67). The term Sq is the usual Einstein- Hilbert 



action (2), with cosmological constant chosen as A = 
for this illustrative case. The term Sbu represents the 
mattor-bumblcbeo coupling, while Sg contains the bum- 
blebee dynamics, including the potential V triggering 
spontaneous Lorentz violation. 

For the classical lagrangian £ bu describing the matter- 
bumblebee coupling, we choose the expression 

Cbu = -m^-{g^, + 2C,2B^B^)ui'u" + (iB^u^. 

(91) 

Here, and ^2 are coupling constants that can vary with 
the particle species. Where needed in what follows, we 
distinguish the coupling constants for a source body S 
and a test particle T by superscripts: (f, (J, i^f, (J- 
Note that the lagrangian (91) could be viewed as the 
point-particle limit of a quantum field theory, in parallel 
with the derivation for the general theory (68). 

For the Lagrange density eCg determining the dynam- 
ics of the bumblebee field, we take 

eC'e = -^eB'^'^B^^-eV + aieB''B''R^,, (92) 

where the field strength is i?^,^ = d^B^ — di,B^. The 
coupling constant ai is sometimes written ui = ^/2k in 
the literature [3, 15]. The potential V has the form 



V = ViBi'B^, ± b^), 



(93) 



where is a real number. Where a definite form is 
needed in the calculations to follow, we adopt for sim- 
plicity the smooth quadratic potential 



V = XiBf'B^ ± by/2. 



(94) 



In any event, the potential is assumed to induce a nonzero 
vacinim expectation value for the bumblebee field, which 
we denote by b^ = {B^) following standard usage, where 
b'^bn = =F&^- Denoting the bumblebee fluctuation about 
the vacuum value by B^, we can expand 



B^ = b^ 



B, 



(95) 



in parallel with Eq. (11). 

A match can be made between the bumblebee action 

(90) and the general action (67) by identifying the various 
coefficient fields for Lorentz violation with specific com- 
binations of the bumblebee field. The term C bu corre- 
sponds to nonzero coefficient fields (aes);* and c^^, given 
by 



C2{B^B,y - jgfii.BaB"). 



(96) 



It is also necessary to introduce an additional scalar field 
fc, defined as 



(97) 
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which normahy can be disregarded in the SME context 
because it is Lorentz invariant. In the presence of k, the 
general action (68) is sUghtly modified, with lagrangian 
now given by the expression 

(98) 

The term Cg yields nonzero coefficient fields s^"^ and u 
in the pure-gravity sector, given by [7] 



(99) 



Note that in this model only a single field underlies 
all the coefficient fields {aes)^, c^i/, k, s^i,, and u. It fol- 
lows that searches for Lorentz violation that are sensitive 
to any one of these coefficient fields could provide infor- 
mation constraining the others, at least in part. This 
special feature of the bumblebee model may not extend 
to models with more complicated field structure. Note 
also that the coefficient field k in the matter sector is a 
species-dependent analogue of the coefficient field u in 
the gravity sector. A nonzero value of k can introduce 
apparent WEP violations. Since these are Lorentz invari- 
ant, the resulting phcinomenology lacks the various time 
dependences that characterize WEP violations resulting 
from Lorentz breaking. 



B. Solving the Model 

Given the action (90) , we can illustrate by direct calcu- 
lation the correspondence between results from the bum- 
blebee model and ones from the general SME-based ap- 
proach developed in Sec. IV. For this purpose, it suffices 
to work at lowest nontrivial order in the couplings (^i , (^2 , 
and cTi. We focus here on observable effects arising from 
the identifications of (aofi)^, c^iy, and k in Eqs. (96) and 
(97). Observable effects involving s^'^ and u as defined 
in Eq. (99) are studied in Ref. [7]. 

The basic goal is to predict effects such as trajectory 
deviations for given values of (ocff);^ and c^,^. As dis- 
cussed in Sec. IV C 3, observability of {aes)^i involves 
nonminimal couplings, so in the present context we can 
expect dominant effects from (flcff)/^ to be proportional 
to the product Ci'^i- Iii contrast, dominant observable 
effects from c^^, are generated directly from (2 ■ 

We remark in passing that the special bumblebee 
model considered here is experimentally viable provided 
the sizes of Ci and C2 are compatible with existing con- 
straints on long-range spin- independent forces [33]. The 
proportionality of (acff)^ to (lai implies the model can 
yield Lorentz- violating effects involving large &^ that are 
detectable only in gravitational experiments [11]. 

In this subsection, working at the appropriate pertur- 
bative order and taking the newtonian limit where useful, 
we obtain and solve the bumblebee equation of motion 



and the modified Einstein equation. These results suffice 
to determine the trajectory equation for a test particle in 
terms of the vacuum value of the bumblebee model. 
Comparison to the general SME-based approach devel- 
oped in Sec. IV yields an explicit match for (aeff),^, c^i,, 
k in terms of the couplings Cij C2, o"! and the vacuum 
value bu- 



1. Bumblebee equation 

Consider first the equation of motion for the bumble- 
bee field, which follows from varying the action Sb- At 
the perturbative order of interest, this equation takes the 
form 



2V'K - 2axbi'R^, + Cf + . . . , (100) 



where 



dr u^S'^ix — x') 



(101) 



is the source 4-current at the relevant order. In Eq. (100), 

the prime on V denotes a derivative with respect to the 
argument, while the ellipsis indicates that soiirce terms 
proportional to exist but provide no obs(;rvable con- 
tributions to the order at which we work. Adopting the 
smooth quadratic potential (94) and the expansion (95), 
the bumblebee equation can be written 

-2\Kh''h^Kp - 2aib"Rc,^ + Cfju. (102) 

The idea is to solve this expression for the fluctuations 

B so they can be eliminated from the analysis as needed. 

The solution can be obtained in momentum space with 
the propagator chosen as a suitable Green function [7]. 
The appropriate bumblebee propagator is 



P P baP° 



4Ap2(6„p")2 

(103) 



in momentum space. Note that the additional poles in 
this and following expressions can be understood as a 
consequence of residual gauge freedom to the order at 
which we work [52]. Using this propagator, we flnd the 
solution 

LD Pa 

(104) 

where (B^J^(p), (B^J^(p), and (^^^^(p) are contribu- 

tions to 13^ (p) proportional to ci , Ci , and (^f respectively. 
In the limit of vanishing and Cf , the solution (104) re- 
duces to the known result [7] once the conversion from 

Bfj_ to Bf^ described in Sec. II B is implemented. 
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Explicitly, the quantity {3^-^)^ is given in momentum 
space by [7] 



{Ba,)n{p) = Z2— + + 



+ 



p2 4Xb"Pa P^b^Pa 

2aib"Ra^, 2aip^b"bl^Ra0 



P^ p'^b'^Pa 

For the piece proportional to Ci, we obtain 



. (105) 



p^ p^Pob^ p^Poib"- 
Cfbc^b'^Pf^p'^j, C!p,p"j,. 
p^ippbf^Y 4A(p„6")2' 



(106) 



The remaining term in Eq. (104), which contains contri- 
butions proportional to Cf, is irrelevant to the order at 
which we work. 

We emphasize that the explicit solution (104) for the 
bumblebee fluctuation is obtained by direct calcu- 
lation from the action Sb- This calculation depends 
on knowledge of the bumblebee dynamics as described 
by the Lagrange density (92). In contrast, the general 
SME-based method presented in Sec. IV to obtain an 
arbitrary coefficient fluctuation tx^„,,, replaces the need 
for complete knowledge of Sb with the judicious use of 
perturbation theory, the NG constraint, diffeomorphism 
invariance, and Newton's third law. 



2. Modified Einstein equation 

Varying the action with respect to the metric yields 
the modified Einstein equation. At leading order in B^^ 
and lowest order in h^,y, this equation takes the form 

G^i/ = 2aiK[b°'dad(^B^) + b(^^dad^)B°' - b^^uB^,) 

-Tl^ub'^do^d^B^] + 2Ky'6^6, + k{Tbu)^u, 

(107) 

where contributions from B^ are understood to be lim- 
ited to the appropriate perturbative order. Note that V 
contributes at most through massive modes at this order, 
so it plays no role in the present context. 

The matter-sector contribution {Tbu)hv to the energy- 
momentum tensor can be written 



dr 



mu^u^5'^{x — x'{t)) 
a/1 - 2Clbabfju°'uP ' 



(108) 



This explicit expression is the bumblebee analogue of the 
general form (70). As expected, terms proportional to 
are absent from Eq. (108), confirming that minimal 
couplings cannot generate Lorentz violation of the (aeff)/^ 
type in the modified Einstein equation. In this model, the 
only nonzero observables proportional to Ci arise through 
the bumblebee fluctuations (106). 



Inserting the solution for B^ at the appropriate order, 
the modified Einstein equation (107) can be solved. To 
match the analysis in the general SME-based method of 
Sec. IV the harmonic gauge must be used. At zeroth or- 
der in Lorentz violation, the conventional metric is repro- 
duced. For simplicity in this illustrative model, we limit 
consideration at the next order to the newtonian limit for 
Bp. This avoids possible complications from the residual 
gauge invariance, while permitting a complete match to 
the results of the SME method. 

The solution for the metric can be constructed directly 
from the trace-reversed form of the modified Einstein 
equation (107). Using the NG condition 



6pF = i6''6'^ V 



(109) 



and the bumblebee equation (102), we can write 

R^, = 2aiK[b"dad^^B,) - Cf 6(pi,)] + K{SBu)^,., (HO) 

where {Sbu)imu is the trace-reversed version of the energy- 
momentum tensor (108). In the newtonian limit, the 
first term on the right-hand side is higher-order in time 
derivatives and so is negligible. 

Expanding {Sbu)^^ to the appropriate order, we find 
that the 0(1,1) modifications of the metric fiuctuation 
/loo are given in terms of the bumblebee vacuum value 6^ 
by 



2G]s[bo 



(111) 



The first term arises from the bumblebee fluctuations via 
the nonminimal couplings and is the bumblebee analogue 
of Eq. (88). The second term arises directly from the 



energy-momentum tensor {TBu)tiv of the source S and 
corresponds to Eq. (84). To complete the match to the 
general SME analysis of Sec. IV, it remains to apply this 
result to determine the deviations from geodesic motion 
of a test particle T. 



3. Particle trajectory 

The equation of motion for a test point particle T in the 
presence of the bumblebee field B^ and the metric hfiv 
can be obtained by varying the action Sbu with respect 
to x^. At leading order in the fiuctuations, this yields 



-F ii°'v^ -\-r^ -\-r^ 

^ (0,1) al3"' " ^ -^Ci ^ 



C2 



(112) 



where F^p 



is the linearized Christoffel symbol. The 
terms x'^^ and x'^^ represent contributions to propor- 
tional to C,]- and Cj, respectively. The above equation is 
the bumblebee analogue of the equation of motion (78) 
obtained for the general SME analysis in Sec. IV. 
The explicit form of the quantity x^^ takes the form 



(rCi)%"""''-^5'"'(5.5„-9„S,)zi", (113) 
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where (r,;^)''^,^ contains terms proportional to (^f that 
enter via Lorentz-violating corrections to the metric. The 
piece of contributing to this equation at the relevant 
order can be written as 

{B,,)^ = aih^^b" - \aih^K" + ... (114) 

in harmonic coordinates, where the ellipsis represents 
terms that play no role within our approximations. The 
latter equation is the bumblebee analogue of the SME 
result (87) for the coefficient fluctuation (hchO^i: and as 
expected it yields contributions to the trajectory equa- 
tion proportional to the product cr ■ 
The quantity x^^ can be written 

-2(^n'"'M,Bp + b0daS,)u''u^ 

-2Cjbad^B0u'^u'^uV, (115) 

where (r(:;2)'^af3 contains terms proportional to that 
enter via Lorentz-violating corrections to the metric. In 
this equation, only the first term in Eq. (104) produces a 
relevant contribution to the fluctuation in the present 
context, which matches the SME result (85) for the co- 
efficient fluctuation Cjiij.. 

At this stage, we can verify the conservation of total 
4-momentum of the system of the source S and test body 
T, as described in Sec. IV C 1. Substituting for and 

h^Q^"^ in the trajectory equation reveals the antisymmetry 
under interchange of S and T required to satisfy Newton's 
third law. We can also complete the correspondence be- 
tween the bumblebee model and the general SME-based 
analysis of Sec. IV by making the identifications 

a = 2(71, (ocff)^ = Ci^A»> 
c^. = C2b^b, + lC2Vt..b'', fc=iC2&^ (116) 

which can be obtained by matching Eqs. (Ill) and (114) 
to the SME results (84), (88), and (87). 

VI. EXPERIMENTAL BASICS 

In the remainder of this paper, we apply the theoret- 
ical framework developed above to explore some experi- 
mental prospects for detecting Lorentz violation through 
matter-gravity couplings. As before, we adopt coordi- 
nates satisfying the condition (23), which produces sim- 
plified expressions without the photon-sector coefficients 
{kp)'^ fj^av The primary focus is on signals involving the 
coefficients (a^)/* and (c™)^,y. Certain effects associated 
with the coefficient 's^v in the pure-gravity sector are also 
considered. 

In the present section, we provide some basic infor- 
mation broadly applicable to searches for Lorentz viola- 
tion, including an outline of frame conventions and a dis- 
cussion of sensitivities to coefficient combinations. Each 



subsequent section addresses a particular class of exper- 
imental searches. Section VII examines tests with ordi- 
nary neutral matter in Earth-based laboratories, while 
Sec. VIII studies satellite-based searches with ordinary 
matter. Section IX considers more exotic laboratory and 
satellite-based t(;sts. including ones using charged parti- 
cles, antimatter, and particles beyond the first genera- 
tion. Section X addresses solar-system observations, in- 
cluding lunar and satellite ranging and measurements of 
perihelion precession. Finally, Sec. XI considers signals 
from photon-gravity couplings. 

A. Frames 

A substantial advantage of the SME framework is the 
ability to compare signals for Lorentz violation across 
a wide variety of experiments and observations. To fa- 
cilitate these comparisons, it is useful to report search 
results in a canonical inertial frame. 

In Minkowski spacetime, the canonical frame is a Sun- 
centered celestial-equatorial frame [21], which is approx- 
imately inertial over the time scales of most searches. In 
this frame, the Z axis is aligned with the rotation axis 
of the Earth, while the X ax;is points from the Earth 
to the Sun at the vernal equinox. The origin of the time 
coordinate T is the time when the Earth crosses the Sun- 
centered X axis at the vernal equinox. 

For post-newtonian investigations involving gravita- 
tional effects in the solar system, the canonical frame is 
identified with an asymptotically Minkowski frame that 
is comoving with the rest frame of the solar system and 
that coincides with the canonical Sun-centered frame [7]. 
In this Sun-centered frame, cartesian coordinates are de- 
noted by 

= {T,X-')=^{T,X,Y,Z) (117) 
and are labeled with capital Greek indices. Also, we write 
eE = (eT,ej) (118) 

for the corresponding coordinate basis vectors. 

Various typc;s of observers appear in the analyses be- 
low, including ones at rest in an Earth-centered frame, in 
a laboratory frame, in a satellite frame, and others. The 
corresponding frames are specified as needed in the sec- 
tions that follow. Among the sets of basis vectors having 
generic applicability are one for an observer at rest in the 
Sun-centered frame and a related one for an observer in 
uniform motion relative to the Sun-centered frame. We 
summarize these two sets briefly here. 

For an observer at rest at the point (T, X) in the Sun- 
centered frame, dX'' /dT = 0. Suitable basis vectors are 
denoted as with fi = {t,j), and they can be written 
as [7] 

et = 5^Jl + i/iTT(r,X)+PNO(4)]eT, 

ej = S'jiej - ^hj''{T,X)eK]+S-'jhTj{T,X)eT. 

(119) 
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This basis is orthonormal. 

If the observer is in motion with four-velocity in 
the Sun-centered frame, then an appropriate set of basis 
vectors can be taken as with /t = (t, j), where the 
components (e^)" are identified with the four- velocity, 
(e^)" = u~. This basis is given by [7] 

e£ = 6*^(1 + ^v'^)et + v^ej, 

e. = S'^^v'^Rkjet + S'^^iS''' + ^v''v')Rijek, (120) 

where is the coordinate velocity of the observer in 
the frame (119), and Rjk implements the appropriate 
rotation. This basis is also orthonormal. 

In applying the above equations, the relevant contribu- 
tions to the metric fluctuation and the observer velocity 
can be obtained from the modified Einstein equations and 
from the equation of motion of the observer. Note that 
the results typically depend on coefficients for Lorentz vi- 
olation. Also, some simplifying assumptions can usually 
be adopted without loss of generality. For example, in 
certain laboratory experiments the contributions to the 
metric fluctuation hsA sourced by the energy-momentum 
tensor of the Sun can safely be neglected. 

We remark in passing that the above definition of the 
Sun-centered frame could be sharpened in various ways, 
such as allowing for the precession and nutation of the 
Earth, establishing the vernal equinox via the centroids 
of bodies, and incorporating the motion of the Sun with 
respect to the center of the solar system. Some of these 
effects may allow additional sensitivities to Lorentz vio- 
lation via the resulting time dependence of the standard 
frame. Note also that the notion of parallelism used in 
the Minkowski-spacetimc definition of the Sun-centered 
frame is inapplicable in the context of curved spacetime. 
One way to address this latter issue is to define the Z 
axis so that aligns with the spin ax;is of the Earth after 
Eqs. (119) and (120) with Rjk = and the appropriate 
velocity are applied. For the various searches considered 
in this work, the standard definition of the Sun-centered 
frame suffices. A more complete investigation of these 
issues is of potential interest but lies beyond our present 
scope. 



B. Sensitivities 

In the following sections, we consider the observational 
effects of the coefficients (a^)^, {c^)fiu in the matter 
sector and s^i, in the pure-gravity sector. This subsection 
offers some comments about attainable sensitivities to 
these coefficients. 

Measurement of the coefficients (a^j)^ is of particular 
interest because they are virtually unexplored to date. 
The existence of the field redefinitions described in Sec. 
II C 1 means that observation of effects from {a^ff)^ re- 
quires either flavor-changing physics or gravitational cou- 
plings. At the level of quarks, the fiavor-changing weak 



interactions have been used to provide access to observ- 
ables involving differences of two coefficients {a"')^ with 
w including second- and third-generation quarks [53, 54]. 
Flavor oscillations can also be used to constrain the coef- 
ficients in the neutrino sector, where they form 3x3 
matrices in fiavor space [55]. However, to date gravi- 
tational couplings have been used to obtain sensitivity 
only to limited combinations of the 12 independent com- 
ponents of the SME coefficients (a*^)^, {oP)^, {o")tJ. for 
electrons, protons, and neutrons [11, 12]. These coef- 
ficients are otherwise unconstrained and could be com- 
paratively large, so they offer interesting prospects for 
further investigation in gravitational tests. 

In the present context, we can extend the single bound 
on (a™)^ given in Ref. [11] by taking advantage of the 
result of Sec. II C 1 that (a^)^ always appears at leading 
order with (e^)^ in the combination (a™g)^ given by Eq. 
(16). Using this result immediately yields a constraint 
on three of the independent components of {a^g)^ for 
electrons, protons, and neutrons, given as 

]a(a|ff)T+a(a^ff)T-0.8a«ff)T| < 1x10"" GeV (121) 

at the 90% confidence level. 

In contrast, many of the coefficients {c^)^i/ are read- 
ily observable in nongravitational experiments. Nonethe- 
less, gravitational tests offer additional opportunities to 
achieve sensitivities to (c"')^^, including some compo- 
nents that are unmeasured to date. For electrons, pro- 
tons, and neutrons, there are 27 independent observable 
symmetric coefficients {c^)^v A compilation of existing 
limits on {c^)fi„ for different flavors w is given in Ref. [2]. 

The coefficients s^^ lie in the pure-gravity sector 
of the minimal SME and therefore can be measured 
only in the gravitational context. The corresponding 
post-newtonian corrections to the gravitational field are 
known [7, 56]. Constraints on most of the nine in- 
dependent components of s^^ have been obtained us- 
ing a variety of techniques, including among others 
perihelion-precession studies, lunar laser ranging, atom- 
interferometer gravimetry, and laboratory and space- 
based experiments [4-7]. All these analyses disregard 
matter effects. In this work, we show that Lorentz viola- 
tion in the matter sector can contribute in different ways 
to signals involving the coefficients s^^. 

For all the coefficients (a™ff)j(i, (c^)^i/, Sfj,v, the effects 
of interest here involve gravitational couplings to matter. 
It is therefore reasonable to expect that the best sensitiv- 
ities to Lorentz violation arc associated with couplings to 
dominant gravitational effects. This suggests that tests 
with high sensitivity to Newton gravity are of particular 
interest. As described in Sec. IV B 3, the fiavor depen- 
dence of the coefficients for Lorentz violation implies that 
WEP tests also lie in this category. 

Many of the signals sought in gravity tests require an- 
cillary measurements of time and distance. These typ- 
ically involve matter in some form, and they may in- 
troduce additional Lorentz- violating effects beyond those 
comprising the direct signal of interest. However, most 
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of these additional effects are negligible in the present 
context because the corresponding coefficients are tightly 
constrained via tests in Minkowski spacetime [2] , whereas 
sensitivities in gravitational tests are typically substan- 
tially reduced by the weak gravitational field. Among 
the coefficients of interest in the present work, this is- 
sue is relevant only to {c^)^^ because {a^g)^ and s^^ 
are unobservable in Minkowski-spacetime tests and be- 
cause we adopt the coordinate choice (23) making unob- 
servable the photon-sector coefficients (fcF)"^ai/- Among 
the coefficients (c*")^^ for ordinary matter, the neutron- 
sector coefficients (c")^,y are the least well constrained 
at present. Their effects may therefore be important for 
certain tests, in which case a detailed analysis of the mea- 
surement method may be necessary. 

Another consideration relevant for identifying sensitiv- 
ities in tests with atoms or bulk matter is the role of the 
contributions from binding energy. In some cases, ac- 
counting for these contributions can disentangle eff'ects 
from different coefficients, thereby producing additional 
independent sensitivities. This can occur when coeffi- 
cients from two or more sectors are involved, either di- 
rectly within a WEP test or indirectly via comparison 
of results obtained for different bodies. In the remain- 
der of this subsection, we discuss this possibility for the 
coefficients (a^)^ and (c^)/ii/ in turn. 

Consider first combinations of the coefficients {a'c«)t^- 
Following the discussion in Sec. IV A 2, a body B has 
an effective coefficient {afg)^ given by Eqs. (73) and 
(75). The dimensionless quantity relevant for a test is 
(aj^g )p/m^, and comparisons involving two bodies there- 
fore appear as the difference of two quantities of this 
form. For two neutral bodies involving bound electrons, 
protons, and neutrons, this diff'erence can be expanded 
as follows: 

m\m2 

Here, the numbers of particles of species w for the two 
bodies are , N2, and m'l, m'2 are the binding-energy 
contributions to the masses mi, m2 of the two bodies, as 
defined in Eq. (71). Also, we define 



!' 



V"eff 



(a^ff)^. (123) 



When the contributions from binding energy are ne- 
glected in Eq. (122), the linear combination (a^^^"")^ of 
coefficients becomes the sole observable involving {0^0)^1 
in a comparison of two bodies, with the effect scaled 



by their difference in species content. However, incor- 
porating the binding energy in the analysis introduces 
the last two terms in Eq. (122), revealing that the effects 
of (aeff)^ + (^eff)M a'^d (Sgg)^ vary differently with the 
content of the bodies. This allows the possibility of in- 
dependent measurements of {al.ff)n + (S^g)^ and {a\)^- 
Note that the sensitivity of such measurements is typi- 
cally an order of magnitude less than that of measure- 
ments of (floff'' ")/:* to the appearance of ratios of 
the form m! jm. 

Next, consider combinations of the coefficients (c™)^;,. 
For a body B, the effective coefficient (c^);^^ is a dimen- 
sionless quantity given by Eq. (74). As discussed in Sec. 
IV A 2, nonzero Lorentz- violating contributions from the 
binding energy given by the coefficients (c'^)^^ are ex- 
pected to exist, along with the usual binding-energy con- 
tributions ml to the body mass. It turns out that these 
Lorentz-violating contributions impede the use of bind- 
ing energy to extract additional independent sensitivities 
to combinations of the coefficients (c"')^,y. To see this, 
consider two neutral bodies as before, and expand the 
analogue of Eq. (122) to get 
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where we introduce 



TD ^ —I— Tn^ 

{r+P--),. = (r+^)^,-_^(c")^.. (125) 



When binding- energy effects are neglected, (c*^"*"^"")^^ 
becomes the only observable combination of the coeffi- 
cients {c^)nv in gravitational tests comparing two bod- 
ies. Including the binding- energy terms as in Eq. (124) 
shows that bodies with different species content can ex- 
hibit distinct effects. Although it seems unlikely that 
nonzero effects at order m'/m in a variety of bodies would 
cancel sufficiently well to evade detection altogether, the 
appearance of the unknown coefficients {c'^)ij,v makes it 
infeasible at present to extract unambiguous independent 
measurements on combinations of the coefficients {c^)ij,v 
using binding-energy effects. 
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VII. LABORATORY TESTS 

This section considers some sensitive laboratory tests 
with ordinary iicnitral bulk matter, neutral atoms, and 
neutrons performed on or near the surface of the Earth. 
The basic theory for these tests is developed in Sec. VII A, 
while Sees. VII B, VII C, VII D, and VII E consider sig- 
nals and sensitivities attainable in a variety of terrestrial 
searches. More exotic laboratory tests with charged par- 
ticles, antimatter, and particles beyond the first genera- 
tion arc considered in Sec. IX. 

Terrestrial experiments seeking gravitational Lorentz 
violation using ordinary matter can be classified either 
as gravimeter tests or as WEP tests. In gravimeter tests, 
the basic idea is to seek variations either in the gravita- 
tional force on a test body or in its gravitational accel- 
eration. The corresponding signals originate in the time 
dependence of laboratory coefficients for Lorentz viola- 
tion induced by the rotation of the apparatus and the 
rotation and revolution of the Earth. These signals can 
be interpreted as an effective time variation of the New- 
ton gravitational constant Gn- In WEP tests, the idea is 
to compare either the gravitational force between two dif- 
ferent bodies or their relative gravitational acceleration. 
The corresponding signals, which can be instantaneous 
or time-varying, are sensitive to differences between the 
coefficients associated with different species of matter. 

Lorentz violation can introduce deviations from New- 
ton's second law, so the distinction between force and 
acceleration can be important. This distinction implies 
the two classes of gravimeter and WEP tests can each be 
further subdivided into two categories, force-comparison 
tests and free-fall tests. The basic idea of a free-fall test 
is to search for a time or composition dependence in the 
gravitational acceleration of a freely falling test body by 
monitoring its motion. The idea of a force-comparison 
test is to balance the gravitational force experienced by 
a test body with a second force, investigating changes in 
the equilibrium arising from the time or species depen- 
dence of the laboratory coefficients for Lorentz violation. 
The force comparison can be achieved either by using 
a seesaw arrangement to balance the gravitational forces 
on test bodies of different composition, which constitutes 
a force-comparison WEP test, or by using a nongravi- 
tational force to counter the gravitational force on the 
test body, which represents a force-comparison gravime- 
ter test. 

We thus have four categories of possible laboratory 
tests with ordinary matter. In what follows, each is 
considered in a separate subsection. Free-fall gravime- 
ter tests, including searches with freely falling corner 
cubes and with atom interferometers, are considered in 
Sec. VII B. Force-comparison gravimeter tests using me- 
chanical and superconducting gravimeters are studied in 
VII C. Free-fall WEP tests, which come in a wide vari- 
ety of forms, are considered in Sec. VII D. Finally, force- 
comparison WEP tests are discussed in Sec. VII E, with 
focus on a torsion-balance configuration. 



Table I provides a list of some conventions adopted in 
this section for the analyses of laboratory tests. Many of 
the quantities are self-explanatory. The laboratory speed 
Vl is due to the rotation of the Earth and depends on the 
laboratory colatitude %. The relative time involves 
a convenient choice of origin, measured from any instant 
when the y axis in the laboratory frame and the Y axis of 
the Sun-centered frame coincide. Using instead of the 
canonical time T in the Sun-centered frame introduces a 
phase tp in the analysis. The angle ( is defined in terms of 
the component accelerations a^, a.^ of a test body along 
X, z in the laboratory frame, 

C = tan-^(as/a£). (126) 

At leading order, ( is approximated by the ratio of the 
usual Newton centripetal and gravitational accelerations, 
1^ ~ 10~^. It represents the angular deviation from the 
vertical at the location of the laboratory of a plumb line 
or of a test body in free fall. 

Table I. Notation for laboratory tests. 



Quantity Definition 

R mean Earth-Sun distance 

i?0 mean Earth radius 

O mean Earth orbital frequency 

oj mean sidereal frequency 

uje apparatus rotation frequency 

V© = ni? mean Earth orbital speed 

Vl laboratory rotational speed 

relative time 

T] inclination of Earth orbit 

X laboratory colatitude 

tp — w(T0 — T) phase induced by 



C ~ f^^-Re / sin(2x)5 deviation angle 



A. Theory 

The relevant observables for laboratory tests of Lorentz 
symmetry in gravity are the motions of test bodies rel- 
ative to the Earth and relative to each other. These 
observables can be obtained from the action for a test 
body, evaluated at the appropriate post-newtonian order 
and expressed in laboratory coordinates. 

Consider the action for a test body given in Eq. 
(76), with the gravitational field of the Earth acting as 

(3) 

the source S. The corresponding lagrangian La_c describ- 
ing the motion of the test body T at PN0(3) can be 
constructed by expanding 5„ with B=T. The solution 
for the metric fluctuation /i^^ at this order is obtained 
from the general expressions (81), (84), and (88), with 
the Earth treated as a rigid rotating source S as described 
in Sec. IV A 2. In what follows, we neglect the gravita- 
tional fields of other bodies such as the Sun, although in 
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a more detailed treatment these could be incorporated 
using similar methods. 

For laboratory searches, it is convenient to begin cal- 
culations in an Earth-centered frame with coordinates 
denoted by x'' = {i, x, y, z). At leading order, the spatial 
components of the Earth-centered basis are taken to co- 
incide with those of the Sun-centered frame, and i = T. 
In the Earth-centered frame, we find 
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(127) 



where r = ^^xjxj. 



This expression contains the conven- 



tional Newton kinetic and potential terms for a test body 
T moving in the gravitational field of S, along with a se- 
ries of corrections that depend on the coefficients {a^ff)ji 
and (c^)fiv. Some of these additional terms are motional, 
analogous to centrifugal effects, and some are gravita- 
tional, including ones analogous to gravitomagnetic ef- 
fects. Effects from the Earth's motion about the Sun 
arc implicitly included via the dependence of (a^fj)^ and 
(c^)/ii> on the orbital speed V®. This dependence can 
be made explicit by expressing the coefficients in Sun- 
centered coordinates instead of Earth-centered ones. 

To obtain results applicable to laboratory tests, the re- 
sult (127) must be transformed from the Earth-centered 
frame to the laboratory frame. We denote the laboratory 
coordinates by x^, where the spatial coordinates x^ are 
taken to coincide with the standard SME conventions for 
a laboratory on the surface of the Earth [21]. In the lab- 
oratory, the X axis points South, the y axis points East, 
and the z axis points towards the local zenith. To the 
required post-newtonian order, i = t and the coordinate 
location of the laboratory in the Earth-centered frame 
can be written [7] 

i=RQ (sin X cos(a;eT -|- (f)) , sin x sin(a;eT + (f)),cosx)- 

(128) 

The transformation between the two sets of spatial coor- 
dinates can therefore be written 



^5 



a + R--x^, 



(129) 



where i?; ■ is the relevant rotation between the bases of 



the laboratory and the Earth-centered frames. Note that 
Eq. (129) implies the coefficients (a^g)/! and {c^)p,o in 



the laboratory frame acquire implicit dependences on the 
laboratory speed and on the sidereal frequency oj, 
which arise from the rotation of the Earth. 

The inclusion of the Earth's rotation in the analysis 
implies the laboratory frame is noninertial. The struc- 
ture of the first few terms in Eq. (127) reveals that iner- 
tial forces in the laboratory couple to {c'^)fii>, which can 
result in nongravitational Lorentz-violating effects com- 
parable in size to the gravitational ones of interest. We 
therefore incorporate these nongravitational effects in our 
subsequent analyses. In practice, this means effects pro- 
portional to the centrifugal acceleration lo'^Rq ~ 10"^^ 
must be considered. 

In what follows, we consider effects up to and including 
PN0(3). The leading PN0(3) effects are proportional to 
the speed of the Earth as it revolves about the Sun 
and are of order gV^ « 10"^^, where g = Gntii^ / R% for 
a laboratory on the surface of the Earth. This yields sen- 
sitivity to various components of the coefficients (a^)/n 
and {c^)fiv For some laboratory tests, it is advanta- 
geous to consider also PNO (3) effects proportional to the 
smaller speed Vl of the laboratory due to the rotation of 
the Earth, which are of order gVL ~ lQ~^g. The benefit 
arises in two ways. First, inclusion of the boost Vl intro- 
duces effects proportional to [a^g)^ that vary sidereally 
instead of annually. This offers access to ia^g)f_i for mea- 
surements conducted on comparatively short time scales, 
albeit at a sensitivity reduced by about two orders of 
magnitude. Second, certain laboratory tests have greater 
sensitivity to forces in the x and y directions than to ones 
in the z direction. The inclusion of effects from (a^)/n 
and {c'")^iv that are proportional to Vl can then intro- 
duce new sensitivities or improve existing ones. 

So far, modifications to the trajectory of the test body 
arising from the coefficients s^^ have been disregarded. 
However, it is straightforward to incorporate these in 
the lagrangian at PN0(2) because the coordinate choices 
made here are consistent with those of Ref. [7] at this 
perturbative order. In the laboratory frame, we find the 
PNO (2) contribution from 's^i, to the lagrangian of the 
test body can be written 

= yn^g(^.^^j^-s-~y-l-s^^z-lsfiz). (130) 

(2) 

It turns out that Ll suffices to achieve sensitivity to 
Sf,^ at PN0(3). The point is that the leading PN0(3) 
effects are proportional to P®, while inclusion of effects 
proportional to Vl offers no additional benefit in this 
case for the tests we consider. The coefficients s^^ are 
species independent, so they are unobservable in WEP 
tests. Moreover, inspection of Ll reveals that the coef- 
ficients sp,p already vary at the sidereal frequency through 
the transformation to the Sun frame. 

In the laboratory frame, the PNO (3) lag rangian L/a,c,s 
obtained from Eq. (127) and incorporating effects from 
Sfic, via Eq. (130) is somewhat lengthy in form. As an 
illustration of its structure and implications, we can re- 
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strict attention to its PN0(2) limit L^^c,s- We find 



T 

-m g 



2a 



2a 



(^eff)t + -s(«eff)t 

III m 



(131) 



Varying this result yields the Euler-Lagrange equations of 
motion, which we can express in the form of the modified 
force law 



B. Free-fall gravimeter tests 

In this subsection, we consider laboratory tests that 
monitor the motion of a test body in free fall near the 
surface of the Earth. The equation of motion for the 
test body can be obtained from the PN0(3) lagrangian 

(3) 

La.Cs described in Sec. VII A. Its explicit form is lengthy. 
However, all information relevant for present purposes is 
contained in its solution expressed to the desired pertur- 
bative order. This solution can be written in the form 
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(135) 



At this perturbative order, the inertial and gravitational 
forces acting on the test particle are given by 



Fs = 



m^gszx, 
rn^gszy, 

T 

-m g 



2a 
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+ (c^)tt + (c'^)ft + l%-+^s- 
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m 



^fe = m^(l + (c^),-,05-.^ + 2m^(c^)(3.^) 



(133) 



(134) 



where the test body has initial position Xo and initial ve- 
locity Vg. The quantities of interest in searches for grav- 
itational Lorentz violation are the components a- of the 
acceleration of the test body in laboratory coordinates. 

For purposes of data analysis and reporting sensitiv- 
ities to coefficients for Lorentz violation, it is useful to 
express the acceleration components a- in a form that 
displays explicitly the time variation and the dependence 
on particle species. In free-fall gravimeter tests, the time 
variation appears at frequencies 0, lu, 2uj, lu ±Q, 2lu 
and n, which are collectively labeled as n in what follows. 
The dependence on particle species arises from the com- 
position of the test and source bodies. It is characterized 
by the label w, which ranges over e, p, n for ordinary 
matter. 



Table II. Amplitudes for the acceleration a^ 



Amplitude 



is the effective inertial mass. 

These results reveal the generic feature that the grav- 
itational force Fj acquires tiny corrections both along 
the z direction and perpendicular to it. Also, the re- 
sponse of the test body deviates slightly from the di- 
rection of the applied force because the effective inertial 
mass m--^ depends on the coefficients {c'^)ni>. In princi- 
ple, some of these effects are detectable in sensitive labo- 
ratory tests, and the corresponding signals are discussed 
using PN0(3) results in the following subsections. 

Some coefficients appear in combinations that are chal- 
lenging to separate in laboratory tests. This is true, 
for example, of the coefficients a(aJg)T and {c^)tt- = ^"T'"'((c"')jfx — (c"')rr) sin2x 



Phase 



^0 



m-^ smxcosx[{c^)xx + (c"')rr - 2{c^)zz] 



AZ = 2m»(c-)(xz) cos2x+ fVLa(a-ff)ycosx 
= \VL[a{a^s)Y + 2m-(c-)(Ty)] cosx 
2m-^{l^)^YZ) cos2x - lVLa{a^s)x cosx 
= -IVl Ha^^)x + 2m-'{^\TX)] cosx 



-- 



Consider for simplicity the scenario with only isotropic 
Lorentz violation in the Sun-centered frame, where 



2ui 



m-(c-)(xy) sin2x 



the nonzero coefficients are a{aJg)T and 



3(cT)x;f = 3(cT)^^ = 3{c^)zz. 



TT = 

In the laboratory 

frame, (a^)(- « (^(a^)^ and ~ (c-'^)tt up to boost 

factors. These coefficients therefore cannot be readily 
separated in gravimeter tests, which depend on time vari- 
ations from anisotropic effects. Moreover, inspection of 
the PN0(2) lagrangian (131) reveals that if 3a(ajg)( = 
m^{^)^ then the contributions of a(a^)j and to 
the effective inertial and gravitational masses are identi- 
cal. The combination Q:(a^)T — m'^{c'^)TT/^ therefore 
cannot be readily separated in conventional WEP tests 
either. Note that WEP tests comparing a particle and its 
antiparticle can in principle evade this difficulty because 
the sign of a{aJg)T differs between the two. Another pos- 
sibility would be to compare matter with light, an option 
considered further in Sec. XL 



K+n = -m^^V^icniTX) sinr?cos2x 
^:^+o = -^"V-©[(5")(Ty)sinr? 

- (c'")(TZ)(l - COST?)] cos2x 
^"-n = -m^V(s(c^)^TX) sinr,cos2x 
S"-a = -™'"^^©[(c"')(Ty) sin 7? 

+ (c"')(TZ)(l + cosr;)] cos2x 
-im"'l/e(c'")(Tr)(l - cosr?)sin2x 



b: 
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im»Fe(c-)(Tjf)(l-cosr,) sin2x 
^L-n = im'-ye(c'")(Ty)(l + COST?) sin 2x 
BZ-a = -im-ye(c'")(TX)(l +cosr,)sin2x 



= -m^V(s[{c^)^TY) COST? 

-2(c'")(Tz) sinr?] sin2x 
= m-ye(c'")(TX) sin2x 
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Table III. Amplitudes for the acceleration a.y 



Amplitude Phase 



Co"=m-yi(c-)(Tz)sin2x 
(7-=2m-(c-)(y2) cosx 

- |VLa(a^ff)x + 2m'"FL(c™)(TX) sin^ x V 
C'J = -^Vl [aia:s)x + 2(c'")(tx)] 
= -2rn-(c"')(xz) cosx 

+ 2m-VL(c'")(Tr) sin^ x - ^VLa{a^fi)Y V 

= -i^i [2m-(c-)(^^) + a(a-ff)r] V 

C2»^ = 2m-(c'")(;fy)sinx 2^ 

= -m»((c"')x;f - (cnYY) sinx 2^ 
C;:r+a = m^Ve [(c"')(tz)(1 - cosr?) 

- (c"')(Ty) sin??] cosx V' 

^"+o = (c"')(T;f ) sin cos X V 

C^clT-n = -m'"y® [(C'")(TZ)(1 + COST?) 

+ (c"')(Tr) sin r?] cosx V' 

£'"-n = m-ye(c")(TX) sin t? cosx V' 

C^2"^+a = m'"ye(5")(TX)(l - cosry)sinx 2^^ 

^L+n = m-Ve(c'")(TY)(l - cosr,)sinx 2^^ 

<^2" -n = -"^"^e (c"')(TX) (1 + cos rj) sin x 2^ 

^L-n = -m-ye(2")(Tr)(l + cost?) sinx 2^ 



For the .i component of the acceleration, some calcu- 
lation yields an expression of the form 

as = w^i?© sinx cosx 

r / AT"" A^"" \ 

[ [^"^n + + COsKr + V'n) 

+ [-^^B- + + iS„ j sinKT + V>J . 

(136) 

In this equation, the amplitudes A^, B™, _B^"' con- 
tain the coefficients for Lorentz violation {a^g)^, (c™)^!/ 
and hence depend on particle species. These amplitudes 
and their associated phases arc listed in Table II. The re- 
maining amplitudes j4„, _B„ contain the coefficients s^i^ 
from the gravitational sector, which are independent of 
the composition of the test body. These amplitudes can 
be obtained from the amplitudes A"^, by the sub- 
stitutions to'" -> 1, (c'")ee ^ssH, and (a^)H 0, 
disregarding contributions proportional to Vl ■ 

The ij component of the acceleration can be decom- 
posed similarly. We find 

= E4(;^^» +;#^" + 5Cn]cOsKT + V'n) 

(137) 



The amplitudes C;f, C^™, D^^, depend on particle 
species through the coefficients (aj|^)^, (c"')^^ and are 
listed in Table III, along with the corresponding phases. 
The remaining amplitudes C„, JD„ are obtained from C™, 
using the substitutions to*" — >■ 1, (c'"')i;h — ^seh, 
and (a™g)E: — > 0, disregarding contributions proportional 
to Vl as before. 

For the z component of the acceleration, we obtain 



aj = -g + upR(sSiT?x 

r / TV™ iV'"' \ 

+ E 4 ( ;;^T + + ) cos(a;„T + V^„) 

■^Fn + + ^Fn) sin(a;„r + V'^) . 

(138) 

The amphtudes i;^, f;^, F^"" depend on particle 
species via the coefficients (ajf^)^, (c™)^^. The ampli- 
tudes En, Fn arc independent of species and given in 
terms of the coefficients s^,y. All these amplitudes and 
their phases are provided in Table IV. 

In principle, the results of a free-fall laboratory test 
using any gravimeter can be analyzed with the above 
equations. The dominant effects appear at different fre- 
quencies for different coefficients, so the time scale of data 
taking in a given experiment affects the breadth of its 
reach in coefficient space. Also, each signal frequency can 
be expected to have distinct systematics. For example, 
dominant effects from the coefficients {acs)j oc;cur at the 
annual frequency fi, for which seasonal systematics are 
relevant. Note that all the Lorentz- violating effects can 
be accessed at or near the sidereal frequency lj, although 
in some cases at reduced sensitivity. 

At least two kinds of devices can be classified as free- 
fall gravimeters: falling corner cubes, and matter inter- 
ferometers. Falling corner cubes, which typically are 
sensitive only to the direction of the free-fall motion, 
are used to monitor time variations of the gravitational 
field for geodesy and other geophysical purposes [57]. 
In principle, they are of interest for free-fall gravimeter 
tests of Lorentz violation. However, matter interferome- 
ters presently carry several advantages over falling corner 
cubes in this context. They are slightly more sensitive, 
some types can sense accelerations in more than one di- 
rection, and the composition of the test body can be 
determined more readily. We therefore focus on matter 
interferometers in this subsection, revisiting the use of 
both falling corner cubes and interferometers in the con- 
text of free-fall WEP tests in Sec. VII D. 
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Table IV. Amplitudes for the acceleration aj. 



Amplitude Phase 



= -2a{a^g)T + 2m^{c^)zz cos^ x 

+ {{c^)xx + {c^)yy) sin^X 

E'o^ = -2a(a™ff)T - m^ic^)TT 

= 2m-(c-)(xz) sin2x - IVLa(a^^)Y sinx V 

= -IVl (3a(a- )r + m-(c-)(Ty)) sinx V 

= 2m-(c"')(yz) sin2x + |Via(a™g)x sinx ^ 

= IVl {Ma:s)x + m"'(c-)(Tx)) sinx V 

i;^^ = m»((c-)xx - (cnYY) sin2 x 2^ 

FZ = 2m^{c^)^XY)Sin^X 2^ 

^"+n = -m-y©(c-)(Tx) sin77sin2x V 
^"+n = -^"^e[(c'")(Tr) sinr? 

- (c'" ) (TZ) ( 1 - COS ??)] sin 2x V 

K-n = -m"Ve(c'")(TX) sinr?sin2x V 
F:^_^ = -m^V^[{c^\TY)smrj 

+ (c"')(TZ)(l + COST])] sin 2x V 

i^2" +n = -m"'Ve{cn(TY){l - cosr,) sin^ x 2^ 

^^2^+0 = m™l^e(S")(T;f)(l - cost?) sin^ x 2^ 

^^2" -Q = m-ye(2")(Tn(l + COST?) sin^ x 2^- 

^^2'"^-o = -m«'T/e(^"')(rA-)(l + cosr?)sin2x 2^ 
= 2Fea((a^a.)y cosr/ + (a-.)z smr?) 
-2m'"ye[(c"')(TF) cosTysin^x 

+ 2(c^)(TZ)sin?7COs2x] 
E'^ = 2yea((a™ff)y cosr? + sin,?) 

+ 2m-ye((c"')(ry) cost? + (c"')(tz) sinr/) 

= -2Fea(a-ff)x + 2m™l/e(c™)(T;f) sin^ x 

F/r = -2Fea(a-ff)x - 2m^V^{c^\TX) 
Fo = - \szz cos^ X - i (sxx + syy) s\v? x 

- fsTT 

-Ea, = -lS(xz)Sin2x ip 

^l, = -5S(FZ)Sin2x V 

-^2^ = - i (sxjf - Syy) sin^ x 2^) 

^2.^ = -iS(xy)Sin^X 2^) 

Eui+o. = \ t^e S(rx) sin r? sin 2x '^i' 
F^+n = jV(s[s{TY) sin?? 

- S(TZ)(1 - COST?)] sin 2x V 

i^i^-n = i Ve S(TX) sin ?? sin 2x V 
^lu-n = [s(Ty) sin?? 

+ S(TZ)(1 + COS?/)] sin2x V 

-E2a;+a = 3V®S(Ty)(l - COS??) sin^ x 2^) 

^2c^+n = -il^©S(TX)(l - cos??) sin^ x 2?/^ 

-E2w-n = -il^©S(Ty)(l + cos??)sin^X 2^) 

-F2u-n = 3'^©S(TX)(1 + cos??) sin^ x 2^) 
-Ba = [s(Ty) cos ?7(i sin^ X + 3) 

+ S(TZ) sin ??(cos2 X + 3)] 

Fa = -yeS(TX) (5 sin' X + 3) 



Matter interferometers, which permit quantum- 
mechanical laboratory measurements of the motion of 
falling matter, have attained impressive sensitivities to 
gravitational acceleration [58] and to rotational accelera- 
tions via the Sagnac effect [59] . In the context of gravita- 
tional Lorentz violation, matter intcrferomctry has been 
used to measure combinations of the coefficients s^i, and 
(^F)"^ai/ [5] based on the gravimeter analysis of effects 
from the pure-gravity sector of the SME [7]. Here, we 
extend the latter analysis to include effects from the co- 
efficients (a™g)^, (c'")^y and generalize it to other inter- 
ferometer configurations. 

The basic idea of a matter interferometer is to place the 
matter in a superposition of spatially separated quantum 
states, which may acquire a measurable relative phase. 
In the gravitational tests considered here, the behavior 
of the interferometer is close to the classical limit, and 
a convenient way to perform the analysis is to proceed 
semiclassically via path integration along the classical 
motion [58, 60]. The phase difference between the final 
states can then be viewed as a sum of three contribu- 
tions: the phase difference acquired from the momentum 
transfers used to control the beams, the phase difference 
accumulated from the classical action along the different 
paths, and in some configurations a phase difference com- 
ing from a final separation of the states. It turns out that 
the dominant effect for acceleration sensing is the phase 
difference acquired through momentum transfers. Since 
leading-order Lorentz-violating motional effects appear 
as modified accelerations, the phase difference from the 
momentum transfers is the relevant contribution in the 
present context. 

For definiteness, suppose the interferometer paths 
trace a parallelogram. This includes the limiting case of 
temporal path separation, where the parallelogram has 
zero area. Other shapes could also be analyzed using the 
equation of motion (135). In the Lorentz-invariant case, 
the standard result for the phase shift due to the Earth's 
gravitational field is Acp = k^gr^, where k is the magni- 
tude of the momentum transfer in the beam splitter and 
r is the time of flight between impulses. With Lorentz 
violation present, we find the phase shift A(j) takes a sim- 
ilar form but with the Newton gravitational acceleration 
replaced by the accelerations in Eq. (135), giving 

A^ = k.a-^T^. (f39) 

The signal frequencies associated with Lorentz viola- 
tion can be identified by substitution of the expressions 
(136)-(138) for the acceleration components a--. Note 
that Lorentz-violating effects on the atomic energy lev- 
els could generate additional contributions to the phase 
difference but are already tightly constrained in other ex- 
periments and so can typically be neglected. Note also 
that possible Lorentz-violating effects varying with the 
particle spins, which are described explicitly by the rel- 
ativistic hamiltonian of Sec. IIIB, are disregarded here 
as outside our present scope. A comprehensive investiga- 
tion of their implications for matter interferometry may 
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be of interest [61]. 

Several atom interferometers currently or recently op- 
erating are relevant to free-fall gravimeter searches for 
Lorentz violation. An impressive sensitivity of about 
1 X lO^^"^ to the vertical acceleration was achieved by 
Peters, Chung, and Chu [58]. In another apparatus, a 
differential-acceleration sensitivity of 3 x 10~^g/V^ has 
been demonstrated [62]. An interferometer designed for 
experiments in space is expected to achieve sensitivity 
of about 3 X 10~^g in ground operations [63]. Initial 
sensitivities to accelerations in each direction of about 
6 X lO^^g after 10 minutes of averaging have been at- 
tained in a device using highly parabolic trajectories 
[64] . Recent estimates suggest that future measurements 
of vertical acceleration could achieve sensitivities at the 
level of about IQ-^^g [65]. 

Given this information and the phase shift (139), we 
can use Tables II, III, and IV to obtain crude estimates 
for attainable sensitivities to coefficients for Lorentz vi- 
olation in existing or near-future atom interferometers. 
With present capabilities, sensitivities at the level of 
parts in 10^ could in principle be obtained to combi- 
nations of the coefficients Q:(a^fj),7, J = X,Y,Z and of 
several currently unconstrained components of the coef- 
ficients (c^)eh, including (c")(7 j) for the neutron. The 
relevant signals are associated with the Earth's boost as it 
revolves about the Sun, so they exhibit an annual period- 
icity. The next generation of atom interferometers could 
in principle improve this sensitivity to parts in 10^*^. The 
boost of the laboratory due to the Earth's rotation pro- 
vides sensitivities that have sidereal periodicities instead 
but that are weaker by a factor of about 100. Note that 
the boost suppressions could in principle be avoided for 
certain coefficients, including presently unboimded com- 
binations involving (c")zz, by the use of an interferome- 
ter sensitive to the accelerations a^, a,y that is placed on 
a rotating turntable. Note also that individual sensitivi- 
ties to neutron coefficients can in principle be extracted 
by performing atom interferometry with different neutral 
atoms having distinct proton-to-neutron ratios. Another 
possibility with weaker existing sensitivity includes neu- 
tron interferometry [66], which could provide indepen- 
dent and clean bounds on neutron coefficients. 



given force-comparison gravimeter can be extracted from 
the PN0(3) lagrangian ii'^c.s discussed in Sec. VII A. 
Since macroscopic bodies are involved, the analysis must 
include an assessment of their composition. Note also 
that conventional intuition from Newton's second law 
can be misleading because the effective inertial masses 
depend on the coefficients (c^)^i/, as discussed following 
Eq. (134). 

Superconducting gravimeters have already been pro- 
posed as suitable devices for measuring the gravity-sector 
coefficients s^i, for Lorentz violation [7]. Here, we ex- 
tend this discussion to include effects from the; coefficients 
(a^)^ and (c™)^!,. The analysis proceeds directly from 

the PN0(3) lagrangian La,c,s by noting that the device is 
designed to maintain = 0. The applied force required 
to hold this constraint can be taken as the relevant ob- 
servable and can be written 



Fb' = F£CosC + -FxSinC, 



(140) 



where ( is the deviation angle defined in Eq. (126). To 

maintain consistent counting of small effects, we restrict 
terms independent of velocity to first order in C, and terms 
containing a power of velocity to zeroth order in (. 

Wc find that the relevant contributions to the force F^' 
can be decomposed by frequency as 

Fb' = mTg(l-Ctanx-iO-m'^5 



1\JW ATW 

+ -^E':' + iG„ ) cos(a;„T + V„) 



N'" 



m' 



TTV 



s " 



+ ii/„)sin(a;„r + V„) 



(141) 



where constant effects that are unobservable in supercon- 
ducting gravimeters are neglected. In this expression, the 
amplitudes GJ" , and their phases are given in Table 
V, while E'^ and Fl^"" are listed in Table IV. The remain- 
ing amplitudes G„ and ff„ can be expressed in terms of 
amplitudes given in Tables II and IV as 



C. Force-comparison gravimeter tests 



Gn = ^nC + En 



Hn=Br,C + Fr,. (142) 



Another class of gravimeter tests is based on the idea 
of countering the gravitational force with an appropri- 
ate electromagnetic force. Force-comparison gravimeter 
tests can be performed with gravimeters based on sys- 
tems of springs and masses [57] and with superconduct- 
ing gravimeters [57, 67, 68]. At present, the latter de- 
vices have sensitivities competitive with those of exist- 
ing atom interferometers. Certain experiments studying 
short-range gravity may also offer relevant sensitivities 
[69]. 

The signals for gravitational Lorentz violation in a 



The frequency decomposition (141) can be examined 

to extract crude estimates of attainable sensitivities to 
Lorentz violation. In this way, we estimate that the 
presently unbounded coeflacients Q!(a^)j and {c^)(tj) 
could be measured at the level of parts in 10^ using ex- 
isting data from superconducting gravimeters [68]. Im- 
proved sensitivities are likely to be attainable in a dedi- 
cated experiment of this type. 
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Table V. Amplitudes for the force F^i. 



Amplitude 



Phase 



GZ = 2m^acnixz) 



\VLaia:s)Y sinx - 2m"'VL(c-)(TF) sinx ^ 



if- = 2m-C(5")(yz) 



\VLaia:f,)x sinx + 2m"'yL(c"')(Tx) sinx ^ 



IXX 



H. 



2w 



2m-C(c'")(XY) 



GZ = 2V^a{{a:s)Y cosr? + (a-ff)z sinjj) 

+ 2m^VQ[{c^)(^TY) cosr/sin^x 

+ 2(c"')(Tz) sinrycos^ x] 
= 2V^a{a:^i,)x + 2m^V^ic^)(TX) sin^ X 



2V; 
2^} 



For tests of short-range gravity and certain other ap- 
plications, it is useful to consider the standard case of 
two point masses mi and m2 at coordinate locations xi 
and X2- With this setup, the modified Newton potential 
V at PN0(2) in the laboratory frame can be obtained 
from L^Xs- We find 

Gjvmimz r 2a ^ 2a ^ ^ 

V = — =rTlH («eff)t"l (aeff)t 

\xi - X2\ '- mi ^ ""^'^ m2 

where x = (xi — X2)/\xi — X2\- This modified poten- 
tial exhibits the usual inverse-distance dependence, and 
it generalizes Eq. (137) of Ref. [7]. The corresponding 
modified Newton force typically has a component per- 
pendicular to the unit vector x, while obtaining the ac- 
celerations requires determining also the effective iner- 
tial masses. As usual, any motion of the; masses relative 
to the Sun-centered frame implies time dependence of 
the laboratory- frame coefficients. In principle, the above 
modified Newton potential could be used in conjunction 
with integration or finite-element methods to determine 
the effects of the coefficients (a^)/:^, {c^)nv, Sfj,v on the 
behavior of two interacting bodies. 



D. Free-fall WEP tests 

In this subsection, we consider WEP tests in which 
signals for Lorcntz violation can be sought by monitoring 
the relative motion of two freely falling bodies of different 
composition. Typical free-fall WEP tests are sensitive to 
motion along the direction of the net acceleration a,^i. 
This acceleration is the combination 



&z' = Siz cos C + a* sin ( 



(144) 



of the component accelerations (136) and (138), weighted 
by the deviation angle ( given in Eq. (126). In what 
follows, terms containing both a boost factor and a factor 
of C are treated as higher order and negligible, as in the 
previous subsection. 



The relevant observable for free-fall WEP tests is the 
relative position Az' of two test bodies 1 and 2 in a given 
drop. It can be written as 

Az' = [{vo)l, - {vo)%] i + i (4, - a?,) P. (145) 

This relative position varies with the canonical time T. 
Decomposing by frequency yields the expression 



Az' 



^ \ mi m,2 



X [7^ cos(a;„T + V„) + J^T sin(a;„r + 



(146) 

The amplitudes and can be expressed as 

= A-C + El, J- = B-C + (147) 

where A^, are listed in Table II and E'l, F,^ are 
given in Table IV, along with the associated phases t\}n- 
In Eq. (146), the quantities iV™ and are the numbers 
of particles of type w appearing in the test bodies 1 and 
2, respectively, while mi and m2 are the corresponding 
conventional masses. 

The frequency decomposition (146) of the signal (145) 
can be used to provide rough estimates of attainable sen- 
sitivities to Lorentz violation in existing or near-future 
free-fall WEP tests. We combine values for the fractional 
acceleration sensitivity Aa/a discussed in Refs. [65, 70- 
75] with the result (146) to compile some estimates in Ta- 
ble VI. In this table, the first row lists the fractional accel- 
eration sensitivity, while each of the other rows concerns 
a particular combination of coefficients. For brevity, in 
the first column we adopt the notations 

(aeff )f-i-z = (ocff)y cos r/ -I- (aoff)z sin 77, 

(c-)q = {T^)xx + {T^)yy-2{c^)zz. (148) 

along with those introduced in Eqs. (123) and (125). We 
follow common procedure in the literature [2] by taking 
(c^)tt, (c'")jfx — (c^)ff, and (c'")q as the relevant inde- 
pendent combinations of the traceless coefficients (c'")se. 

Each column in Table VI lists estimated attainable sen- 
sitivities on the moduli of various quantities in specified 
types of free-fall WEP test, expressed to the nearest or- 
der of magnitude. Values listed with neither brackets 
nor braces are limits based on published data that are 
implied by our present analysis. Values shown in brack- 
ets arc our estimate of sensitivities that could in principle 
be obtained from a suitable reanalysis of existing data. 
Values shown in braces represent our estimate of sensi- 
tivities attainable using data from future tests. 

The second column of the table concerns free-fall WEP 
tests using falling corner cubes [70, 71]. In the second 
entry of this column, we present a single bound on the 
time-independent portion of the signal implied by exist- 
ing data. The remainder of this column lists crude es- 
timates of sensitivities that could be attained through 
sidereal and annual analysis of the same data. 
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The third and fourth columns of the table list sensitiv- 
ities from free-fall WEP tests using atom intcrferomctry. 
In the second entry of the third column, we present a sin- 
gle bound extracted from existing data [72]. The fourth 
column concerns proposals for future tests with atom in- 
terferometers [65], based on the idea that the relative 
vertical acceleration of two different atoms may be mea- 
sured using a simultaneous dual-species fountain [76]. 

The remaining columns of the table concern other pro- 
posed free-fall WEP tests. Crude estimates are provided 
of the sensitivities that might be achieved in the Principle 
of Equivalence Measurement (POEM) [73], via balloon 
drops in the General Relativity Accuracy Test (GReAT) 
[74] , and using the Bremen drop tower [75] . 

In the table, the estimates for the coefficients listed in 
the second and third rows and for (c")q in the penul- 
timate row all arise from the time-independent compo- 
nent of the data. A nonzero signal for any of these mea- 
surements would therefore be challenging to distinguish 
from other potential sources of WEP violation. Note 
that obtaining the independent sensitivities in the third 
and penultimate rows requires combining data taken in 



free-fall WEP tests performed at different colatitudes x- 

Independent sensitivities can also be achieved via other 
techniques. One possibility is positronium intcrferomctry 
[77], which via comparison with ordinary matter could 
yield a bound on different linear combinations of (o^)t 
and {c^)tt- Also, some independent measurements can 
be extracted by combining results from free-fall WEP 
tests with those from the force-comparison WEP tests 
discussed below. 

The next generation of the POEM experiment [73] is 
the proposed Sounding Rocket POEM (SR-POEM) [78], 
which is a WEP test designed to measure to 10^^^ the rel- 
ative acceleration of freely falling test bodies on a sound- 
ing rocket during certain phases of its flight. Although 
not terrestrial, this experiment can also be analyzed us- 
ing the methods presented here. A competitive sensitiv- 
ity is anticipated for measurements of the combination 
a{alP^^"')T — ™^(c'^"'"^-")tt/3 of isotropic coefficients. 
Obtaining sufficient data to resolve the periodic changes 
necessary for sensitivity to other coefficient combinations 
would be challenging. 



Table VI. Sensitivities for free-fall WEP tests. 





Falling 


Atom 


Tossed 


Balloon 


Drop 


CoefScient 


corner cube 


interferometry 


masses 


drop 
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combination 


[70, 71] 


[72] [65] 


[73] 


[74] 


[75] 


Aa 


10-10 


10-7 {10-17} 


{10-14} 


{10-15} 


{10-12} 


a(aeT"")T - ^mPic^+P-^TT 












+ (icos2x-i)m"(c")Q 


10-« GeV 


10-5 GeV {10-15 GeV} 


{10-12 GeV} 


{10-1-^ GeV} 


{10-10 GeV} 


a(a:+^-")T-imP(c^+^'-")TT 




{10-15 GeV} 


{10-12 GeV} 


{10-13 GeV} 


{10-19 GeV} 




[10^4 Q^v] 


{IQ-ii GeV} 


{10-8 GeV} 


{10-9 GeV} 


{10-6 GeV} 


( — e+p— n\ 


[10--* GeV] 


{IQ-ii GeV} 


{10-8 GeV} 


{10-9 GeV} 


{10-6 GeV} 




[10-2 GeV] 


{10-9 GeV} 


{10-6 GeV} 


{10-7 GeV} 


{10-4 GeV} 




[10-2 GeV] 


{10-9 GeV} 


{10-6 GeV} 


{10-7 GeV} 


{10-4 GeV} 


(c")q 




{10-15} 


{10-12} 


{10-13} 


{10-10} 




[10-4] 


{10-11} 


{io-«} 


{10-9} 


{10-6} 



E. Force-comparison WEP tests 

Typical force-comparison WEP tests can be viewed 
as comparing the motion of two or more bodies joined 
through electromagnetic forces with that predicted by 
an equation of the modified form (132). The predicted 
motion depends on the details of the configuration, so a 
unified analysis for all force-comparison WEP tests is im- 
practical. Here, we consider as an illustration a sensitive 
existing force-comparison WEP test based on a torsion 
pendulum [33, 79]. Exceptional sensitivity to Lorentz vi- 
olation can be achieved using a torsion pendulum with a 
spin- weighted bob [35, 41, 80], but here we treat instead 
a bob with a dipolar composition. We remark in pass- 
ing that another interesting option for force-comparison 
WEP tests is the use of superconducting gravimeters to 



compare gravitational forces on parts of the Earth having 
different compositions [81], although present sensitivities 
to coefficients for Lorentz violation are likely to be some- 
what weaker. 

For the torsion pendulum, a simple model of the bob 
is a dumbbell viewed as a rod with test bodies 1 and 2 
placed on each end. The two test bodies are composed 
of different materials, and the bob is suspended by a tor- 
sion fiber attached at the midpoint between them. The 
resulting pendulum is typically rotated in the laboratory 
to improve the modulation of the signal. The relevant ob- 
servable in such tests is the twist angle ^(T) of the torsion 
fiber. This angle can be calculated from the Newton sec- 
ond law as modified by the presence of Lorentz violation. 

At zcroth order in Lorentz violation, the pendulum 
hangs at an angle C, from the local vertical in the labo- 
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ratory given by Eq. (126). Lorcntz-violating corrections 
to C exist, but these make no contribution to the signal 
to the order at which we work. This angle represents 
the equilibrium position for the swing mode of the pen- 
dulum. Lorentz-violating modifications to this position 
could drive small excitations of the swing mode, but ex- 
periments tuned to the torsion mode are typically com- 
paratively insensitive to other modes [79]. 

The orientation of the bob about the axis perpendicu- 
lar to both the torsion fiber and to the dumbbell dipole 
moment can also be considered. This is the equilibrium 
position for the wobble mode of the pendulum. For sim- 



plicity, we assume here that the bob is suspended at a 
point P equidistant between the centers of mass of the 
test bodies and that the test bodies are constructed to 
ensure the dumbbell is perpendicular to the torsion fiber. 
In the absence of Lorentz violation, this implies equality 
of the two masses mi and TO2. However, in the presence 
of Lorentz violation, mi and m2 differ at leading order in 
the coefficients for Lorentz violation. As the pendulum 
rotates, this difference could shift the dumbbell orienta- 
tion away from its equilibrium point and generate small 
excitations of the wobble mode about P. Again, experi- 
ments are comparatively insensitive to this mode. 



Table VIL Amplitudes for torsion-pendulum tests. 



Amplitude 



Phase a„ 



Lui.-Q 



-\ {jn^{T^)xx + m^{c^)YY - 2a(a^a.)T) (l + 



sin X ) sin 2% 
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W^V"® [(c"')(TY) sin 77- (c"')(j,2)(l - cos 77)] sin^x 



Lco^+uj+ii = im™V© [(c"')(TY) sin 77 - (c™)(tz)(1 - cos 77)] sin^ x 

Ku,e+u;-n = -kni^V(s{c^)(^TX) sin?7sin2 x 

Luj.+w-a = ImTVis [{c^)(TY) sin 77 -|- (c^)(tz)(1 + cost?)] sin^ x 

K^^-^+n = -\m'^VQ){'ff")(TX) sin77sin^x 

-fwe-o^-n = -\rn^VQ{c'")^TX) sin77sin^x 

T _ _ ] 

-^We— u;— f2 — " 
Lue+2u+Q. 

Koj,+2w-n = -5m"'ye(c^)(TY)(l + COST?) sin x(l + cosx) 
Lu,,+2u,-u = -5r?2'"ye(c"')(TX)(l + cos77)sinx(l + cosx) 
Kcj^-2uj+n = -5m"'V©(c"')(TY)(l + cost;) sinx(cosx - 1) 
Luje-2w+Q = 5m"'y©(c™)(Tjf)(l + cos77)sinx(cosx- 1) 
Kue-2uj-n = |"^"'^e(c^)(TY)(l - cos 77) sin x(cos x - 1) 
Lu,,-2u-n = |772"'ye(c^)(TX)(l - COST?) sin x(l - cosx) 



i7n"'ye(c"')(TY) (1 - cos r?) sin x(l + cos x) 
i77T."'y©(c"')(Tx)(l - cost;) sinx(l + cosx) 



4>-tp 

4> + 2%l) 
(p + 2ijj 

(p + ip 
^ + V 
(/) + V 

(t> + tp 

^- V 

(ji-Tp 
(/> + 27/> 

+ 2^ 
</) + 2^ 
(l) + 2iP 
(t)-2xP 
(f)-2'4) 
(l)-2tp 
(f>-2il) 



To analyze the torsion mode, it is convenient to express 
the relevant contributions to the difference rrii — t«2 of 
the test-body masses in a form displaying the dependence 
on particle species. This form can be obtained from Eq. 



(133), giving 

mi - 7712 = - ^ (iVf - N^) (2a«ff)r + m^ic^M • 

W 

(149) 
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At leading order in Lorentz violation, the oscillations of 
the system are determined by the second-order differen- 
tial equation 



27/ 



_de 



k9 



(150) 



where 7 is the torsional damping constant and k is the 
torsional spring constant. The moment of inertia / can 
be taken as 



I = (mi + m2)rl + 1', 



(151) 



where ro is the distance from P to the test bodies and /' is 

the moment of inertia of the remaining matter comprising 
the dumbbell. The torque r includes Lorentz-violating 
effects and is determined by the forces on the test bodies 
calculated from the PN0(3) lagrangian La.Ls- 

The damping term in Eq. (150) ensures that free oscil- 
lations vanish in the steady state. The time dependence 
of the steady-state solution is therefore determined by the 
rotations of the pendulum relative to the Sun-centered 
frame. Neglecting possible torques other than those im- 

(3) 

plied by La.Ls, the steady-state solution can be written 
in frequency-decomposed form as 



6{T) 



n.w 



X [Kn sin(a;„T + ^„ + a„) 
+L„ cos(a;„T + /3„ + a„)] , (152) 

where (3n = '2juJn/ (wg — w^) and where a„ is a phase 
fixing the relationship between the time coordinate in 
the turntable frame and the Sun-centered time T. The 
amplitudes K^, L„ and the phase a„ are given in Table 
VII. With the exception of the first row in the table, these 
signals for Lorentz violation are distinguished from other 
potential sources of WEP violation by their characteristic 
time dependence. 

In the above analysis, the assumption of a steady- 
state solution implies the pendulum motion is governed 
by leading-order Lorentz violation, while the torque r 
is taken as the only relevant source of Lorentz viola- 
tion. Note that Lorentz-violating contributions to the 
moment of inertia / can be neglected here because they 
enter only at higher order. These are a manifestation 
of the angular-momentum nonconservation that accom- 
panies Lorentz violation, and they are analogous to the 
Lorentz-violating contributions to the effective inertial 
mass in the modified Newton second law (132). In prin- 
ciple, the rotation of the apparatus in the laboratory in- 
troduces similar effects proportional to cjk and uilrQ. 
These may be comparable in magnitude to effects listed 
in Table VII that are suppressed by co^R^V^, but they 
offer no additional advantage in terms of sensitivity and 
so are disregarded here. 

The above analysis can be used to extract constraints 
on Lorentz violation from the results of the torsion- 
pendulum WEP tests reported in Refs. [33, 79]. The 



attained sensitivity to the differential acceleration of Be 
and Ti test bodies at the level of 10^^^ ms^^ [33] is the 
experimental basis for our limit (121), which extends an 
earlier bound [11] to include the coefficients Cj^. In the 
remainder of this subsection, we revisit this issue to in- 
corporate the slightly weaker constraints from torsion- 
pendulum experiments using Al, Be, Cu, and Si test bod- 
ies [79], and we consider implications of nonzero (c'")^!/. 

First, suppose (5"")^^ = 0. Inclusion of data from tests 
with different materials permits the extraction of some 
independent sensitivities to neutron coefficients and to 
combinations of electron and proton coefficients. This 
treatment relies on differences in binding energy be- 
tween the materials involved, so the signal sensitivity of 
10~^^ ms~^ relative to w^ii© ~ 3 x 10~^ ms~^ is sup- 
pressed both by the typical material-dependence factor 
of (iVf - N^) ~ 10-2 appearing in Eq. (152) and by 
another order of magnitude from the binding-energy dif- 
ference. By combining available Be-Ti and Al-Be data 
[79] we obtain the estimated bounds 



l«(«eT)T| < 10-'° GeV, 
|a(«eff)T| < 10"'° GeV, 



(153) 



vahd for (c"')^^ = 0. 

If instead nonzero coefficients {c^)ij,v are present, then 
we obtain the estimated bound 



|a(a:+^-")T-im^(c'=+^-")TT- 



'g"^"(c")Q| < 10-" GeV. 

(154) 

The contributions due to the spatial neutron coefficient 
(c")q cannot be disentangled from those due to the tem- 
poral components at this order in the analysis. However, 
this separation becomes feasible when the result (154) 
is combined with the limit achieved via free-fall WEP 
tests given in row 2 of Table VI. We thereby obtain the 
constraints 



l"(«:ff''"")T-3m^(c'^+^"")TT| < 10-«GeV 

l/_n,\ I ^ 1,^-8 



< 10" 



(155) 



As discussed following Eq. (125), the possibility of 
(c"')^!/-type Lorentz-violating effects in the binding en- 
ergy impedes its direct use in extracting independent 
sensitivities to \a{alg^)T—'^m^{c^~^P)TT\ and \a{a^Q)T — 
imf(c")TT|. 

In addition to the constraints (153), (154), and (155), 
other new bounds could be placed on the moduli of cer- 
tain coefficients for Lorentz violation by reanalyzing the 
time dependence of the data obtained in the experiments 
of Refs. [33, 79] using the result (152). Crude estimates 
of these sensitivities are given in Table VIII. These are 
obtained disregarding binding-energy considerations but 
making the strong assumption that all relevant frequen- 
cies in Table VII can be studied in the data. Allowing for 
binding- energy effects could yield independent sensitivi- 
ties to the neutron coefficients and to a combination of 
proton and electron coefficients, both reduced by roughly 
a factor of 10. 



34 



Table VIII. Sensitivities for torsion-pendulum tests. 



Coefficient 



Sensitivity 



'-{TJ) 
/—e-\-v—n\ 

oi{%s )y+z 

/—e+p—n\ 
/ — e+p— n \ 



[10-7] 
[10-8 GeV] 
[10-7 GeV] 

[10-8 GcV] 
[10-7 GeV] 



VIII. SATELLITE-BASED WEP TESTS 

Space-based platforms offer certain advantages in tests 
of gravity [82] and searches for Lorentz violation [83]. 
The long free-fall times that may be attainable on 
a drag-free spacecraft make satellite-based WEP tests 
particularly attractive. Several proposals are in an 
advanced stage of development, including the Micro- 
Satellite a trainee Compensee pour 1' Observation du 
Principe d'Equivalcncc (MicroSCOPE) [84], the Satellite 
Test of the Equivalence Principle (STEP) [85], and the 
Galileo Galilei (GG) mission [86] . A WEP reach similar 
to that of STEP has also been suggested for the Grand 
Unification and Gravity Explorer (GaUGE) mission [87] . 

The basic idea underlying these missions is to mon- 
itor the relative motion of test bodies made of differ- 
ent materials as they orbit the Earth in a satellite. In 
the presence of nonzero coefficients for Lorentz violation 
(o^)jii and {c^)iMv, the orbits of the test bodies become 
material dependent. In this section, we determine the re- 
sulting apparent WEP violations and then obtain crude 
estimates of the sensitivities to (aeir)A' ^^'^ attain- 
able in MicroSCOPE, STEP, and GG. 



A. Theory 

The basic observable for a satellite-based WEP test is 
the differential local acceleration between the test bod- 
ies. The typical design goal is to achieve excellent sensi- 
tivity to one or two components of this acceleration. For 
present purposes, we can idealize the situation as a pair 
of test bodies aboard a satellite traveling in a circular or- 
bit. In what follows, we allow for the possibility that the 
test bodies are also rotating about an axis perpendicular 
both to the direction of motion of the satellite and to the 
direction of acceleration sensitivity. 

Some notation relevant for our analysis of satellite- 
based WEP tests is summarized in Table IX. Paralleling 
the analysis of terrestrial experiments in Sec. VII A, it 
is convenient to introduce an Earth-centered frame with 
coordinates ~ {t,x,y,z), chosen so that t = T and 
so that the spatial components match those of the Sun- 
centered frame at leading post-newtonian order. The 
Earth-centered coordinates can be related to the Sun- 
centered ones as discussed in Sec. VIA. The angles ^i. 



^2 in the table arc defined relative to the basis vectors of 
the Earth-centered frame. The notation for properties of 
the test masses 1 and 2 follows that of Sec. VI B. 

Table IX. Notation for satellite-based WEP tests. 

Quant ity Dcliuit ion 



^1 

62 



mean Earth radius 
mean Earth orbital speed 
Earth-satellite separation 
satellite orbital frequency 
satellite rotational frequency 
inclination of satellite orbit 
longitude of satellite-orbit node 
phase fixing satellite location at T = 
phase fixing satellite orientation at T = 



Establishing the signal arising from nonzero coeffi- 
cients for Lorentz violation requires the transformation 
from the Sun-centered frame to a frame comoving with 
the satellite. The satellite frame serves as the equiva- 
lent of the laboratory frame for terrestrial searches. We 
denote coordinates in the satellite frame by x'^. 

Since the satellite orbit is inclined relative to the Earth- 
centered frame, it is also useful to introduce an interme- 
diate frame aligned with the satellite orbit and hence 
rotated with respect to the Earth-centcircd frame. The 
intermediate coordinates are denoted by x** . The rota- 
tion transformation from to can be written as the 
matrix 



/ cos ^2 — COS ^1 sin ^2 sin ^1 sin ^2 
sin ^2 COS ^ 1 cos ^2 — sin ^1 cos ^2 

c COS^I 



v 







sinCi 



\ 



(156) 



using the angles ^1 and ^2 defined in Table IX. 

The connection between the satellite coordinates and 
the Earth- centered coordinates can be written 



x^ = Ri'' {R^^'x- 



k'i J 



„k'\ 



(157) 



Here, x^ is the world line of the satellite in the in- 
termediate coordinate system. This world line can be 
parametrized as 



[r cos{wsT + ei),r sm{wsT + 0i),O), (158) 



where r is the magnitude of the Earth-satellite separa- 
tion. The satellite therefore orbits in the x'-y' plane. 
Also, in Eq. (157) the rotation i?2 ' of the satellite is 
given by the matrix 



R: 



k'i 




— sm{uJrT - 

COs{uJrT -i 





62) 

92) 



(159) 



The axis of the satellite rotation is therefore along z. 
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Table X. Amplitudes for satellite-based WEP tests. 



Amplitude 



phase 



= m'^ruJs 

5r 



vR^a cos ^2 



(a-a.)xcosa + raisin ^: 



Qu>,-CJr = 



+ [«s)x COS a 4 

) cos^i COS ^2 - (c"')(TZ) sin ^2 



sin^ 



Q2U1, 



{c^)(TY) sinCi + {c'")(TX) cosCi 
(c"')(TJf) sin^i cos ^2 - (c'")(TF) 
(flciiOx sin^i - {a^s)Y cos^i 
cos ^2 cos 2^i(c"^')(xy) + sin ^2 sin ^1 (c"')(y 2) 
+i sin 2^1 cose2((c"')yy - (c'")xx) + sin^ cosei(c'")(xz) 
(cos2 a cos2 a - sin^ a + i sin2 6)((c"')xx - (c'")rY) + ^ 
- cosa sin2a(c"')(yz) +sina sin2a(c'")(xz) +sin2a(l 
(cos2 a sin2 a + i cos2 a + \) {{^)xx - {^)yy) 
-i sin2 a {{lf')xx + {c}")yy - 2{c^)zz) + 2{1^)yy 

+ sin2a (l-cos^a) (c"')(;fy) -sinasin2a(c'")(x2) +cosasin2a(c'")(y2)] - MKs)t 

(c )(TX) cosa + (c"^)(Ty) sma 

a cosa - (c"')(TX) sina cosa + (c"')(tz) sina 



a((c"')xx + (c"')YF-2(c"')zz) 



+ cos2 a)(c"')(xy) 



(a^g.)x cosa + «ff)y sina 



5r 



(c"')(Ty) cos 

- (a^gOxsina - (a^fiOi-cosaj 
(cos^ a ^ sin^ a cos^ a ^ cost; cos a cos 2a — sin 77 sin a cosa) (c"')(tx) 
+ sina sina (cosa - cosry) {c^){tz) 



(cosa + cosa cos^a ^ sin 77 sina ^ 2cos?7Cosa cosa) sina(c"')(Ty) 
a cos a - cos cos a (1 + cos^ a)) sina(c"')(TX) 



(2 COS a cos a ^ sin 77 sin ; 

— (cos 2a cosa ^ sin 77 cos a sina cosa + cos 77(1 — cos^ a sin^ a)) (c"')(Ty) 



] sin a - cos rj cos a ) sin a (c"" ) (tz) 



— (cos a — sin rj 

Pn+^,->^^ = m^V(s [ (1 - sin2 a sin^a) (c'")(tx) + | sin 2a sin^ a(c'")(TY) - 1 sina sin 2a(c"')(TZ) 

„ T A /— 7/J N 



i (cos 77 sin 2a sin^ a — sin 77 sin a sin 2a) {c^)(tx) 

1 - sin^ a cos^ a) cos 77) {c^){ty) 



+ (i sin 77 cos a sin 2a + siy--"'/yv^ Ai yj 

+ (sinr7sin^a + | cos 77 cos a sin 2a) ^)(tz) + aV© [(a^)z sinT? + (a^fj)r cost? 
(1 - sin^ a sin^ a)(c^)(TX) + 5 sin 2a sin' a(c'")(Tr) - 5 sina sin 2a (c'")(tz)^ 
-aye(«^ff)x 



<3n-u-,+u-r = rri^Vfs \ (sin 77 sina sin2a - cosr7sin2a sin' a) (c^)(tx) 

- (I sin 77 cosa sin 2a + cos 77(1 - cos' a sin' a)) ^'"){ty) 

- (sinr7sin'a + 5 cos 77 cos a sin 2a) (c"')(tz) + aV'e («e 
(cos' a ~ sin' a cos' a + sin 77 cos a sin a + cos 77 cos 2a c 



(a^jf)z sin 77 + (a^fj)r cos 77 
2acosa)(c"')(TX) 



Qf2— CUs— Wr- 



(cos' a ~ sin' a cos' a + sin 77 cos a sin a + cos 77 cos 2a cos a) {(^){Ty. 
+ ( i sin 2a ( 1 + cos' a ) + sin 77 sin a sin a + cos 77 sin 2a cos a) (c"' ) {ty) 
+ (i sin 2a + cos 77 sin a) sina(c"')(TZ) 



(i sin 2a + cos 77 sin a) sina(c"')(TZ) 
- (sin 2a cos a + ^ sin 77 sin a sin 2a + ^ cos 77 sina (1 + cos' a)) (c'")(tx) 
+ (cos 2a cos a + ^ sin 77 cos a sin a — cos 77 (sin' a — cos' a cos' a ) ) (c^ ) (i 
+ (cosa sina + sin 77 sin' a + 5 cos 77 cosa sin 2a) (c"')(tz) 



f 1 — W2 

W\ — O2 

201 — 62 



n -I- C2 



71 - C2 



71 - t72 



-Qx + 62 



-Ol + 02 



-01 — 02 



71 - (72 
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For our purposes, it sufSccs to obtain explicitly the 
local differential acceleration Aa"^ of the test bodies in 
the X direction. We have 

Aa- = ^ = AaSdai + Aa£v + .... (160) 

The first term on the right-hand side of this expression 
is the conventional Newton tidal term. It takes the form 

+wl + \ujf)/^x. (161) 

The second term in Eq. (160) contains Lorentz-violating 
contributions to the differential acceleration. It can be 
written 

'^>£v = ~?i:(f -^) M 

X {Pn sin(a;„T + a„) + Qn cos(a;„T + Q!„)^ . 

The amplitudes Pm, Qm and the corresponding phases 
are provided in Table X. Finally, the ellipsis in Eq. 
(160) represents higher-order gcneral-relativistic correc- 
tions and Lorentz-violating effects at the same post- 
newtonian order as Aa^^g^j. The latter are typically of 
lesser interest. If desired, the differential acceleration 
Aa^ along y can be obtained by performing the transfor- 
mation oJrT oJrT - Tr/2 on Eq. (160). 

B. MicroSCOPE and STEP 

Within our idealized scenario, MicroSCOPE [84] and 
STEP [85] can be analyzed in parallel. Each apparatus 
consists of a pair of cylindrical test bodies made of differ- 
ent material but having a common symmetry axis. The 
test bodies are free to move along this axis. In satellite 
coordinates, this direction lies along x and is perpendic- 
ular both to the direction of motion of the satellite and 
to the axis of the satellite rotation. 

One prosaic origin of relative motion of the test bodies 
along the x direction could be the influence of tidal forces 
on a misalignment of the two centers of mass, which 
would lead to the acceleration Aa^j^^j in Eq. (160). This 
can be separated from the acceleration due to WEP vi- 
olations stemming from Lorentz-invariant sources, which 
enters with the characteristic frequency lUs — LUr- Here, 
we are interested in a WEP- violating acceleration Aaf^y 
arising from the coefficients (aejf);^ and (c*^ )^,/ for Lorentz 
violation. This can be distinguished from both the above 
effects through careful separation of the frequencies as- 
sociated with the amplitudes in Table X, except for the 
amplitude Qus-^j,-- 

The sensitivity goals of MicroSCOPE and STEP are 
Aa/ro;^ < lO^^''' and Aa/rw^ < 10~^*, respectively. 
These sensitivities and the results in Table X can be 



used to obtain rough estimates of the reach of these ex- 
periments for studies of Lorentz violation. For this pur- 
pose, we take the quantity /mi — N2 /m^ appearing 
in Eq. (162) to be of order 10~^ GeV^^, which is the 
best available value with the Pt-Ir, Be, and Nb test bod- 
ies presently proposed for STEP. Note that the bounds 
scale linearly with this difference, so a careful choice of 
test-body material can maximize sensitivity to Lorentz 
violation. Moreover, combining results for different test 
materials can yield additional independent sensitivities. 
Note also that the experimental reach may vary with the 
choice of orbit. For dcfiniteness, we suppose the sines 
and cosines of ^1 and ^2 are of order one. 

Our crude estimates for attainable sensitivities to the 
moduH of (a^fj)^ and (c"')^^ for MicroSCOPE and STEP 
are presented in Table XL In each row, the listed sensi- 
tivities are obtained under the assumption that all coef- 
ficients vanish except those appearing in the first entry. 
The key factor underlying the difference in reach for the 
various coefficient combinations is the boost entering the 
relevant amplitude in Table X. Amplitudes containing Vq 
are suppressed by roughly 10~^, while those containing 
rujs are suppressed by about 10^^. As before, the braces 
indicate the estimated sensitivities involve data from fu- 
ture tests. 



C. Galileo Galilei 

Certain design features of GG [86] differ from those 
of MicroSCOPE and STEP in ways that are significant 
for studies of Lorentz violation. Although GG also uses 
coaxial cylindrical test bodies, it is sensitive to accelera- 
tions in the plane perpendicular to the axis of the cylin- 
ders. Also, the cylinders are rotated about their axis at 
a comparatively high frequency of about 2 Hz. 

In applying the generic analysis of Sec. VIII A to GG, 
it is convenient to take the cylinder axes to lie along 
z. The experiment is then sensitive to accelerations in 
the x-y plane. The differential acceleration Aa^ along 
X is given in Eq. (160), while Aa^ can be obtained by 
adjusting the phase 62- 

The sensitivity goal of GG is Aa/rwf < 10"^'^. In 
Table XI, we present rough estimates of the correspond- 
ing reach for measurements of the coefficients (aj|^ )^ and 
{(?")tiv for Lorentz violation, obtained using the result 
(162). The values for GG in the table are based on the 
same assumptions as those discussed above for Micro- 
SCOPE and STEP. This includes the material-dependent 
factor, with the proposed materials for the GG test bod- 
ies being Be and Cu. The boost factors leading to the 
varying sensitivities for GG listed in the table are also 
of the same order of magnitude as for the other satellite 
experiments. 

We remark in passing that the comparatively high ro- 
tation rate for the GG cylinders could introduce addi- 
tional Lorentz-violating effects. Typically, the presence 
of nonzero (c™)^!/ introduces modifications to the effec- 
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tivc moment of inertia of a body. This can affect the dy- 
namical balance of the system, which can lead to observ- 
able signals. For example, potential effects of this type 
on the timing of pulsar signals have been used to con- 
strain some combinations of (c")^^ [88]. In the present 
context, the observable signals could include a material- 
dependent Lorentz- violating wobble varying at the satel- 
lite frequency and at the Earth's orbital frequency. It is 
conceivable that these Lorentz-violating effects could be 



detected by the GG apparatus that senses the test-body 
location. Notice that the signals would be independent 
of gravity. They may be detectable using sophisticated 
terrestrial dynamical-balancing equipment, perhaps in- 
cluding that used in the Galileo Galilei on the Ground 
(GGG) experiment [89]. The investigation of these ef- 
fects represents an interesting open question for future 
work. 



Table XI. S<-usi( ivii i<-s for sat.<-llit.<--bas(-(l WEP lorils. 



Coefficient MicroSCOPE GG STEP 

a(a^^P"")T- imP(c'^+f-")TT {lO-^''* GcV} {IQ-^^ GeV} {lO-i^ GeV} 

aials^'lx {10-9 GeV} {10"" GeV} {IQ-i^ GeV} 

a(o^+^'-")y+z {10-9 GeV} {10"" GeV} {IQ-^^ GeV} 

a{al+P-")Y {10-7 GeV} {10-^ GeV} {IQ-i" GeV} 

a{al+P~")z {10-7 GeV} {10-^ GeV} {IQ-io GeV} 

(c")q {10-13} {10-15} {10-16} 

(C")(TJ) {10-9} {10-11} {10-1^ 



IX. EXOTIC GRAVITATIONAL TESTS 

In this section, we offer a few remarks about some 
gravitational searches for Lorentz violation using mate- 
rial test bodies other than neutral bulk matter, neutral 
atoms, or neutrons. These more exotic searches typically 
present unique experimental challenges, but they could 
provide access to combinations of coefficients for Lorentz 
violation that are awkward or impossible to isolate and 
measure in other searches discussed in this paper. Here, 
we briefly consider tests with electrons and ions, studies 
with antihydrogen, and experiments using particles from 
the second and third generation of the SM. 



A. Tests with electrons and ions 

Measurements of the gravitational acceleration of 
charged matter remain of definite theoretical interest be- 
cause the WEP and other foundational aspects of gravity 
are comparatively poorly tested in this regime. In this 
subsection, we consider possible signals from studies of 
charged electrons or ions. Given the experimental chal- 
lenges of these tests and their limited attainable sensi- 
tivities, we restrict attention here to effects from (a^g)^, 
setting other coefficients to zero for simplicity. 

In the context of searches for Lorentz violation, grav- 
itational tests with charged matter offer unique access 
to the coefficients (a^g)^. For example, measurements of 
this kind can disentangle coefficients for Lorentz violation 
in the proton and electron sectors. They can also detect 
certain countcrshadcd effects that are otherwise invisi- 
ble. In particular, some models have coefficients (a^)/i 
proportional to electric charge, which would evade de- 
tection in searches with neutral test bodies [11]. This 
possibility is a natural consequence for theories in which 



the photon modes are interpreted as NG bosons from 
spontaneous Lorentz breaking and in which (a^)^ re- 
mains physically observable, such as nonminimally cou- 
pled bumblebee electrodynamics [14]. 

One candidate technique to measure gravitational ef- 
fects from the coefficients (a^fj);^ is charged-particle in- 
terferometry. Electron interferometry has been used to 
measure the Sagnac effect at the 30% level [90] , while ion 
interferometry is under investigation as a practical tool 
for sensitive tests of Coulomb's law [91]. In the present 
context, electron or ion interferometry offers an interest- 
ing alternative prospect to the free-fall tests with neu- 
tral matter discussed in Sec. VII. For a given geometry, 
the observed phase shift can be determined using the; 
methods of Sec. VII B. In the limit of interest here, the 
vertical acceleration a^^ of the electron or ion T in the 
gravitational field of the Earth S is given at PN0(2) by 

As before, the PN0(3) version of this acceleration can 
be frequency decomposed relative to the Sun-centered 
frame, with the corresponding amplitudes depending on 
the coefficients {a^g)ij, as given in Table IV. 

In principle, a charged-particle interferometer can be 
used for free-fall gravimeter tests of the type discussed 
in Sec. VII B or for free-fah WEP tests as in Sec. VII D. 
A free-fall gravimeter test is insensitive to {a^g)T and 
has only boost-suppressed signals from (a^)j, so a sub- 
stantial improvement over the existing reach of charged- 
matter interferometers would be required to achieve a 
sensitivity compatible with perturbative consistency. In 
contrast, a free-fall WEP test is directly sensitive to 
(^eff)T t)ut requires a simultaneous measurement with 
two test bodies. One option along these lines could be a 
direct comparison with neutral matter via a falling corner 
cube or an atom interferometer. 
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Another approach to gravitational tests with charged 
matter is to study the motion of charged particles in a 
vertical metallic drift tube. This setup is accompanied by 
gravitationally induced electric forces caused by the sag- 
ging of the tube [92] , along with a variety of challenging 
systematics. An experiment of this type with electrons 
[93] confirmed that the gravitational forces on the elec- 
trons in the tube and on the electrons within the metal 
are comparable to about 10%. An analogous experiment 
involving cold antiprotons [94] was designed to achieve a 
sensitivity of 0.1% to the gravitational acceleration [95]. 
These measurements are all experimentally challenging, 
and their interpretation is theoretically subtle [96]. 

In the present context of gravitational Lorentz viola- 
tion involving the coefficients (a^)^, intuition for the 
theoretical implications of ballistic tests of this type can 
be gained by considering the idealized case and working 
at PN0(2). We suppose a test particle T of charge 
moves along the symmetry axis of a vertical cylindrical 
metallic drift tube with body comprised of a lattice of 
ions of type / and conduction electrons e. Disregarding 
applied fields, stray fields, and various systematics, the 
overall conventional force on the particle T is the sum 
of the direct gravitational force on T from the Earth S 
and the net force on T from the electromagnetic field 
arising from the gravitationally induced sagging of the 
tube. The presence of nonzero coefficients (a^)^ intro- 
duces corrections to both these forces. At PN0(2), the 
gravitational force on T is given by 

(-Fgrav)^ = -rn^g-2ga{ajfi){-2ga-^{a^fi)v 

(164) 

while the vertical component of the force on T from the 
gravitationally induced electric field is 

+7^ (m'g + 2ga{ais)i + 2ffa^(afff),-^ . 

(165) 

In these expressions, m*^ and are the masses of an 
electron and an ion in the tube lattice, respectively, while 
e and are the corresponding charges. The factor 7 is 
a constant, set by the properties of the metal lattice. In 
Eq. (165), the first three terms arise from the sagging 
of the electrons in the tube walls, while the last three 
are proportional to the dilation derivative of the work 
function for the metal and originate in the longitudinal 
compression of the lattice. These expressions reduce to 
standard ones when the coefficients {a^g)fj. vanish. 

Although the expressions (164) and (165) hold in an 
idealized situation, they suffice to demonstrate in princi- 
ple that experiments of this type are sensitive to nonzero 
coefficients (aJJ^)^, even when these coefficients are un- 
detectable with neutral matter. This is also true if the 



particle T is an electron, when the sum of the forces 
(^grav)z ^iiid (-Fcm)i Icavcs Only the last three terms in 
Eq. (165). In practice, however, the reported reach of 
drift-tube experiments to date is insufficient to achieve 
useful sensitivity in gravimeter tests. A WEP test relat- 
ing a drift-tube setup to an independent gravimeter may 
be of more interest. For a given experiment, specific sen- 
sitivities can be estimated using the analyses presented 
in Sec. VII. 

A third methodology for investigating gravitational 
Lorentz violation from nonzero coefficients (a™^)^ could 
conceivably be to adopt as the gravimeter a device that 
wholly confines charged particles. For example, a single 
charged particle can be trapped for long periods using 
a Penning trap [97]. Measuring gravitational effects in 
this way is ambitious, as can be appreciated from the 
size of the quantity \m'^g/e\ ^ 6 x lO"""^^ V/m. Nonethe- 
less, the feasibility of gravitational measurements with 
trapped antiprotons at a sensitivity of about 1% has been 
suggested, using a gravity-induced shift of radial orbits 
[98] . This would also lead to sensitivity to the coefficients 
{a^ff)fi via an analysis similar to those discussed above. 



B. Tests with antimatter 

The study of antimatter offers another realm in which 
to search for Lorentz and CPT violation. Antihydro- 
gen has been detected [99, 100] and produced in copi- 
ous amounts [101], while prospects for studies of trapped 
cold antihydrogen are excellent [102] . Antihydrogen spec- 
troscopy could yield special sensitivity to nongravita- 
tional SME coefficients for Lorentz and CPT violation 
[103], and the experiment for Atomic Spectroscopy And 
Collisions Using Slow Antiprotons (ASACUSA) expects 
to achieve sensitivities of parts in 10^^ to the predicted 
shifts in hyperfine transitions [104]. 

To study the interaction of gravity and antimatter, var- 
ious ideas for measuring the gravitational acceleration 
of antihydrogen have been advanced. Among them are 
methods involving trapped antihydrogen [105], antihy- 
drogen interferometry [106], antihydrogen free fall from 
an antiion trap [107], and tests in space [108]. One ap- 
proved project, the Antimatter Experiment: Gravity, In- 
terferometry, Spectroscopy (AEGIS) [109], has an intero- 
ferometric design with an initial sensitivity goal of 1% to 
the gravitational acceleration of antihydrogen. 

In the context of gravitational Lorentz and CPT vi- 
olation, these experiments offer the prospect of special 
sensitivities to the coefficients (a^)^ and (c"*)^!/. The 
key point is that a CPT transformation has the net ef- 
fect of reversing the sign of (a^)/* while leaving (0""')^^ 
unchanged. As a result, experiments with antihydrogen 
could in principle observe distinctive and novel behav- 
iors. Moreover, when compared with similar measure- 
ments on hydrogen, the results would offer the opportu- 
nity for clean separation of effects. For instance, free-fall 
WEP tests comparing hydrogen and antihydrogen could 
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yield independent sensitivity to {c'^~^^)tt- In general, the 
theoretical treatment of prospective free-fall gravimeter 
or WEP tests with antihydrogen follows the same path 
as described in Sec. VII, except with the sign of {a^g)^ 
reversed throughout. 

The literature contains numerous attempts to place 
indirect limits on the possibility of unconventional 
antimatter-gravity interactions, many of which are re- 
viewed and critiqued in Ref. [110]. In the present con- 
text, the SME offers a general field-theoretic approach 
that can elucidate aspects of this issue and provide new 
insights about possible limitations on effects. We next 
present an explicit toy model that evades some previous 
indirect limits on large unconventional effects in antihy- 
drogen. 

For simplicity, we choose to work within the isotropic 
limit of the SME. In any specified inertial frame O, a sub- 
set of Lorentz- violating operators in the SME Lagrange 
density preserves rotational symmetry. Setting the coef- 
ficients of all other operators to zero produces an inter- 
esting limiting case. The frame O then becomes a pre- 
ferred frame, since the rotation invariance is broken in 
any frame O' boosted with respect to O. Physical effects 
of Lorentz violation are then isotropic in O but not in O' . 
The frame O could in principle be identified as the rest 
frame U of the cosmic microwave background (CMB), 
the Sun-centered frame S, or any other desired choice. 
Isotropic models of this type are sometimes called 'fried- 
chicken' models because of their popularity and simplic- 
ity. 

In Minkowski spacetime, toy isotropic models can be 
used to show that Lorentz- and CPT-violating effects 
could in principle be substantially larger in antihydro- 
gen than hydrogen. One example is the isotropic 'in- 
visible' model (IIM) [111], which is defined in the CMB 
frame U and yields effects challenging to see in searches 
with ordinary matter. Denoting coordinates in U by 
{T' , X' ,Y' , Z'), the IIM assumes the only nonzero co- 
efficients for Lorentz violation are (6p)t' and isotropic 
{d^)s'S' obeying the simple condition 

(6P)t' = kmP{(F)T'T' (166) 

for a suitable choice of constant k. In the Sun-centered 
frame S, this one-parameter model generates nonzero 
coefficients and {<P)jt- The dominant signals in 

terrestrial experiments with hydrogen appear in the hy- 
perfine structure and involve the combination — 
mP{dP)jT, which vanishes for suitable k. These experi- 
ments can therefore detect only effects suppressed by at 
least one power of the boost of the Earth around the Sun, 
which is about 10~^ and requires an experiment sensitive 
to annual modulations. In contrast, the dominant effects 
in experiments with antihydrogen involve the combina- 
tion (b^) J + mP (cP) jT , which produces unsuppressed sig- 
nals in the hyperfine structure. The IIM thus provides a 
toy field-theoretic scenario in which observable effects in 
antihydrogen arc at least 10,000 times greater than those 
in hydrogen or other nonrelativistic neutral matter. 



The IIM involves spin-dependent operators for Lorentz 
and CPT violation in Minkowski spacetime. In this 
work, the focus is on the gravitational couplings of 
spin-independent operators with coefficients (a^)/^ and 
{c^)fii^- At the end of Sec. VIII A, we remark on the dif- 
ficulty of observing with matter any signals depending 
on the combination a{aJg)T — rn^{'5^)TT/^ of isotropic 
coefficients. Here, we consider some implications for an- 
timatter gravity of a specific toy model, the isotropic 
'parachute' model (IPM) , in which unobserved combina- 
tions of this type provide the dominant source of Lorentz- 
violating effects and could yield significant a priori differ- 
ences in the gravitational accelerations of hydrogen and 
antihydrogen. 

To construct the IPM, consider the Lagrange density of 
the SME in the Sun-centered frame S, with nonzero co- 
efficients restricted to (aJfff)T and isotropic {c^jss- Fol- 
lowing the derivation in the early sections of this work, 
we can extract the PN0(3) effective classical lagrangian 
for a test particle T moving in the gravitational field of 
a source S. This can be written in the suggestive form 

LiPM = \mjv^ + (167) 

where is the effective inertial mass of T, while mJ 
and rrig are the effective gravitational masses of T and 
S, respectively. All these effective masses are defined in 
terms of the coefficients (a^)T, (c'")tt for Lorentz vio- 
lation and the body masses of Eq. (71). We find 

W 

= m^ + ^((Af'" + 7V'^)m'"(c'")TT 

W 

+ 2a{N^ -N^)ia^s)T), (168) 

where B is either T or S. These expressions adopt the 
notation N'^ and N"^ for the number of particles and an- 
tiparticles of type w, respectively, while as before m"' is 
the mass of a particle of type w. Note that for a given 
body the passive and active gravitational masses are 
identical, reflecting the preservation of Newton's third 
law in the model. 

For electrons, protons, and neutrons, the IPM is de- 
fined by the three conditions 

a«ff)T = k^^^ic^TT, (169) 

where w ranges over e, p, n. Since there are three inde- 
pendent conditions on six real parameters, this produces 
a three-parameter IPM. The condition (169) ensures that 
for a matter body B the effective inertial and gravita- 
tional masses are equal, 

mf = rUg (matter), (170) 

and hence no Lorentz-violating effects appear in gravita- 
tional tests to PN0(3) using ordinary matter. However, 
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for an antimatter test body T this condition fails, 

mj ^ mj (antimatter), (171) 

so observable signals arise in comparisons between the 
gravitational responses of matter and antimatter or be- 
tween different types of antimatter. Ensuring the valid- 
ity of perturbation theory requires that the coefficients 
a{a^fi)T = {c^)tt/^ are perturbatively small relative 
to . With theoretically conceivable values perhaps 
even as large as CSm"", the gravitational accelerations of 
hydrogen and antihydrogen might differ at the 50% level. 

Rather than a serious effort at a realistic theory, the 
IPM is constructed as a simplistic playground within 
which to explore field-theoretic limitations on unconven- 
tional properties of antimatter and antihydrogen. In the 
next few paragraphs we treat it as such, briefly address- 
ing some concerns about unconventional signals in this 
context. 

One issue is whether energy remains conserved when 
matter and antimatter have different gravitational in- 
teractions [112]. For the analysis of the SME in the 
present work, this issue is moot because an explicit 
conserved energy-momentum tensor exists. As an il- 
lustration, consider the gcdanken experiment in which 
a particle-antiparticle pair is lowered in a gravitational 
field, converted to a photon pair, raised to the original lo- 
cation, and finally reconverted to the particle-antiparticle 
pair. In generic scenarios the particle, antiparticle, and 
photons each provide different contributions to the en- 
ergy and so problems can arise. However, in the IPM 
these complications are avoided. The photons make no 
contribution because they arc conventional, partly via 
the coordinate choice (23). The particle and antiparticle 
do contribute to the energy via the coefficient (a^)T, but 
the two contributions cancel. Contributions involving the 
coefficient (c"')tt exist and combine during the lowering 
procedure, but the definition (70) of the conserved energy 
also contains (c^)rT and so the net change remains zero 
at the end of the experiment. The resolution of this and 
other illustrative scenarios is less transparent when more 
nonzero coefficients for Lorentz violation are present, but 
the existence of a conserved energy-momentum tensor en- 
sures that no contradictions arise. 

Another attempt to argue against the possibility of an 
anomalous antimatter response to gravity is based on the 
large binding-energy content of baryons, atoms, and bulk 
matter [113]. For hydrogen and antihydrogen, a modern 
version of the argument could proceed by first noting 
that the quarks in hydrogen contain only about 10% of 
the mass with most of the remainder contained in the 
gluon and sea binding, and then concluding that since 
the binding forces are comparable for hydrogen and an- 
tihydrogen their gravitational response cannot differ by 
more than about 10%. This type of reasoning implic- 
itly assumes that the gravitational response of a body is 
determined by its mass and hence also by its binding en- 
ergy. However, as shown gencrically in Sec. VIII A, the 
coefficient (a^)T in the IPM leads to a correction to the 



gravitational force that is independent of mass and can 
vary with flavor. Indeed, the binding forces are largely 
conventional in the IPM, and the gravitational responses 
of hydrogen and antihydrogen are primarily determined 
by the flavor content of the valence particles. It is even 
conceivable in principle that a large gravity effect could 
be associated purely with the positron, as occurs in the 
IPM when only (a^ffOT is nonzero and satisfies the con- 
dition (169). A careful treatment of this issue in the 
IPM would require consideration of radiative effects in- 
volving (a^fj)T, (c'")tt, and other SME coefficients for 
Lorentz violation [18, 114], perhaps imposing the condi- 
tion (169) only after renormalization. In any case, the 
essential points illustrated with the IPM remain valid: 
the gravitational response of a body can be independent 
of mass, can vary with flavor, and can differ between 
particles and antiparticles. 

The gravitational response of antimatter could in 
principle also be restricted by the results of experi- 
ments studying kaons [115] and other neutral- meson sys- 
tems, which are natural interferometers mixing strong- 
interaction particle and antiparticle eigenstates via weak- 
interaction effects. When analyzed in the context of 
the SME in Minkowski spacetime, neutral-meson mix- 
ing places tight constraints on certain differences of the 
coefficients (ojf^)^ for w ranging over several quark fla- 
vors [53, 54]. However, these constraints have no dom- 
inant implications for leptons or for baryons, which in- 
volve three valence quarks rather than a quark and an 
antiquark as in mesons. Moreover, the neutral-meson 
constraints necessarily involve valence s, c, and b quarks, 
which are largely irrelevant for protons and neutrons. In 
the presence of gravitational interactions, the same line 
of reasoning holds, with the flavor dependence of Lorentz 
and CPT violation leading to the conclusion that the 
IPM evades restrictions from meson oscillations. 

We can also use the IPM to illustrate a type of con- 
straint on more realistic model building arising from 
the extensive searches for Lorentz and CPT violation in 
Minkowski spacetime. The key point is that the mixing 
of Lorentz- violating operators under rotations and boosts 
can imply indirect limits on some coefficients. In the 
IPM, for example, the coefficient (a^ff)T is unobservable 
in Minkowski spacetime, as discussed in Sec. II C, but 
certain nongravitational experiments could in principle 
obtain boost-suppressed sensitivity to (c™)tt for some 
w via measurements of the coefficients {c^)jk- As one 
illustration, a measurement with a Cs-Rb double foun- 
tain clock over a total of five weeks in the spring and fall 
of 2005 achieved a sensitivity of parts in 10^^ on some 
combinations of the coefficients {'c?')jk [116]. This sug- 
gests that continuing an experiment of this type over a 
longer period could attain parts in 10^^ on the coefficient 
(c^)tt by analyzing the data allowing for the Earth's or- 
bital boost ~ 10~^. Similarly, a careful analysis of 
multiple searches for Lorentz violation involving the elec- 
tron sector could be used to measure (c*^)tt at the level 
of parts in 10^^ [117]. Although these types of nongrav- 
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itational studies remain to be performed, they could in 
principle place experimental limits on the magnitude of 
the anomalous gravitational response of aiitihydrogen in 
the IPM and possibly also in more realistic models. We 
remark in passing that these kinds of constraints nonethe- 
less leave considerable room for realistic model building, 
in particular when operators of arbitrary dimension are 
incorporated in the framework [118]. 



C. Tests with matter beyond the first generation 

Most studies of fermion-gravity couplings to date have 
involved particles from the first generation of the SM. 
However, the SME coefficients for Lorentz and CPT 
violation can differ between sectors, so investigations 
of higher-generation matter-gravity couplings are of in- 
dependent interest. Since fermion masses and hence 
fermion-gravity couplings typically increase with the gen- 
eration, it is conceivable that an unconventional gravita- 
tional coupling may be more readily identified in grav- 
itational tests with higher-generation matter. Compar- 
atively few results exist for the coefficients (aJfgOp and 
(c'")/ij/ for particles w beyond the first generation [2], so 
there is considerable room for measurements of effects 
involving gravity couplings. 

The comparatively long lifetime of the muon makes it 
an interesting candidate for gravitational tests of Lorentz 
violation with second-generation particles. Several muon 
coefficients for Lorentz and CPT violation have already 
been measured [119], but the sensitivities are largely lim- 
ited to spin-dependent effects. Measurements of Lorentz- 
violating gravitational couplings of the muon could be 
achieved via muonium interferometry, with an estimated 
initial reach of 10% [120]. Interferometry with muonic 
hydrogen may also be possible [121]. In principle, these 
experiments could yield first measurements of some com- 
ponents of the coefficients ia^g)fj. and (ji^)nu in the muon 
sector. In particular, free-fall WEP tests using muo- 
nium interferometry to search for Lorentz and CPT vio- 
lation off'er the prospect of direct sensitivity to the co- 
efficients (agff)T and {c^)tt- In contrast, performing 
free-fall gravimetric tests with muonium interferometry 
is unlikely to be useful in the near future because the 
dominant signals appear at annual frequencies and are 
suppressed by the boost V®. 

Consider for definiteness a free-fall WEP experiment 
comparing the gravitational acceleration of muonium 
with that of neutral matter N . Muonium is a bound 
system containing an antimuon and an electron, so its 
spin-independent Lorentz-violating gravitational proper- 
ties are determined by the coefficients — (aeff)A" {'^^)t^i'7 
(^cfiOp (c'^)mi^- Following the line of reasoning in Sec. 
VI B, we find that the dominant observable combination 
of coefficients for CPT-odd effects in a free-fall WEP ex- 
periment is 



where is the mass of N and (a^f is its effective co- 
efficient for Lorentz and CPT violation. Assuming TV is 
composed of first-generation particles, the existing con- 
straints on coefficients [2] imply that for most models it 
is a good approximation to neglect all but the first term 
on the right-hand side of this equation. For CPT-even ef- 
fects, the relevant observable combination of coefficients 
is 



(173) 



-(Ogg)^ for the muon coefficient for Lorentz 



N 



(172) 



Again, only the first term is likely to be significant in 
practice. Similar expressions hold for muonic hydrogen, 
with the replacements e — )■ p for the superscripts and 

(«eff)A' 
violation. 

Searches for Lorentz-violating gravitational couplings 
of other second- and third-generation particles could also 
be countenanced. The typically short lifetimes of these 
particles can in principle be overcome by boosting, so ac- 
celerator experiments are likely to provide the best lab- 
oratory prospects. Studying the gravitational infall of 
particles of extraterrestrial origin in the context of free- 
fall WEP searches might be a source of additional con- 
straints. 

The physical mixing of uncharged particles of differ- 
ent ffavors w offers an interesting alternative method to 
achieve sensitivity to the coefficients (a^)^. Examples 
already yielding SME constraints on Lorentz and CPT 
violation include the interferometric oscillations of neu- 
tral mesons [53, 54] and of neutrinos [55]. Particle mixing 
implies nondiagonal terms in the propagator matrix, so 
field redefinitions of the type (14) cannot be used to re- 
move the coefficients (a^f^)^ from the theory. Differences 
between the coefficients {li^g)^ then become observable 
even in Minkowski spacetime, offering sensitivity to ef- 
fects that would otherwise be undetectable. For instance, 
the recent observation of anomalous CP-violating effects 
in S-meson oscillations [122] could originate in one or 
more nonzero coefficients {o^gjT for Lorentz and CPT 
violation in the quark sector, since these control CP- 
odd but T-cven operators that contribute to the effective 
hamiltonian for the mixing [54] . Spin- independent CPT- 
odd Lorentz violation involving coefficients such as (a'")^ 
could also underlie the observed baryon asymmetry in the 
Universe [123]. 

In terms of the perturbative counting scheme of Sec. 
II D, the existing SME studies using neutral-meson and 
neutrino oscillations lie at 0(1,0). Incorporating leading- 
order gravitational couplings along the lines in this pa- 
per would introduce 0(1,1) oscillation effects, including 
species-dependent modifications of the meson or neu- 
trino trajectories with characteristic time dependences 
similar to the WEP-violating effects discussed in Sec. 
VII D. Possible 0(1,0) contributions to the oscillations 
can be distinguished from 0(1,1) ones via the depen- 
dences on energy, baseline, flavor, and time. The advent 
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of neutrino-oscillation experiments with long and very 
long baselines of order 100-1000 km and corresponding 
changes in gravitational potential along the beams may 
offer particularly interesting options for free-fall WEP 
tests of Lorentz and CPT violation of this type. A 
detailed consideration of these possibilities would be a 
worthwhile subject for future investigation. 



X. SOLAR-SYSTEM TESTS 

Studies of the motion of bodies within the solar sys- 
tem provide an important source of information about 
gravitational couplings to matter. In this section, we in- 
vestigate the effects of nonzero coefficients (aoff )^ and 
c^^ for Lorentz violation in two solar-system contexts: 
lunar and satellite laser ranging, and perihelion preces- 
sion. The analysis here neglects effects that act merely 
to scale the mass of the gravitational source. These are 
unobservable using solar-system observations alone, but 
they may be detectable in combined measurements using 
photon tests. This latter issue is revisited in Sec. XI. 



A. Lunar and satellite laser ranging 

Lunar and satellite laser ranging provides a sensitive 
test of gravitational physics. The relevant orbital per- 
turbations to the motion of a satellite orbiting the Earth 



that arise from nonzero Lorentz violation in the pure- 
gravity sector of the minimal SME have been established 
[7] and used to constrain some of the coefficients 's^^ [4] . 
Here, we seek to extend these results to include dominant 
effects from nonzero coefficients {aeg)^ and c^^. 

Where possible in this subsection, we follow the con- 
ventions of Ref. [7] . A summary of our notation is given 
in Table XII. The flavor dependence of the matter-gravity 
couplings leads to composition-dependent factors in some 
of the equations to follow. To simplify these expressions, 
it is useful to define the eight combinations 
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Table XII. Notation for laser-ranging tests. 



Quantity 



Definition 
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m2 
N^ 
M = 
5m 



■ + TO2 

TO2 — mi 



m. 

M: 



o 



N, 




r( 

ri 

= 

RJ = 
^® = 
vJ = 



r{ -ri = {x,y,z) 
■ {mir'i + m2r2)/M 
- ^/GnMq/R3 
v(-v( = dr-' /dT 
■■ {miVi + m2V2)/M 



satellite mass 

number of particles of species w in the satellite 
Earth mass 

number of particles of species w in the Earth 
total Earth-satellite mass 
Earth-satellite mass difference 
mass of the nth perturbing body 
Sun mass 

number of particles of species w in the Sun 
satellite position 
Earth position 

Earth-satellite separation, of magnitude r = \ fi — r2| 
position of Newton center of mass for Earth-satellite system 
mean Earth orbital frequency 
relative Earth-satellite velocity 

velocity of Newton center of mass for Earth-satellite system 



The primary observable in laser-ranging tests is the 
coordinate acceleration a^g of the relative Earth-satellite 
separation. Working in the Sun-centered frame, we can 
obtain this acceleration from the equation of motion (78). 
The relevant contributions to the coefficient and metric 



fluctuations from {aes)fj. and c^,y can be found in Sec. 
IV C, while those from Sj^^, are given in Ref. [7]. 

Incorporating perturbative effects of other bodies in- 
cluding the Sun, the coordinate acceleration can be writ- 
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ten 



-12^.1 



^LV + — 



(175) 



The first three terms in this expression involve effects 
independent of Lorentz violation. They represent the 
acceleration due to the Newton gravitational field of 
the Earth-satellite system, the Newton tidal quadrupole 
term, and the quadrupole moment of the Earth, respec- 
tively. Their explicit form is given in Ref. [7]. The 
leading Lorcntz-violating contributions to the accelera- 
tion are represented by the fourth term aLy. This term 
can itself be split into four pieces, 



The first two terms are the ones of interest in the present 
work and are discussed below. The last two depend on 
the coefficient s^^, with a^gg arising from the Earth- 
satellite system and a^tj^ai involving perturbations due 
to other bodies. The explicit form of these two quantities 
is provided in Ref. [7]. 



The term a^^^^ -gg in Eq. (176) provides the Lorentz- 
violating acceleration of the Earth-satellite system from 
the matter-gravity couplings (cioff )^ and c^,^. It takes the 
form 



^aeffjC.ES 



W 

- 2n^m^{c^\TK)V''r' + 2n^m^{c'")(TK)v''r' + 2n^wy'n-"' {c^\TK)VLr 

- 2(n^ - 2n^)m^rj'''{c^)(^TK)VLT 



(177) 



In principle, a^^^^ ^ j^g also acquires contributions proportional to R^ijO, but these are neglected here because they are 
typically suppressed compared to effects proportional to V'' and . 

In Eq. (176), the term a^^^^ c tidal contains the Lorentz-violating tidal acceleration involving (aeff);^ and 0^^, which 
arises from perturbing bodies. When the satellite is taken as the Moon, the dominant tidal contributions are due to 
the Sun and can be written 



tidal = E { (2"«ff)T + m^{c-)TT) - ^n5"a«ff)T] (Sr^RLR' - r^) 

w ^ 



M 



R' 



m© m© \ M J \ J 

2<ry'^^14(c-)(TK) + 2<77^^^;^(c-)(TK) 

+ n^F^(c«')(TL) + n^t;^(c'")(Ti)] {j.r''' RmR^ - r^) }• (178) 



- 2 



M 



If instead the satellite is artificial, then there are tidal 
effects from both the Sun and the Moon. However, these 
are suppressed relative to the Earth-satellite acceleration 
(177). 



The Lorentz-violating coordinate accelerations given 
by Eqs. (177) and (178) exhibit some interesting fea- 
tures. The first two terms in Eq. (177) and the first 
term in Eq. (178) are composition-dependent scalings of 
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the corresponding Newton accelerations. These terms 
are therefore detectable only by comparison to results 
obtained using satellites of different compositions. Also, 
unlike the contributions a^ti^^i obtained in Ref. [7], here 
the tidal acceleration (178) from the Sun on the Moon- 
Earth system involves nontrivial WEP violations because 
the Moon and the Earth fall differently towards the Sun 
when the coefficients (a^g)/j and (c"')^,y are nonzero. It is 
also interesting to note that the tidal acceleration (178) 
contains contributions at PN0(2) that are independent 
of r-^ and hence are enhanced at this order relative to 
other contributions by a factor of R/r. This too is a 
consequence of the WEP violations arising from (aeff),^. 
Similar terms appear at PN0(3) as well. 

A typical experiment measures the time of flight for 
laser photons to travel from the Earth to a reflector on 
the satellite and back. To analyze the results, the laser- 
ranging data can be fitted by incorporating Eq. (175) and 
other conventional perturbing effects into an appropriate 
modeling code. An alternative approach is to perform an 
analytical perturbative expansion along the lines of the 
one performed for the contributions in Ref. [7] and 
then match to the data. This latter method is adopted 
in Ref. [4] to constrain combinations of the coefficients 

In the present context, we can obtain crude estimates 
of sensitivities to (a^)/n and attainable in lunar 

laser ranging via either of these procedures, by using 
term- by-term comparison of the accelerations (177) and 
(178) to the accelerations a^_ES ^"^^ ^tidai obtained for 
the coefficient s^iv in Ref. [7] . With the precision already 
achieved in lunar laser ranging [124], we thereby find es- 
timated sensitivities at parts in 10^*^ to combinations of 
^^){JK) and sjK, and parts in 10^ to various combina- 
tions of (ajfff)j, {(?"){TJ) and stj- Actual measurements 
at roughly these levels can be expected to result from 
a reanalysis of existing data. A significant further im- 
provement is likely to be possible using data from the 
Apache Point Observatory Lunar Laser-Ranging Opera- 
tion (APOLLO) [125]. Assuming that millimeter rang- 
ing is achieved as expected and disregarding probable 
substantially improved statistics, we anticipate competi- 
tive estimated sensitivities of 10~^ GeV on various com- 
binations of a{a'^g)x and a{a'^^)Y+Zi and a sensitivity 
of 10^^ on (c")(7-j), where the notation of Eq. (148) is 
used and the dependence on the coefficients s^v has been 
omitted for simplicity. 

Ranging to artificial satellites with orbit orientations 
different from that of the Moon can yield sensitivity to 
additional independent linear combinations of (a^)p and 
{c^)tiu- Typically, the reach of satellite ranging is ex- 
pected to bc! about an order of magnitude less than lunar 
laser ranging. Other possibilities for gravitational tests of 
Lorentz violation include ranging to objects orbiting bod- 
ies other than the Earth. For example, the time variation 
G^dGj^ /dt of the Newton gravitational constant has 
been constrained by ranging data to the Viking landers 
on Mars, to the Mariner 9 spacecraft orbiting Mars, and 



to other bodies including the Moon [126]. These studies 
primarily seek secular changes in the gravitational force. 
Although secular changes in coupling constants can re- 
sult from Lorentz violation [127], the signals of interest in 
the present context are periodic. Reanalysis of existing 
data to seek periodic effects in Gl^}dGisT/dt would yield 
sensitivities to Lorentz violation estimated to be some- 
what less than lunar laser ranging but involving different 
combinations of (a^)/j and (c^)^;y. 

We conclude this subsection with some comments 
about the coordinate location J?-^ of the center of mass 
of an Earth-satellite system. Boost invariance normally 
ensures this location is fixed, but the presence of Lorentz 
violation means it can be time dependent, although the 
effect may be unobservable via laser ranging. Neglecting 
the effects of other bodies and working at PN0(3), the 
Lorentz- violating contributions to the equation of motion 
for the center of mass of the Earth-satellite system can 
be written as the sum 

R' D + Ri + Ri (179) 

of contributions from (a^),^, ic^)^vi and s^. Explicitly, 
we find 

^^^^^^^2G.nW-a^,)...r- ^ (180) 

w 

which contains only PN0(3) effects involving the internal 
motion of the system. The second term in Eq. (179) is 



w 
X 
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+r?-^^(c'")(T^)(n^yL+n>r)r^ 
+(c'")(TK)(n^y^ + n^^;^)r-^ 

+{^\TK){n^V' + ntv')r''\ (181) 

The first two terms are at PN0(2) and reflect the modi- 
fication of the effective Newton inertial mass in the pres- 
ence of nonzero c^v, while the remaining terms are at 
PN0(3). The ones proportional to arise as a result 
of the system boost in the Sun-centered frame, and those 
proportional to f are due to the internal motion of the 
system. The last term in Eq. (179) is 



^ = 



J„K 



3?7 StkVlt — stkv r 



(182) 

which again consists only of PN0(3) effects proportional 
to the internal motion of the system. Note that all these 
contributions introduce an oscillatory motion for the cen- 
ter of mass, and their presence is required by momentum 
conservation. 
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B. Perihelion precession 

The presence of nonzero coefficients (acir)/j, and c^i^ for 
Lorcntz violation leads to corrections to the motion of a 
test body in a gravitational field. These corrections can 
be calculated from the equation of motion (78) and from 
the expressions for the coefficient and metric fluctuations 
given in Sec. JVC. In this subsection, we determine the 
effect of nonzero {aes)fj. and c^i/ on the perihelion pre- 
cession for planetary orbits. We follow the treatment of 
Ref. [7], which obtains the perihelion shift arising from 
nonzero SME coefficients s^^ . Our notation matches that 
of Table XII and Eq. (174) in Sec. X A, with the labels 1 
and 2 representing the planet and Sun, respectively. 

The derivation of the perihelion precession used here 
relies on the method of osculating elements [128], in 
which the instantaneous motion of the planet is treated 
as part of an ellipse. The ellipse is characterized using 
the standard Kepler orbital elements, and the motion of 
the planet is described by specifying them as a function 
of time. The relevant orbital elements in the present case 
are the angle oj between the line of ascending nodes and 
the semimajor axis of the ellipse, the longitude fl of the 
ascending node, and the inclination i with respect to the 
ecliptic. These are specified in the reference coordinate 
system, which can be taken as the Sun-centered frame 
for the planetary orbits considered here. More generally, 
the reference frame is related to the Sun-centered frame 
by a rotation and possibly a boost, as discussed in Sec. 
V E 5 of Ref. [7] . The physical quantity relevant for the 
perihelion precession is the change per period Aw of the 
perihelion angle u) with respect to the equinox. In terms 
of the basic orbital elements, w can be expressed as 



n cosz. 



(183) 



For the cases of interest here, the angle i can be assumed 
small. 

The secular changes in the orbital elements arising 
from {aes)fj, and Cfj,^, can be obtained by considering the 
relative acceleration of the planet and the Sun, which has 
the form 
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The unperturbed elhpse is given as the solution rb of the 
Kepler-type equation 



df^ 



^ ^ [m + 2n^a(a:fff)o + n^m^{c^U 



(185) 



This shows that the frequency n and semimajor axis a of 
the unperturbed elliptic motion are related according to 



= Gjv^ [M+2nX<ff)o+nfm"'(c'")oo]. (186) 



Note that the right-hand side of this expression depends 

on the composition of the planet and the Sun. 

The orientation of the orbit can be specified using three 
unit vectors fc, P, and Q. The first is chosen perpendic- 
ular to the orbit, the second points from the focus to the 
perihelion, and the third completes the orthonormal set. 
Their explicit form in terms of orbital elements is given 
in Eq. (116) of Ref. [7]. In terms of this basis set, the 
unperturbed elliptical orbit can be expressed as 
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1 + e cos / 



(Pcos/-|-Qsin/), (187) 



where e is the eccentricity and / is the true anomaly. 

The perturbing acceleration a'^ consists of the terms 
in Eq. (184) that are absent from Eq. (185), 
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The time dependence of the orbital elements can be ex- 
tracted from this equation via the method of osculating 
elements. The general procedure is to insert the unper- 
turbed solution (187) for r into the expression (188) for 
a'^ , to project the result as desired, and to integrate over 
the true anomaly. 

To obtain the perihelion precession, the final results 
for the orbital elements uj and 17 must be combined ac- 
cording to Eq. (183). After some calculation, we obtain 
the expression 
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for the shift in the perihelion per orbit. Here, the sub- 
scripts P and Q on the coefficients for Lorentz violation 
indicate projections along the directions P and Q, re- 
spectively. The qua ntity e is the eccentricity function, 
defined by e = 1 — \/l — e^. 

The result (189) reveals that the perihelion precession 
depends on the orbit orientation through the projections 
of the coefficients {aes)n and c^^ along the directions P 
and Q. Also, the factors scaling the coefficients in Eq. 
(189) vary with the composition of the orbiting body. 
This means that the orbits of different planets or, more 
generally, different satellites are affected by different lin- 
ear combinations of coefficients for Lorentz violation. It 
is therefore valuable to consider data from multiple sys- 
tems so that independent measurements can be obtained. 
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To illustrate the sensitivities that can be achieved, we 
consider explicitly the perihelion precessions of Mercury 
and of the Earth. Substituting the relevant orbital data 
for the two planets into Eq. (189) in turn, taking the plan- 
etary mass as small compared to the solar mass m© , and 
incorporating the results for the coefficients s^,y obtained 
in Eq. (190) of Ref. [7], we find the overall perihelion 
shifts u)^ of Mercury and u>^ of the Earth are given in 
units of arc seconds per century C by the expressions 

7 X 107"_ 



+ 



C 

w 

2 X 10^" 



/ /V 

5:(3xio-3-j«,)^-^(.-)^ 



We 



C 



c 



niQ - TOe 



r)e). 

(190) 



The combinations of coefficients for Lorentz violation ap- 
pearing in these equations are defined as 



[(5"')QQ-r)pp] -6 xl0-3(c-)(0Q), 
[(c"')QQ-(c-)pp]-5xl0-2(c-)(0Q), 
(spp - sqq) - 6 X 10"-^S(oQ), 



(spp - sqq) - 5 X 10 ^S(oQ). 



(191) 



Note that the subscripts P, Q here represent projections 
that differ for Mercury and the Earth. 

The chemical composition of the Sun is believed to 
be over 70% hydrogen and about 27% helium by mass 
[129]. The factors in Eq. (190) that depend on the so- 
lar composition can therefore be estimated as Nq/ttiq = 
Nl/niQ ~ 0.9 GeV-i and iVg/m© ~ 0.1 GeV"!. As 
can be seen from Eq. (191), these factors suffice for plac- 
ing approximate bounds on the coefficients (a^)/n from 
knowledge of the perihelion precessions. The composi- 
tion of Mercury is believed to be about 70% iron and 
about 30% rocky material [129], so the analogous ratios 
for Mercury are roughly N^/m^ = N^/m^ ~ 0.4 GeV""^ 

and N^/m^ ~ 0.6 GeV~^. For the Earth, using Ref. 

[46] and following the discussion of Sec. IV A 2, we find 

iV|j/m© = A^^m© « N^/rriQ ~ 0.5 GeV'^. However, 
for the approximate bounds obtained below on the co- 
efficients c^y, it suffices that the composition-dependent 
factors for the planets are of order 10~^ GeV"'^. 

We are now in a position to place constraints on some 
combinations of the coefficients (a^)^, c^j,, and s^^i/ by 
adopting the established error bars in the existing data 
for perihelion shifts. These error bars are 0.043" for 



Mercury and 0.4" C'^ for the Earth [8, 130]. Taking the 
error bars to be upper bounds on the perihelion shifts in 
Eq. (190), we obtain the order-of-magnitude constraints 

-10-''(c^)^ - 10-^(cf )^ - 10-i(c")^l < 10-9 GeV, 

[3© + lO-^[iat«)e + (a^ff)e] + 10~'«ff)e 
-10-^(c«)© - 10-i(cf )e - 10-i(c")©l < 10-« GeV. 

(192) 

Assuming a model with nonzero coefficients {a^f{)ij, only, 
this yields the approximate constraints 

Ka^ff)^ + (a^ffO^^ + O.Ua'eftO^I < lO"*^ GeV, 
m«)e + (Ksh + 0-l(ae"ff)®l < 10-*^ GeV. (193) 

Similarly, assuming a model with nonzero coefficients 
only and making use of existing limits on c^^ for 
protons and electrons [2], we obtain the approximate con- 
straints 



5")^l 



< 10" 



(5")el<10- 



(194) 



A careful reanalysis of the existing data for multiple bod- 
ies in the solar system could yield sharper sensitivities. 

The result (193) represents first constraints on the spa- 
tial coefficients a(a^) j. A sense of the maximal attained 
sensitivity to the nine components in Q:(a^)j can be ob- 
tained by taking each component in turn to be the only 
nonzero one. Extracting these sensitivities requires the 
explicit form of the vectors and Q©. The relevant or- 
bital elements in heliocentric coordinates arc uj^ ^ 29°, 
~ 48°, ~ 7° and ~ 103°, VL^ = 0° = 0°. 

Converting to the Sun-centered frame using a counter- 
clockwise rotation by 77 ~ 23.5° about the X axis yields 

~ -0.97ex-f 0.15ey + 0.18ez, 
Q© ~ -0.97ex-0.21ey-0.10ez. (195) 

Taking each component Q;(a^) j as the only nonzero coef- 
ficient in turn yields the order-of-magnitude sensitivities 

HKs)x\ < 10-^ GeV, 

\aiats)Y\MKs)z\ < 10-^ GeV, 

Hals)x\ < 10-^ GeV, 

\a{a^MHaPf,)z\ < IQ-' GeV, 

|a(a^ff)x| < 10-'^ GeV, 

\a{a^s)Y\M<s)z\ < 10-4 GeV. (196) 

These results are the maximal sensitivities achieved to 
date on the coefficients a(a^)j. 
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XI. PHOTON TESTS 



The wave 4- vector can be linearized as 



In this penultimate section, we consider searches for 
gravitational Lorentz violation involving the trajectories 
of photons. With the coordinate choice (23) adopted in 
this work, photons follow null geodesies. The signals of 
interest therefore arise from the modifications to the met- 
ric, which are associated with Lorentz-violating matter- 
gravity couplings of the source body and in certain cases 
also of the clocks and rods used for measurements. 

Photon tests for Lorentz violation involving the coef- 
ficients Sfjiy ill the purci-gravity sector of the SME have 
been studied in Refs. [7, 131]. Here, this analysis is ex- 
tended to include the matter-sector coefficients (aJJ^)^ 
and i(f")fiij. The treatment and notation of Rcf. [131] is 
adopted where possible. Some quantities relevant for the 
analysis are listed in Table XIII. 

Table XIII. Notation for photon tests. 



)uaut,it \- 



Dcliuitiou 



Xe = {tE,rE) 

te 

a;^ = {tp,rp) 
rp 

R = fp - te 
R = R/R 
R= \R\ 
V =i' 
b 



R^rp ■ R 



m 

Ip = fp ■ R 
-Ie = fh-R 

TE 

Tp 

Ue = dx'^/dTE 
= dx'^ /drp 
V = dvE I dt 
w = dfp/dt 
Ue 

Up 



coordinates of event E 
magnitude of 
coordinates of event P 
magnitude of fp 
zeroth-order light trajectory 
unit vector along R 
magnitude of R 
impact-parameter vector 
magnitude of b 
mass of source body 
A at P 
A at £ 

proper time of E 
proper time of P 
4- velocity of E 
4-velocity of P 
3-vclocity of E 
3-vclocity of P 
frequency at E 
frequency at P 



In what follows, we consider various effects on a light 
signal as it travels from an emission event E to a space- 
time point P located near a massive body. The light path 



can be specified parametrically as x^ 



'(A), where A 



is the path parameter. The wave 4-vector p^^ of the ray 
tangent to the path is 



dxi^ 
rfA"' 



(197) 



and it obeys the conditions 



= — r'' 



dx^" 
'dX 

p'^P^gi.u = 0. 



(198) 



: P"" + Spl", 



(199) 



where the first term is the zeroth-order wave vector and 
the second term contains gravitational corrections. Our 
interest here lies in the 0(1,1) contributions to Jp'^. 
The basic procedure is to insert the modifications (84) 
and (88) of the metric arising from matter-sector effects 
into the general expressions obtained in Ref. [131]. We 
consider in turn Lorentz-violating contributions to the 
Shapiro time delay, to the gravitational Doppler shift, to 
the gravitational redshift, and to the null redshift, and 
we compare the results to the effective mass of a gravita- 
tional source as measured by orbital tests. We also offer 
some comments about the implications of the results for 
various experiments. 



A. Shapiro time delay 

In this subsection, we obtain the Lorentz-violating 
modifications to the Shapiro time delay of a light signal 
as it passes from a source to a detector in the presence of 
a massive body such as the Sun. The one-way time delay 
tp — tE can be determined by integrating 6p^ along the 
path and applying the null condition, 



tp — tE — R 



h^j.^jfp'^dX. 



(200) 



Inserting the Lorentz-violating metric modifications (84) 
and (88) and integrating, we find the delay can be written 

in the form 

tp — ts = R+{tp — tE)GR 

+ {ip-tE)a^,,.-c + {tp-tE)-s. (201) 

Here, R is the zeroth-order time difference. The second 
term is the standard OR contribution, which at 0(0,1) 
and PN0(2) takes the form 



{tp - tE)GR = 2Givm^ In 



rE + rp + R 
rE + rp — R 



(202) 



The third term of Eq. (201) consists of contributions 
from Lorentz-violating matter-gravity couplings associ- 
ated with the source body S. At 0(1,1) and PN0(2), 
these contributions are 



{tp — tE)a^ 



\ I I t lit r 



X In 



Giva((«eff)o + (a^ff),i?^)(l^ + I^ 



my m 
rs + rp + R 



rE + rp 



/'77S 



R 



h 



rE rp 



-Gjva(afff)fe6'= 



rE - rp 
rErp 



(203) 
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The final term in Eq. (201) arises from gravitational 
Lorcntz violation involving the coefficient s^j, and is given 
in Ref. [131]. 

In typical time-delay measurements, an observer emits 
a light signal at E that is reflected at the spacetime point 
P and subsequently detected by the observer at E' . The 
round-trip coordinate travel time At, which is related to 
the measured proper time Ate by the factor drE/dt, can 
be written to 0(1,1) and PN0(2) as 

At = 2R{l + v^-v-R) + {At)GR + {At)a^,,^^ + {Ath. 

(204) 

The zeroth-ordcr term in this expression incorporates 
Lorentz-violating corrections to the trajectory of the 
emitter, which here can depend on particle species. These 
can in principle be determined by modeling the relevant 
orbits along the lines of the treatment in Sec. X and Ref. 
[7]. The second term in Eq. (204) contains the leading 
GR corrections, 



AtcR = AGNmP In 



te +rp + R 
TE + rp - R 



(205) 



The third term in Eq. (204) contains the leading contri- 
butions from nonzero {aes)n and Cij,v, 



AU 



AGntu^ 



-(4f)o + (c^) 



00 



X In 



TE + rp 



R 



-2GNa{als)o 



\rE +rp - R 

Ie 



Te 



lp\ 
rp) 



(206) 



The last term of Eq. (204) contains corrections involving 
the coefficient and is given in Ref. [131]. Note that 
contributions from the coefficients (aeff)j and soj cancel 
in the round-trip expression, a result that can be traced 
to the parity-odd nature of the corresponding Lorentz- 
violating operators. Note also that the time-delay signal 
changes over two relevant time scales, the conjunction 
time h/v and the typically longer orbital time r/v, which 
enables separation of the zeroth-order and gravitational 
effects. 

The dominant Lorentz-violating corrections to At are 
proportional to the logarithm in Eq. (206). The primary 
effect of the Lorentz-violating matter-gravity couplings is 
therefore to scale the factor of Gniti^ in the usual GR 
time delay (205). The scaling can be interpreted as an 
effective value {GnM)td for the source body relevant for 
time-delay tests. 



(GnM) 



TD 



Gntti^ (l + ^(afff)o + {c^)oo + soo) 



(207) 

This scaling is unobservable in time-delay tests alone. 
However, we show in what follows that other tests can 
yield different effective values of Gjvrn^, so suitable com- 
parisons can reveal signals for Lorcntz violation. This 
prospect is considered in Sec. XI E below. 



B. Gravitational Doppler shift 

When light passes near a massive body, it suffers a 
frequency shift as well as a time delay. In this and the 
subsequent subsections, we consider the corrections to 
the frequency shift due to the matter-sector coefficients 
(oeff)^ and c^„. 

The relevant quantity is the ratio of frequencies ob- 
served at the two events E and P, 



Up^ 

ve 



{u''p^,)e 

At PN0(3), this can be written as 



i^p 
Ve 



l-tZ;-i? 
l-v-R 



1+ — 



(208) 



(209) 



Here, the term labeled g contains gravitational effects 
involving both the Doppler shift and the redshift, 



(210) 



This subsection treats the gravitational Doppler shift, 
while the redshift effects are discussed in the next sub- 
section. 

Corrections to the gravitational Doppler shift 
{vp/ve)'ds depending on the coefficients (aeff)^, c^;/, 
and Spjy can be obtained by inserting into Eq. (31) of 
Ref. [131] the modifications to the metric from Eqs. (84) 
and (88), along with those due to s^j, given in Ref. [7]. 
Near conjunction, we find that the dominant effects take 
the form 



(l + ^(a,Sff)o + ^(aeU-^^- 



00 



dh 
di 



m 

i^p 



(211) 



where the last term contains the contributions from s^u 
found in Ref. [131]. 

Typical searches measure the round-trip frequency 

shift, 



dv 



SGMm^ 



(l-H^(af,)o + (c^)oo)| 

(212) 



DS 



+ 



b 

V 



DS,: 



Note that the effects from parity-odd operators again 

cancel. The coefficients (aoff)o, cqo, and sqo associated 
with isotropic Lorentz violation in the chosen inertial 
frame act to scale the factor Gnttl^ in the usual expres- 
sion for the gravitational Doppler shift, leading to an 
effective value (G'ArM)Ds given by 



(GjvM)ds = Gntu^ (l 



m' 



S («efT)o 



(c')oo 



Soo , 



(213) 
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The scaling (213) is unobservable in Dopplcr-shift tests 
alone. This result for (G'jv-M)ds is identical in form to 
that of the time-delay value (GjvAf)TD in Eq. (207). 



C. Gravitational redshift 

The Lorentz-violating contributions to the term 
{vp/ve)ks in Eq- (210) for the gravitational redshift can 
be viewed as subdominant to the time delay or Doppler 
shift because they occur at the slow time scale. How- 
ever, in dedicated redshift measurements, the Lorentz- 
violating gravitational redshift can appear as the dom- 
inant effect. In this subsection, we discuss Lorentz- 
violating modifications to the usual gravitational redshift 
and effects in nuU-redshift tests. 

To place in context the results in this subsection, we 
note that clocks can be used to perform three distinct 
types of gravitational tests that are often convolved in 
the literature under the term 'redshift tests.' The first 
type, which measures the traditional gravitational red- 
shift, involves two clocks held at different gravitational 
potentials whose frequency is compared using light or 
some other signal passing between them. This type of 
test is discussed in Sec. XI C 1 below. The second type of 
test is called a miU-rcdshift test, and it involves monitor- 
ing the frequencies of two clocks of different composition 
as they move together through the gravitational poten- 
tial. This is discussed in Sec. XI C 2. The third kind of 
test involves synchronizing two clocks and then moving 
one of them around a closed path in the gravitational po- 
tential. The signal in this case is the accumiilated phase 
difference between the clocks. An example of this 'twin- 
paradox' redshift test is the free-fall gravimeter measure- 
ment with interferometers discussed in Sec. VII D. These 
three kinds of tests produce related signals in GR. How- 
ever, they can yield distinct sensitivities in a more gen- 
eral context such as the SME, as is demonstrated in what 
follows. 



1. Modified redshift 

The term {vp/ve)ks in Eq. (210) for the gravitational 
redshift can be understood as the product 



(214) 



of the factors relating proper and coordinate times for 
the clocks at the two points E and P. Each factor is de- 
termined by the dispersion relation for the corresponding 
clock, which depends on coefficients for Lorentz violation 
via its material composition and on the Lorentz violation 
associated with the gravitational field. 

For simplicity in what follows, we assume the sending 
and receiving clocks arc identical. This eliminates the 
need to consider 0(1,0) effects, which have been sought in 



numerous clock-comparison experiments performed with 
both clocks at the same gravitational potential [2]. To 
the order at which we work, the redshift can then be 
expanded as 



(0,1) 



(1,1)S 



'^p 



(1,1)T 



(215) 

where the term at 0(0,1) is the conventional redshift, the 
term at 0(1,1) labeled by S contains Lorentz-violating 
corrections from the gravitational source, and the last 
term labeled by T involves 0(1,1) contributions from the 
clocks. For our present purposes, it suffices to work at 
PN0(2). 

For an ideal clock, the Lorentz-violating contributions 
to the first two terms in Eq. (215) can be calculated by 
inserting into the usual redshift equation the modifica- 
tions (84) and (88) to the metric from the coefficients 
{aes)ti and c^^, along with the corrections from Ref. [7] 
involving the coefficients s^^. This gives 



(0,1) 



vp 



1,1)S 



/.(0,l)x /I.(l,l)\ 

(^00 )e- {ho6 )e 



1 



(216) 



Expanding to PN0(2) and keeping leading-order terms 
in Lorentz violation, we obtain the conventional PN0(2) 
result, 



(0,1) 



together with the correction 



iyp\ 
^eJ-hs 



(1,1)S 



Vp_ 



(217) 



RS,s 



The last term contains the contributions from s^u given 
in Ref. [131]. 

For the remaining term in Eq. (215), the situation is 
more complicated because the clock frequency must be 
calculated directly and typically depends on the struc- 
ture and composition of the clock. Moreover, although 
our interest is at 0(1,1), all three of the perturbative con- 
tributions 0(1,0), 0(0,1), 0(1,1) must be treated due to 
the appearance of cross terms in the calculation. For 
convenience, we can express the last term in Eq. (215) in 
the form 



vp \ 
"eJrs 



(1,1)T 



= GnITI ^clock 



(218) 



where ^ciock is a function of the coefficients for Lorentz 
violation associated with the clock. If the clock's ticking 
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rate is set by its inertial properties, as is the case for most 
atomic clocks, then ^ciock can be expected to depend on 
the coefficients (c"')^,y. If the clock's ticking rate depends 
intrinsically on the local gravitational acceleration, as oc- 
curs for a pendulum clock, then ^ciock can be expected to 
depend on the coefficients (a^^)^. In general, the value 
of ^ciock can depend on both sets of coefficients, 

^clock = ^clock((o^ff)/:t, {c^)iJ,v) ■ (219) 

The key point is that different clocks have different ^ciock 
according to the details of their construction and flavor 
content. 

Combining the above results, we see that the domi- 
nant Lorentz-violating effects for the gravitational red- 
shift can be represented as an effective value {GnM)-r,s 
implementing a scaling of Gniti^, in parallel with the 
results for the time delay and the gravitational Doppler 
shift. We obtain 

(GjvM)rs = GjvmS(l + |j(afff)o + (c^)oo + |soo 

+^ciock). (220) 

This represents an unobservable scaling in any particu- 
lar redshift test, but comparing redshift tests performed 
with different clocks could yield access to differences in 
Cciock- Moreover, the result for (Gjv-M)Rg differs from 
both the time-delay value (GjvM)td in Eq. (207) and the 
Doppler-shift value {GnM)-£,s in Eq. (213), so comparing 
results from different tests could yield independent sen- 
sitivities to (affj)o that are inaccessible in other searches 
with ordinary matter. This prospect is considered in Sec. 
XI E. 

We remark in passing that for certain special mod- 
els the observable redshift effects in {GnM)yis niay be 
hidden in WEP tests. A simple example is provided by 
the isotropic parachute model discussed in Sec. IX B. By 
virtue of Eq. (169), the effective inertial and gravitational 
masses in this model are equal for a test body made of 
ordinary matter, so no signals are observable in WEP 
tests. However, the presence of nonzero (c'")/ij/ implies 
a nonzero rescaling of (GjvAf)RS) which is observable by 
comparing to {GnM)t'd or {GnM)-£,s- Signals from this 
model could also arise in the nuU-redshift tests discussed 
in Sec. XI C 2 below. 

We conclude this subsection with an illustrative cal- 
culation of ^ciock for a simplified clock based on tran- 
sitions between the Bohr levels of hydrogen, for which 

we determine {i'p/i^e)'rs ^^ and ^ciock = Cn.Bohr assum- 
ing both the clocks and the gravitational source are at 
rest. This calculation is straightforward due to the spher- 
ical symmetry and the zero velocity, and also because 
a simple match exists between the zeroth-order hamil- 
tonian h^^'^^ and the kinetic contributions to the sum 
/j(i.o) ^ /j(o.i) -|- of the perturbative corrections pre- 
sented in Sec. Ill C. By matching these expressions, we 
find that the kinetic portion of the hamiltonian in the 
presence of gravity and Lorentz violation can be obtained 



from the zeroth-order one by the following simple replace- 
ments for the proton and electron mass: 

^ ^ ^(l-i/^oo + i(o^)oo + f(c^)ooM. 

^ ^(l-|/i00 + |(c^)00 + f (C^)OOM- 

(221) 

Also, the source term in the Maxwell equations is cor- 
rected by the vierbein determinant e, and the result 
can be obtained by a simple replacement for the proton 
charge, 

qP^^^ qP{l - /loo). (222) 
e 

It follows that the calculation of interest can be directly 

performed by implementing the above replacements in 
the standard result for the Bohr energy levels. This yields 

E ^ e(i- i/ioo 

+ ;;^FT^K(5^)oo +m^(c^)oo](i - i/^oo)). 

(223) 

The modification (218) to the gravitational redshift is 
therefore given by 

UJhs = " 3(m. + me) ("^'(^)°°+"^^^^)°°) 
X ( (224) 

V feTp J 

when the clock transitions are those of the Bohr levels of 
hydrogen. This implies the result 

^H.Bohr = - o, ,1 ,d riiF{^U + m^(c^)oo). (225) 

The value of ^ciock for a realistic clock can be obtained 
via calculation if the hamiltonian describing the clock is 
known. 

2. Null redshift 

Another Lorentz-violating signal can be accessed by 
comparing two clocks of different types as they explore 
the gravitational potential together. This type of mea- 
surement is called a null-redshift test [132]. 

Consider comparing the frequencies of two clocks A 
and B having different values ^ciock — and ^ciock = 
that are located at a point P with gravitational potential 
Hqq. The frequency ratio is given by 

{^y=[^+hiu-M](^^y (226) 
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where the superscript (0) denotes a frequency at a hy- 
pothetical zero gravitational potential hf^,y = 0. This 
frequency ratio depends inseparably on the potential /igo 
at point P and the ratio in zero potential. 

When the same two clocks are moved to a point Q at 
potential Hqq, the frequency ratio takes a new value. If 
the values and differ, then so do the frequency 
ratios at P and Q. The ratio of frequency ratios then 
shifts away from 1 and is given by 

{^y {'^y=^-'2iu- ^B){hg - h^,). (227) 

The shift is an observable, and it depends on the differ- 
ence A^AB = £,A — ^B and also on the potential difference 
between P and Q. 

For a gravitational source with hoo = ^GNm^/r at 
PN0(2), we obtain 

(^)7^)^.1-G.^s^a.^^^^. (228) 

Unlike the other Lorentz-violating photon effects dis- 
cussed here, all of which represent scalings of G^m^, 
this result is a qualitative change from conventional grav- 
ity. It is also strictly a gravitational effect, vanishing in 
Minkowski spacetime. 

Since the shift varies with spacetime position, it ex- 
hibits features analogous to violations of local position in- 
variance, which have been the subject of numerous stud- 
ies [8]. In the present case, these features arise from the 
Lorentz-violating flavor dependence of the clock mate- 
rial. Note also that the observable (228) contains the 
same information as the result of two separate redshift 
tests performed with diflFerent clocks but the same gravi- 
tational source. This can be verified by inspection of the 
effective value (G7v-^'^)rs in Eq- (220). Some relevant 
experiments are described in Sec. XI E. 

D. Comparison to effective orbital mass 

The preceding subsections reveal that the Lorentz- 
violating contributions to the Shapiro time delay, the 
gravitational Doppler shift, and the gravitational redshift 
arc all controlled by the effective value of Gat M for the 
gravitational source. In the context of the solar-system 
tests discussed in Sec. X, rescalings of GnM also occur 
but can be disregarded as unobservable. Here, we deter- 
mine the effective value of GjvM relevant to observations 
of orbiting bodies, {GnM)ob- 

For Lorentz violation involving the coefficients {aes)^ 
and Cf_iu, the secular changes in the orbital elements for 
the trajectory of an orbiting body are given by Eq. (184). 
The analogous result for the coefficients s^i, is given in 
Eq. (162) of Ref. [7]. Inspecting these equations, we can 
deduce the effective reduced mass of the source and test 
bodies and hence extract the effective value {GnM)ob- 



Making no additional assumptions about the masses of 
the source and test bodies, we flnd 

(GjvM)oB = GjvM(l-h^(afff)o + ^(aJff)o + fsoo 




To obtain an expression that is more readily compa- 
rable to the effective values of (GjyM) measured in pho- 
ton tests, we note that m/^ <^ under typical circum- 
stances. The above result then reduces to 

(GjvM)oB = Gjvrn^(l-h^(afff)o + (c^)oo + |soo 

+ ^(aJff)o-i(c%o). (230) 

This expression for (GjvM)ob contains a linear combina- 
tion of coefficients for Lorentz violation that is indepen- 
dent of the three combinations {GnM)'yy), {Gj^M)y)s, 
and {GnM)-rs obtained for photon tests. Some com- 
ments about tests with this result are provided in the 
next subsection. 



E. Experiments 

The above subsections show that each type of photon 
test of Lorentz symmetry is sensitive to an effective value 
of GnM that contains a combination of coefficients for 
Lorentz violation. The time-delay value {GnM)^!) is 
given by Eq. (207), and it depends on the coefficients 
Q:(agg)o, (c^)oo, and sqo- The gravitational Doppler shift 
involves the value (GArM)Ds in Eq. (213) and involves 
the same combination of the three coefficients. The value 
(GjvAf)Rs for the gravitational redshift is given by Eq. 
(220), which contains a different combination of coeffi- 
cients and varies also with ^ciock- All three of these pho- 
ton tests yield sensitivities differing from those in orbital 
tests, which involve the value (GjvM)ob in Eq. (230) 
that depends also on the test-body coefficients Q;(a^)o 
and {0^)00 ■ Note that no qualitatively new signals are in- 
volved in any of these cases, since the effects are merely 
scalings of established physics. In contrast, the nuU- 
rcdshift observable given in Eq. (227), which depends 
on the difference of clock quantities AS,ab, represents 
a qualitative departure from conventional gravitational 
physics. 

Comparisons of the time-delay value {GnM)td or 
the Dopplcr-shift value (GjvAf)DS to the redshift value 
(GjvM)rs for the same source body can be used to ob- 
tain sensitivity to combinations of coefficients for Lorentz 
violation. High-quality data for the time delay and the 
gravitational Doppler shift have been obtained by track- 
ing the Cassini spacecraft [133] in the gravitational field 
of the Sun. Proposed missions such as the Astrodynam- 
ical Space Test of Relativity using Optical Devices (AS- 
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TROD) [134], the Mercury Orbiter Radio-scienee Exper- 
iment (MORE) [135], the Search for Anomalous Gravita- 
tion using Atomic Sensors (SAGAS) [136], and the Solar 
System Odyssey (SSO) [137] have the potential to im- 
prove these measurements using the Sun as the gravita- 
tional source, while the Beyond Einstein Advanced Co- 
herent Optical Network (BEACON) [138] could sharpen 
results using the Earth as the gravitational source. An- 
other relevant recent proposal involves the use of very- 
long-baseline interferometry (VLBI) [139] to measure the 
deflection of radio waves from distant sources by solar- 
system objects. The sensitivity of this measurement to 
Lorentz violation is likely to be comparatively weaker but 
may be offset by the enhanced access to independent co- 
efBcient combinations offered by multiple measurements 
and perhaps by access to anisotropic effects involving spa- 
tial components of s^,^. 

Redshift tests permit sensitivities to effects controlled 
by ^ciock- These can be isolated either by comparing sep- 
arate redshift tests performed with different clocks in 
the same gravitational source or more directly by null- 
redshift tests, in which the signal depends on the differ- 
ence A^^s between two clocks A, B and vanishes in the 
absence of gravity. The results of some investigations of 
local position invariance can be reinterpreted as measure- 
ments of A^^B. For example, a recent Earth-based test 
comparing a hydrogen maser with a Cs fountain [140] 
obtained a sensitivity that corresponds to the bound 

ICh-^CsI < (0.1±1.4) X 10-^ (231) 

while another comparing a hydrogen maser with a cryo- 
genic sapphire oscillator [141] yields the measurement 

- -ecso = (-2.7 ± 1.4) X 10-^ (232) 

These resiilts offer a benchmark for currently attainable 
sensitivities to ^ciock. The two experiments involve dif- 
ferent clocks and hence likely different sensitivities to the 
coefficients {c^)^v Calculating the specific constraints 
on (c™)^iy and possibly other coefficients for Lorentz vi- 
olation from these and other tests is an interesting open 
project. Note that Earth-based searches of this type typ- 
ically take advantage of the annual and diurnal variations 
in the gravitational potential of the Sun as experienced 
in the laboratory. In searches using the annual variation, 
it is challenging and perhaps impossible to disentangle 
gravitational effects of nonzero ^ciock from other Lorentz- 
violating effects in Minkowski spacetime. However, di- 
urnal searches can distinguish the two types of effects 
because the Minkowski-spacetime signals occur at the 
sidereal frcqu(;n(;y instead. Note also that other clock- 
comparison tests normally viewed as sensitive to SME 
coefficients in Minkowski spacetime may also have sensi- 
tivity to ^ciock- One intriguing possibility is that suitable 
choices of clocks could separate effects from (c*)tt and 
(c^)tT) which would then lead to independent sensitivi- 
ties to (aefj)T and (agjf)T, a result otherwise challenging 
to achieve. 



Satellites carrying two different clocks offer interesting 
prospects for improved null-rcdshift searches for Lorentz 
violation. Since the attainable sensitivities improve with 
the gravitational potential difference according to Eq. 
(227), it is desirable to acquire elliptical orbits. The 
Space-Time Asymmetry Research (STAR) program [142] 
presently under development proposes to compare two 
different clocks on a satellite traveling in an elliptical or- 
bit. This mission could improve sensitivities to ^ciock by 
an order of magnitude or more relative to ground-based 
tests. Improved sensitivities may also be possible by com- 
paring clocks aboard the proposed SAGAS spacecraft. 
Note also that experiments in highly elliptical orbits can 
be expected to have increased sensitivity to anisotropic 
effects on the redshift produced by sj/c . 

Provided effects due to ^ciock are excluded, either 
through independent experiments or by using a clock 
with ^ciock = 0, then the dependence of {Gm^)ks on 
(a^ff)T implies that measurements of the gravitational 
redshift can be compared with other photon tests per- 
formed with the same gravity source to obtain indepen- 
dent sensitivities to (aJJ^ )T. The Gravity Probe A (GPA) 
mission [143], which used the Earth as the gravity source, 
confirmed the conventional gravitational redshift to parts 
in 10^. This result could eventually be combined with 
proposed time-delay or Doppler-shift measurements of 
the BEACON type to yield sensitivity to the coefficient 
(agfj)T for the Earth. Improved tests of the gravitational 
redshift are also proposed for the Atomic Clock Ensem- 
ble in Space (ACES) [144], SAGAS, and STAR missions. 
With the Sun as the gravity source instead, the Galileo 
space probe obtained sensitivity to deviations for the 
gravitational redshift at the level of parts in 10^ [145]. 
Given knowledge of ^ciock for the Galileo clock, this re- 
sult could be combined with the Cassini results to yield 
sensitivity to the coefficient (afff)T for the Sun. If feasi- 
ble, a redshift test performed directly with Cassini would 
be of interest in this respect. Other gravitational sources 
could also be used. For example, the gravitational red- 
shift was measured to parts in 10^ using Saturn as the 
source during the flyby of Voyager [146]. Time-delay or 
Doppler-shift data could therefore permit sensitivity to 
the coefficient (acfr)T for Saturn. 

Combinations of photon tests with measurements of 
the effective orbital mass are also of interest. In the 
limit of zero matter-sector Lorentz violation, the result 
(230) for {GnM)ob has been combined with Eq. (207) 
for (Gjv.^)TD to extract sensitivity to stt [131]. How- 
ever, with nonzero coefficients (aeff);^ and c^,y, the effec- 
tive value (G'ArM)oB involves properties of the test body 
as well as the source. Note that these appear in the famil- 
iar combination a(aJg)T — m'^(c"^)TT/3 discussed in Sees. 
VII and IX B. The WEP tests considered in this work 
constrain the degree to which this combination can differ 
between neutrons and neutral combinations of electrons 
and protons, but only through the indirect arguments 
involving binding energy described in Sec. VI B. In con- 
trast, comparing GjvM factors for measurements with 
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orbiting bodies and for photon tests offers the opportu- 
nity to obtain direct sensitivity to Q;(agg)T — to'^(c")tt/3 
and a{alg^)T — to'^(c®+^)tt/3. Comparisons with sen- 
sitive gravimetcrs may also be of interest in this respect. 

We remark in passing that proposed sensitive experi- 
ments to measure gravitational light bending, including 
the Laser Astrometric Test of Relativity (LATOR) [147] 
and the Space Interfcrometry Mission (SIM) [148], are 
likely to have signals affected by Lorcntz violation. The 
attainable sensitivities can be expected to be similar to 
those discussed above, but the analysis of this possibility 
lies beyond our present scope. 



XII. SUMMARY 

This work studies the gravitational couplings of matter 
in the presence of Lorcntz violation. The framework for 
the investigation is the fermion sector of the gravitation- 
ally coupled minimal SME in a post-newtonian expan- 
sion. Our primary goal is to develop a suitable method- 
ology for searches for Lorentz and CPT violation that 
exploit the couplings of matter to gravity, incorporating 
in particular effects that are challenging or impossible to 
detect in Minkowski spacetime. 

Section II presents the basic formalism for the work. 
The action for the gravity-matter system is given in Sec. 
II A, and the linearization procedure is outlined in Sec. 
II B. Some types of Lorentz violation are unobservable 
in principle. This issue is discussed in Sec. II C, which 
also fixes the coordinate choice (23) used in this work. 
The metric and coefficient fields for Lorcntz violation can 
fluctuate about their background values, and the corre- 
sponding interactions must be incorporated in analyses 
of experiments. In Sec. II D, we develop general per- 
turbative techniques to analyze these fluctuations. Two 
notions of perturbative order are introduced. One is de- 
noted 0(m, n) and tracks the orders in Lorentz violation 
and in gravity, while the other is denoted PNO(p) and 
tracks the post-newtonian order. The goal of this work is 
to investigate dominant terms involving Lorentz violation 
in gravity, which are at 0(1,1). 

Section III studies the quantum theory of the gravity- 
matter system. Starting from the field-theoretic action, 
we construct the relativistic quantum mechanics in the 
presence of gravitational fluctuations and Lorentz viola- 
tion. Formulating the quantum theory for matter in the 
presence of gravitational fluctuations is a standard chal- 
lenge. In Sec. Ill A, we present a solution to this problem 
via a fleld redefinition, which yields a hamiltonian that 
is hermitian with respect to the usual scalar product for 
wave functions. We then use this procedure in Sec. Ill B 
to extract the explicit form of the relativistic hamiltonian 
involving all coefficients for Lorcntz violation in the min- 
imal QED extension. The result forms the appropriate 
starting point for general investigations of Lorentz and 
CPT violation in matter-gravity couplings. To maintain 
a reasonable scope in this work, we subsequently special- 



ize our focus to the study of spin-independent Lorentz- 
violating effects, which arc governed by the coefficient 
fields (flcff)^, c^i^ and the metric ffuctuation h^^. The 
nonrelativistic quantum hamiltonian for this case is ob- 
tained in Sec. Ill C using the standard Foldy-Wouthuysen 
procedure. 

Measurements of gravity-matter couplings typically 
are performed at the classical level. Section IV con- 
structs the classical theory associated with the quantum- 
mechanical dynamics of matter involving nonzero (oeff)^, 
c^^, and /i^y. The behavior of test and source bodies in 
the presence of Lorcntz violation is the subject of Sec. 
IV A. Working from the action for a point particle, we 
provide expressions for the mass and for the effective co- 
efficients for Lorcntz violation for a test or source body, 
along with the effective action (76) describing the dy- 
namics of the body. These results enable the derivation in 
Sec. IV B of the modified Einstein equation and the equa- 
tion (78) for the trajectory of a test particle. To apply 
this equation in practice requires knowledge of the coeffi- 
cient and metric fluctuations. In Sec. IV C, we develop a 
systematic methodology for calculating this information 
in perturbation theory and obtain general expressions for 
the coefficient and metric fiuctuations to 0(1,1) in terms 
of various gravitational potentials and the background 
coefficient values (aoff)/^ and c^j^. 

To illustrate the application of the general formalism, 
we consider in Sec. V a specific class of bumblebee mod- 
els, which are theories with a vector field driving spon- 
taneous Lorentz breaking. The action for the bumble- 
bee field Bfj^ is given in Sec. VA, where a match at the 
field-theoretic level to the general formalism of earlier 
sections is made and the coefficient fields (ccff)^ and c^;/ 
are identified in terms of and the metric. In Sec. VB, 
we explicitly solve the model at the relevant order in per- 
turbation theory, extract the modified Einstein equation, 
and derive the equation for the trajectory of a test parti- 
cle. The results are shown to match those obtained using 
the general formalism developed in the earlier sections. 

The largest portion of the paper is devoted to a dis- 
cussion of experiments and observations that can achieve 
sensitivity to the coefficients {aeff)i^ and c^^. Section 
VI presents some general material broadly applicable to 
searches for Lorentz violation. Various choices of refer- 
ence frame and their relationship to the canonical Sun- 
centered frame are discussed in Sec. VIA. Attainable 
sensitivities to the coefficients (aes)^ and c^y in any mea- 
surement procedure are constrained by certain generic 
features. Section VI B considers somc^ of thc^sc. includ- 
ing the role of binding energy in impeding or aiding the 
analysis of WEP tests for signals of Lorentz violation. 

A major class of searches for Lorentz violation involves 
laboratory tests with ordinary neutral bulk matter, neu- 
tral atoms, and neutrons. Section VII treats this topic. 
The PN0(3) lagrangian describing the dynamics of a test 
body moving near the surface of the Earth in the presence 
of Lorentz violation is considered in Sec. VII A. Expres- 
sions are given in an Earth-centered frame and the trans- 
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formation to the laboratory frame is outlined. The re- 
sulting description of laboratory signals for gravitational 
Lorentz violation includes effects from the matter-sector 
coefRcients (Seff)^ and Cj^j/ and ones from the gravity- 
sector coefficients s^j, obtained in Ref. [7]. It reveals 
that the gravitational force acquires tiny corrections both 
along and perpendicular to the usual free-fall trajectory 
near the surface of the Earth, while the effective iner- 
tial mass of a test body becomes a direction-dependent 
quantity. These effects can be sought in numerous lab- 
oratory experiments. Since the standard relationship 
between force and acceleration is modified, it is useful 
to distinguish tests measuring gravitational acceleration 



from ones comparing forces. In Sec. VII B, we consider 
free-fall gravimeter tests such as falling corner cubes and 
atom interferometry. Force-comparison gravimeter tests 
using equipment such as superconducting gravimeters are 
studied in Sec. VII C. An important potential signal for 
gravitational Lorentz violation arises from the flavor de- 
pendence of the effects, which implies signals in WEP 
tests. A variety of free-fall WEP tests is considered in 
Sec. VII D, while force-comparison WEP tests with a tor- 
sion pendulum are treated in Sec. VII E. For all the tests 
considered, the possible signals for Lorentz violation are 
decomposed according to their time dependence, and es- 
timates of the attainable sensitivities are obtained. 



Table XIV. Summary of actual and attainable sensitivities in past or present tests. 
Coefficient Gravimeter Free-fall Force-comparison Solar 

combinations WEP WEP system 
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i)m"(c")Q 




10-8 GeV 




/ — e+p — \ 


imJ'(r+P-")TT 
















10-" GeV 



[10-6 GeV] 
[10-6 GeV] 



(«e'ff)^ + («:ff)^ + 0.1(a^ff)^ ... ... ... 10-6 GeVt 

(a^a.)e + (a:ff)e + 0.1«ff)e ... ... ... 10-6 GcVt 

(c")(Tj) [10-7] [10-4] [10-7] [10-6] 

(c")q ... 10-8t 10-8t 

(c")^ •■• 10-«t 

(c")e ... ... ... 10-7t 



Section VIII considers satellite-based WEP tests, 
which offer interesting prospects for improved sensitiv- 
ities to Lorentz violation. In this context, the signal for 
Lorentz violation is an anomalous time variation of the 
relative local acceleration between two test bodies of dif- 
fering composition located on the satellite. We derive 
the frequency decomposition of the signal for Lorentz vi- 



olation, and we consider idealized scenarios for several 
proposed satellite-based WEP tests. Based on the de- 
sign reach of the missions, we estimate the sensitivities 
that could be achieved to various combinations of the 
matter-sector coeflacients (aefi)/^ and Cp,v 

Studies of the gravitational couplings of charged par- 
ticles, antimatter, and second- and third-generation par- 
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tides present distinct experimental challenges but can 
yield sensitivities to Lorentz and CPT violation that are 
otherwise difficult or impossible to achieve. Section IX 
addresses some of these possibilities, including charged- 
particle interferometry, ballistic tests with charged par- 
ticles, gravitational experiments with antihydrogen, and 
signals in muonium free fall. For antihydrogen experi- 
ments, simple toy models are introduced to illustrate as- 
pects of their discovery potential and to address attempts 
to place indirect limits on possible effects. 

Traditional tests of gravity couplings to matter include 
observations of the motion of bodies within the solar 
system. Section X contains a discussion of the signals 
accessible via lunar and satellite laser ranging and via 
measurements of the precession of the perihelion of or- 
biting bodies. A reanalysis of existing data from lunar 
laser ranging could yield interesting sensitivities to some 
combinations of the matter-sector coefficients (aeflf)^ and 
c^i,. We use the established advance of the perihelion 
for Mercury and for the Earth to obtain constraints on 
combinations of {aes)fi, C/j.^, and 



The interaction of photons with gravity offers a differ- 
ent arena in which to seek Lorentz and CPT violation. 
Section XI is devoted to this topic. We consider signals 
arising in measurements of the photon time delay, studies 
of the gravitational Doppler and redshifts, and compar- 
isons of the behaviors of photons and massive bodies. A 
variety of existing and proposed experiments on space- 
craft offer interesting prospects for these measurements. 

Tables XIV and XV collect estimated sensitivities to 
the matter-sector coefficients and {c^)ij,„ obtained 

from many of the measurements discussed in this work. 
These tables disregard possible effects from the pure- 
gravity coefficients s^^ that could in principle be relevant 
to solar-system tests. Table XIV concerns existing data, 
while Table XV tabulates future prospects. One result 
omitted from these tables is the generalization (121) of 
the constraint obtained in Ref. [11] using data from force- 
comparison WEP tests with a torsion pendulum. In Sec. 
VII E, multiple data sets are combined to separate this 
constraint into the two limits (153), and both of these 
are included in Table XIV instead. 



Table XV. Summary of attainable sensitivities in future tests. 

CoefScient Free-fall Free-fall Satellite Solar 

combinations gravimeter WEP WEP system 



(^iKs^)x 10-10 GeV 10-10 GeV IQ-" GcV 

aial+P)Y+z 10-10 GeV lO-^o QeV IQ-" GcV 

a{al^P)Y 10-« GeV IQ-^ GeV 10"^ GeV 

aial+P)z 10-« GeV IQ-^ GeV IQ-^ GeV 

a{al+P)T ... 10-" GeVt 10-15 GeVt 

a(a^ft)x 10-10 GeV lO'io GeV IQ-n GcV 

a{a^g)Y+z 10-10 GeV IQ-io GeV IQ-n GeV 

a{a^e)Y IQ-^ GeV IQ-^ GeV IQ-^ GeV 

a{a^g)z 10-^ GeV IQ-^ GeV 10-^ GeV 

a(a^ff)T . . . 10-1" GeVf lO-i^ GeV^ 

a{al+''-")x 10-10 GeV IQ-" GeV lO'i^ GeV 10"^ GeV 

a(a^+P-")y+z 10-10 GeV IQ-n GeV IQ-i^ GeV 10-^ GeV 

a(a^^^'-")F 10-8 GcV lO^^ GcV lO-io GeV 

a{al+''-")z 10-8 GeV lO-o GeV IQ-io GeV 

- lmP{c''+P-'')TT ■ ■ ■ 10-15 GeV IQ-i^ GeV 

(c")(Tj) 10-10 10-11 10-12 10-7 

(c")q ... 10-15 10-16 



The formalism and the analytical results for gravita- 
tional signals of Lorentz violation presented in this work 
apply to the nonzero matter-sector coefficients (a™g)yj, 
(c^)/ii^ and in some cases also to the gravity-sector coef- 
ficients Sf^u . Comparatively little is known about the co- 
efficients (ajfg:)^, and scenarios exist in which they could 
be countershaded, having large values while still escaping 
notice in searches to date [11]. However, nongravitational 



measurements have already yielded impressive sensitivi- 
ties to various components of (c^)ju,y [2]. The estimated 
attainable sensitivities to (c"')^^ derived in this work are 
therefore primarily restricted to components of (c")^,^, 
for which existing constraints are weaker. Tables XIV 
and XV reflect these facts, containing mostly entries for 
combinations of the coefficients (a^)/i along with some 
results for (c*")^,^. 
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Tabic XIV summarizes actual sensitivities or estimated 
attainable ones using data from past or present measure- 
ments. The table is based on the ealeulations presented 
in this work and includes only sensitivities below parts 
in 10^. Each entry in the first column of this table rep- 
resents a linear combination of coefficients that is acces- 
sible in principle via existing searches. Each of the other 
four columns contains our estimates for sensitivities that 
could be achieved in the listed class of tests, expressed 
to the nearest order of magnitude. Values in these four 
columns that are shown without brackets represent order- 
of-magnitude sensitivities implied by our present analysis 
to the modulus of the coefficient combination displayed. 
Values appearing in brackets in the table represent our es- 
timate of sensitivities that could in principle be obtained 
from a suitable reanalysis of existing data. An obelisk 
(f) following a value indicates a limit attainable under 
the assumption that either {a^g)^ or {c^)nu is negligibly 
small or vanishes. A diesis (J) indicates a sensitivity that 
is attained by combining data from two different classes 
of experiments, and this sensitivity is placed in each of 
the two corresponding columns in the table. 

Table XV contains future attainable sensitivities to the 
moduli of various combinations of the matter-sector co- 
efficients and (c^)^i/, as estimated in previous sec- 
tions of this work. The structure of this table is similar 
to that of Table XIV. The listed entries are based on the 
best design reach and are given to the nearest order of 
magnitude. For each class of search, we assume enough 
measurements have been performed to achieve the maxi- 
mum number of independent sensitivities. The reader is 
cautioned that for certain coefficients a single measure- 



ment cannot attain the indicated sensitivity, but instead 
only a linear combination of coefficients with multipliers 
controlled by composition and orientation factors. Note 
that elsewhere in this work the convention is to display 
values of future sensitivities in braces, but this conven- 
tion is suppressed in Table XV because all entries are of 
this type. Note also that further improvements in the- 
oretical techniques and experimental design in all types 
of searches, including ones not listed in Table XV such 
as exotic gravitational tests or photon tests, are expected 
to yield additional interesting prospects for future attain- 
able sensitivities. 

Taken together, Tables XIV and XV reveal excellent 
prospects for using matter-gravity couplings to seek ef- 
fects of Lorentz violation. The opportunities for mea- 
suring the countershaded coefficients {a^fi)ij, at sensitive 
levels are of particular interest in this context, as these 
coefficients typically cannot be detected in nongravita- 
tional searches. Indeed, the spatial components of {a^g)n 
remain essentially unconstrained to date. The tests pro- 
posed here can be performed with existing or near-future 
technology, and they offer a promising new arena for 
searches for signals from the Planck scale. 
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